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SECTION  I 


INTRODUCTION 

Hvperqolic  liquid  rocket  propellants  are  toxic,  and  accidental  releases 
durinq  handlinq  and  storaqe  may  endanqer  the  civilian  population.  The  U.S. 

Air  Force  developed  procedures  to  predict  toxic  hazard  corridors  that  would 
result  from  the  release  of  any  individual  propellant  material.  These 
corridors  result  from  the  evaporation  and  subsequent  dispersion  of  the  vapors. 

A  different  situation  arises  when  fuel  and  oxidizer  come  into  contact  as 
the  result  of  an  accidental  release.  In  this  case,  a  violent  explosion  can 
occur  that  will  send  the  combustion  products  and,  possibly,  some  of  the 
unreacted  Drooellants  hiqh  into  the  air.  When  this  occurs,  the  atmospheric 
dispersion  beqlns  at  the  niqh  elevation  and  cannot  be  determined  uslnq  the 
same  calculation  method  as  that  for  a  qround  release. 

Phase  I  of  this  Droqram  (Source  Characterization),  addressed  the  reactions 
between  nitroqen  tetroxide  and  hydrazines.  The  reaction  products  and  the  heat 
released  by  the  reaction  were  quantified  durinq  that  work.  That  information 
was  used  to  determine  the  combustion  time  and  the  helqht  the  resultinq 
fireball  achieves  before  it  cools  and  passively  disperses  with  the  air. 

Phase  II  of  the  Dro.iect  (Atmospheric  Dispersion  .odelinq),  which  is 
described  in  this  report,  investiqated  the  chemical  and  physical  interactions 
of  the  combustion  products  with  air  and  the  disoersion  of  these  products  In 
the  environment. 

The  chemical  interaction  of  reaction  products  of  nitroqen  tetroxide 
combined  with  a  hydrazine  was  determined  by  evaluatinq  the  literature  on  the 
topic  and  by  laboratory  experimental  work.  This  information  showed  which 
products,  both  toxic  and  nont.oxic,  could  be  expected  to  form  and  dissipate  as 
the  cloud  moves  away  from  the  accident  location.  Furthermore,  the  chemical 
kinetics  of  these  reactions  were  determined  so  the  rate  of  appearance  and 
disappearance  of  each  constituent  could  be  assessed. 
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The  atmospheric  dispersion  of  the  products  was  described  by  considerinq 
the  coupled  affects  of  convection,  dispersion,  and  chemical  reaction.  A 
mathematical  description  of  this  transport  process,  a  Hvperqolic  Accidental 
Release  Model  (HARM),  was  developed  by  H.F.  Cramer  Co.  durinq  this  project. 
They  combined  the  source  characterization  alqorithm  developed  at  Martin 
Marietta  (Prince  198?  ,  1 983 )  ,  the  houyant  cloud  models  by  Rriqqs  (1970), 

and  the  dispersion  model  for  rocket  launches  by  Rjorklund  (1982)  to  describe 
this  phenomenon. 

This  report  summarizes  two  interim  reports:  "Atmospheric  nispersion  of 
Hvoerqolic  Liquid  Rocket  Fuels,  Phase  II:  Atmospheric  Dispersion  Modelinq"by 
Stephen  Prince  and  "User's  Manual  for  the  Hyperqolic  Accidental  Release  Model 
(HARMl  Computer  °roqram"  by  C.R.  Rowman,  W.R.  Harqraves,  .l.R.  Rjorklund,  R.K. 

Dumbauld,  and  J.F.  Rafferty.  Complete  copies  of  these  reports  are  included 
as  Appendices  A  and  B  respectively. 


SECTION  II 

INTERACTIONS,  REACTIONS  AND  ANALYSES 


Pl.t  IMF /ATMOSPHERE  INTERACTIONS 

The  hyperqolic  rocket  prooellants  and  the  products  of  combustion  formed 
after  an  accidental  release  of  these  materials  can  be  toxic.  The  results  of 
work  done  during  this  portion  of  the  project  provide  chemical  reaction  rate 
data  for  the  interaction  of  hydrazine  (N?H4),  unsymmetrical 
dimethyl  hydrazine  (IIOMH),  nitroqen  tetroxide  (NTO),  and  combustion  products 
with  air  and  tho  reactive  components  in  the  air  which  include  ozone,  oxyqen, 
water  vapor,  and  carbon  dioxide.  These  reaction  rate  data,  when  coupled  with 
the  dispersion  characteristics  of  the  effluent  plume,  can  be  used  to  assess 
the  environmental  impact  of  the  vapors  and  the  control  necessary  to  avoid 
injuries  to  personnel  in  the  area  of  the  hyperqolic  release.  A  summary  of  the 
important  reactions  that  can  take  place  in  this  situation  will  be  qiven  here. 

A  more  thorouqh  discussion  of  the  chemistry  is  presented  in  the  Phase  I  report 
and  a  more  extensive  discussion  of  the  chemical  reaction  rates  Is  presented  in 
the  Interim  Deport:  Task  1  "Dlume/Atmosoheric  Interaction"  included  as 
Appendix  A. 

N?ha  REACTIONS 

At  high  N^  concentrations  in  air,  the  major  reaction  pathway  for  its 
depletion  is  with  molecular  oxvqen, 

+  Opr* — ^*.N  7  +  ?HpO. 

As  the  concentration  decreases  to  the  part  per  million  level  the  major 
reactions  are  with  the  atmospheric  pollution  contaminants;  hydroxide  radical 

(oh),  ozone  ^3),  and  nitroqen  dioxide  (NO?).  These  reactions  are  about 
an  order  of  magnitude  faster  than  the  ones  with  oxvqen  and  Drevail  as  lonq  as 
the  rate  of  ohotox idation  in  the  atmosphere  is  great  enouqh  to  replenish  the 
reactants. 

The  products  of  these  reactions,  in  qeneral,  have  low  toxicity.  These 
products  are  N? ,  h?0,  NH3 , 
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IIHMH  is  reactive  with  the  same  atmospheric  components  as  is  N7H4.  At 
hiqh  HHMH  cone .ntration,  the  reaction  occurs  with  molecular  oxyqen  and  forms 
formaldehyde  dimethvlhydrazone  (FDH)  as  a  major  product, 

HJHMH  +  ?0?  — ^.?FHH  +  N?  +  4H?0. 

At  lower  MHMH  concentrations  when  the  atmosoheric  contaminants'  concentrations 
are  aooroachinq  those  of  the  hydrazines,  they  aqain  are  the  major  contributor 
to  the  decomposition.  The  products  of  these  reactions  include 
nitrosodiinethylamine  (NHMA,  a  known  carcinooen)  and  other  substituted 
tetrazines  and  nitramine  products  that  are  toxic.  The  known  products  from  air 
oxidization  of  IIHMH  are  formaldehyde  dimethylhydrazone  (FHH),  ammonia, 
dimethylamine,  NHMA,  di azomethane,  nitrous  oxide,  methane,  carhon  dioxide, 
formaldehyde,  formaldehyde  monomethyl  hydrazone  (FHM),  and 
tetramethyltetrazine  (TMT). 

n itpo^fn  ninxTHF  RfArrinNS 

Nitroqen  dioxide  (NO?)  reacts  in  two  major  ways  in  the  atmosphere. 
Atmospheric  water  vaoor  combines  with  NO?  to  produce  nitric  acid  (HNO3) 
vapor  and  at  the  reaction  site  the  concentration  of  these  materials  may  exceed 
the  vaoor  pressure  of  the  water  solution  and  form  liquid  nucleation  sites 
which  will  lead  to  a  mist.  Nitroqen  dioxide  is  also  photochemical ly  active 
and  enters  into  the  photolysis  reaction  train  that  produces  ozone  (O3). 

This  series  of  chemical  reactions  results  in  other  oxidized  nitroqen  species 
includinq  NO,  NO?,  no?,  and  N?0^  and  upon  interaction  with  other 
atmospheric  contaminants,  photochemical  smoq  is  produced. 

RFAFTION  FXPFR1MFNTS 
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Supplemental  chemical  rate  data  were  ohtained  at  Martin  Marietta  on 
reactions  between  the  reactive  components  In  air  (oxvqen,  carhon  dioxide,  and 
water  vapor)  and  vaporized  propellant  soecies.  These  data  were  used,  alonq 
with  published  information, to  develop  rate  expressions  compatible 
with  the  dispersion  model inq  theory.  The  propel  1  ant- atmospheric  reaction 
oxoerime  tal  test  condition  matrix  usnd  for  this  work  is  shown  in  Table  1. 
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TARLF  1  .  TFST  MATRIX  ATMOSPHFRTC  RFACTION  OF  HYPERROLS 


Test 

No. 

Propellant 

Snecies 

°rop. 

Cone. 

Reactant  Ras* 
Composition 

Sampl inq 

Time,  Hrs 

1 

Hydrazine 

0.5* 

Hel ium 

0,0. 5,1, 7, 4 

P 

Hydrazine 

0.5* 

PI*  0? 

in  He 

0,0.25 

1 

Hydrazine 

0.5* 

PI*  0? 

+  1*  H?0  in 

He 

0,0.25 

4 

Hydrazine 

0.5* 

21*  n? 

+  0.05*  CO? 

in  He 

0,0.25 

5 

Hydrazine 

0.5* 

PI*  0? 

+  0.05*  C0? 

+  1*  H?0  in 

He 

0,0.25 

6 

1 JHMH 

0.5* 

Helium 

0,5,12,21,45,70 

7 

liMOH 

0.5* 

71*  0? 

in  He 

0,5,  IP, 21, 45, 70 

8 

IIOMH 

0.5* 

21*  0? 

+  1*  H?0  in 

He 

0,5,12,21 ,45,70 

q 

IIOMH 

0.5* 

21*  0? 

+  0.05*  C0? 

in  He 

0,5,12,21,45,70 

10 

IIOMH 

0.5* 

21*  0? 

+  0.05*  C0? 

+  1*  H?0  in 

He 

0,5,12,21,45,70 

n 

N0? 

0.P5* 

Hel ium 

0,1,3,7,24 

ip 

N0? 

0.P5* 

PI*  o2 

in  He 

0,1,3,7,24 

13 

N0? 

0.P5* 

PI*  0? 

+  1*  H20  in 

He 

0,1,3,7,24 

14 

N0? 

0.P5* 

21*  0? 

+  0.05*  C0? 

in  He 

0,1,3,7,24 

*  1*  H^o  in  Helium  Corresponds  to  50  Percent  Relative  Humid  Air  at  25°C 

(77°f) 


Test  temoerature  were  ambient  (?1-?3°C)  for  amine  fuel  tests  and 
elevated  (35°C)  for  nitrogen  dioxide  tests. 


ratf  stiioy  smmmary 


A  summary  of  the  chemical  rate  constant  information  assembled  durinq  this 
study  is  presented  in  Table  7.  This  information  was  used  to  formulate  the 
rate  of  chemical  consumotion  in  the  HARM  model  for  accidental  releases  of 
propel lant. 

TARLF  ?  .  ATMOSPHERIC  REACTION  RATE  CONSTANTS  FOR 

firfrall  oispfrsion  mooflino 


Fireball 

Species 

Concentration 

in 

Atmosphere 

Meteoroloqical 

Condition 

Rate 

Constant 

Half- 

Life 

Major 

Reaction 

Products 

Percent 

Yield 

Hydrazine 

0-1  ppm 

Ory  or  Humid 

Air 

0.?  pom*1 

7  Hr* 

HpOp> 

NA 

Min* 1 

Hydrazine 

Above  1  oom 

Dry  Air 

O.OS4  Hr'1 

10.8  Hr 

N£ ,HpO 

100 

Hydrazine  Above  1  Dnm 

Humid  Air 

0.141  Hr*1 

4.0  Hr 

Np  ,H20 

100 

NOMH 

0-1  ppm 

Ory  or  Humid 

Air 

0.15  ppm'1 

0.7  Hr* 

NOMA 

50 

Min* 1 

UOMH 

Above  1  npm 

Ory  Air 

0.0007  Hr'1 

PRO  Hr 

FOH** 

57 

MOMH 

Above  1  Dnm 

Humid  Air 

0.015  Hr*1 

45  Hr 

FPH** 

57 

N0? 

10  com 

Humid  Air 

f.5X1041? 

5  Hr 

HN07(q) 

100 

mole'^sec" 1 

N0? 

SO  opm 

Humid  Air 

5.5X10419 

0.?  Hr 

HNO^(q) 

100 

mole*?sec* 1 

N0? 

Above  SO  opm 

Humid  Air 

5.5X1041? 

0.7  Hr 

HN03(q) 

nr, 

mole-?  sec* 1 

N0o 

All  Cone. 

Ory  Air 

NO  RFACTION 

*  Clean,  Unpolluted  Atmosphere 

**  NfiMA  Formed  in  0.’  °ercent  Yield  (.00?  Moles  NOMA  Found  Per  1  Mole  UDMH 
Reacted) 

NO  =  Not  Available 


OVERALL  ENERGY  CONSIDERATIONS 


The  combustion  products  contained  in  the  "fireball"  that  leave  a 
hyoerqolic  accident  are  cooled  by  air  entrainment  and  heat  transfer  to  the 
environment.  The  assumption  was  made  that  the  maior  cooling  affect  after  the 
effluent  cloud  has  stabilized  is  caused  uy  air  entrainment.  On  cooling,  the 
temperature  cloud  can  fall  below  the  saturation  point  of  some  of  the 
constituents  at  hich  concentrations. 

The  calculations  necessary  to  determine  if  there  is  a  possibility  for 
condensation  are  shown  in  the  reoort  in  Appendix  A.  The  results  of  these 
calculations  revealed  that  even  when  there  is  excess  fuel  (a  situation  that 
creates  the  most  likelv  chance  for  condensation  to  occur  because  the 
concentration  of  the  low  vapor  pressure  material  is  hiqhest  and  the 
temperature  is  lowest),  no  possibility  for  condensation  exists. 

The  overall  results  obtained  durinq  this  part  of  the  project  furnished  the 
data  needed  to  incorporate  the  effects  of  chemical  reactions  in  the  environment 
into  the  overall  dispersion  model  for  propellant  releases. 


SECTION  III 

HYPERGOl.IC  ACCIDENTAL  RELEASE  MODEL  (HARM) 


The  MvDerqol  ic  Accidental  Release  Model  (HARM)  was  developed  t.o  predict 
the  transport  and  dispersion  of  the  products  downwind  of  a  hyperqolic  reaction 
between  fuel  and  oxidizer.  This  model  is  an  improvement  over  the  existinq 
method  used  to  assess  the  hazard  of  propellant  releases  (Ocean  Rreeze/Ory 
Gulch  Model)  because  it  can  compute  results  for  elevated  sources,  thermally 
buoyant  clouds,  and  considers  the  effect  of  chemical  reactions  takinq  place  in 
the  cloud.  The  HARM  model  can  predict  qround  level  concentration  of  toxic 
materials  that  result  from  a  catastrophic  reaction  of  hvoerqolic  propellants. 

A  summary  of  the  model  is  oresented  here  and  Appendix  II  is  a  cooy  of  the 
liser^_s  Manual  for  the  Hvoerqolic  Accident  Release  Model  (HARM)  Computer 
Program  by  C.R.  Rnwman,  W.R.  Harqraves,  T.R.  Rjorklund,  R.K.  Oumbauld,  and 
1 . F  .  Rafferty. 


Vhe  HARM  model  is  currently  constructed  to  evaluate  the  environmental 
impact  that  would  result  from  the  accidental  release  and  subsequent  reaction 
of  the  oropel lants  in  a  Titan  II  TCRM  vehicle.  The  thermal  ererqy  released 
from  the  reaction  of  A-RO  and  NTO  heats  the  combustion  products  and  unreacted 
propellants  to  form  a  buoyant  cloud.  The  model  uses  basic  i  .formation  about 
the  cloud  and  its  environment  to  determine  the  concentration  and  dosage  of  the 
constituents  at  around  level  as  a  function  of  time  and  distance  downwind  from 
the  reaction  site.  Rased  on  information  from  Pnase  I  of  this  project.  Source 
Characteri zati on,  and  specific  information  about  the  accident  (quantities  of 
propellant  spilled  and  location  of  the  accident),  the  proqram  calculates  the 
temperature,  buoyancy,  and  composition  of  the  fireball  at  lift-off. 

An  instantaneous  cloud  rise  model  that  has  been  adapted  to  rocket  launches 
is  incoroorated  into  HARM  to  determine  the  stablization  height  and  rise  time 
for  the  buoyant  cloud  that  has  resulted  from  the  combustion.  The  results  of 
these  calculations  determine  the  initial  conditions  for  the  next  phase  of  the 
r.oinouta t i on .  These  conditions  are  very  important  because  they  fix  the  region 
in  the  atmosphere  wher°  disnersion  ho. jins  and,  based  on  the  local  conditions, 
wi’l  det=»r>n  i  no  the  transport  process. 

id 


The  dispersion  of  the  stabilized  cloud  is  determined  by  dividinq  the 
atmosphere  into  sections  that  comply  with  rawinsonde  measurements  made  at 
prescribed  times  to  qlve  the  wind  direction  and  soeed  and  temperature  profile 
for  the  atmosohere.  The  dispersion  model  is  derived  assuminq  that  the  source 
of  material  to  be  transported  can  be  represented  by  a  *1nite  verlcal  line 
source  in  each  of  the  rawlnsonde-def ined  layers.  The  diameter  of  the  cloud  at 
the  stabilized  heiqht  Is  equal  to  the  diameter  used  In  the  cloud-rise  model. 
The  concentration  of  source  material  In  each  layer  is  uniform  in  the  vertical 
direction  and  Gaussian  alonq  the  wind.  The  assumption  is  also  made  that  the 
concentration  at  one  cloud  radius  from  the  center  line  of  the  cloud  is  0.1 
times  the  maximum  concentration.  The  concentration  variation  in  the  vertical 
direction  is  uniform  in  each  lave*"  but  Gauss  ion,  layer  by  layer.over  the  height 
.if  the  cloud. 

The  mixinq  and  expansion  of  the  cloud, as  the  material  moves  downwind  are 
calculated,  based  on  the  standard  deviation  of  the  vertical  and  horizontal 
wind  direction  over  a  10-minute  time  Interval.  Standard,  tested  correlations 
are  available  to  determine  the  deqree  of  mixinq  that  results  and  are  given  in 
the  renort  in  Appendix  B. 

The  computer  proaram  that  calculates  the  concentration  and  dosaqe  (time- 
inteqrated  concentration)  uses  a  rectanqular  coordinate  system  with  the  origin 
at  qround  level  under  the  cloud  center  at  the  time  of  stabilization  with  the  x 
axis  directed  alonq  the  mean  wind  velocity  vector  and  the  y  axis  orthoqonal  to 
the  wind  vector.  All  points  on  the  qround  are  then  related  hack  to  the  point 
of  the  accident,  and  expressed  in  a  polar  coordinate  system  tbrouqh  a 
coordinate  transformation. 

The  calculation  routine  in  the  model  also  accounts  for  the  consumotion  of 
chemical  soecies  and  the  generation  of  others  usinq  the  reaction  rate  data 
presented  earlier.  Standard  quantitive  pais...... .  v.  '.Iteration  of  the 

constituents  by  precipitation  are  also  included. 

Section  S  of  the  HARM  reoort  (Apoendix  B)  qives  complete  execution  and 
data  oreparation  instructions  for  the  HARM  computer  ornqram.  It  also 
describes  the  format  available  to  plot  the  meteoroloqical  data  and  the  output 
concentration  and  dosaqe  data.  Format  information  for  the  input  rawinsonde 
data  is  also  described. 
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1.0  INTRODUCTION 


This  document  and  the  experimental  work  described  herein  is  intended  to 
provide  time  dependent  gas  phase  reaction  rate  data  between  hypergolic  rocket 
propellants  (hydrazine,  unsymmetrical  dimethylhydrazine  or  UDMH  and  nitrogen 
tetroxide  or  NTO)  with  the  reactive  atmospheric  components  of  air  (which 
include  ozone,  oxygen,  water  vapor,  and  carbon  dioxide).  This  reaction  rate 
(kinetic)  data  will  be  used  to  modify  an  existing  atmospheric  plume  dispersion 
model  to  incorporate  time-dependent  atmospheric  reactions  into  the  dispersion 
pattern  of  an  instant  release  source  of  a  hypergolic  fireball.  A  computer 
program  will  be  developed  to  model  this  instant  release/continuously  reacting 
hypergolic  plume  which  uses  atmospheric  meteorological  data  and  predicts  the 
ground  level  concentrations  of  toxic  materials  that  result  from  a  catastrophic 
reaction  of  hypergolic  propellants.  Additional  information  on  the 
identification  of  the  gas-phase  reaction  products  between  the  hypergolic 
propellants  and  air  were  obtained  during  this  task  period.  The  purpose  of 
defining  atmospheric  reaction  rates  is  to  assess  the  potential  environmental 
impact  of  vaporized  propellants  and  the  corresponding  degree  of  control 
required  on  their  releases  into  the  atmospheric  following  a  '.ypergolic 
fireball  explosion. 
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2.0  COMPILATION  OF  EXISTING  ATMOSPHERIC  REACTION  DATA 


2.1  Hydrazine-Air  Reaction 


Since  hydrazines  do  not  photolyze  in  the  actinic  region  (A*- 290 

/> 

nmp  the  major  atmospheric  sinks  for  hydrazine  decomposition  include 
hydroxide  radical,  ozone,  and  nitrogen  dioxide,  all  of  which  are  produced 
in  the  N0X  photooxidation  system  in  polluted  and  urban  atmospheres,  and 
by  molecular  oxygen  in  clean  unpolluted  air.  In  general,  the  decay  of 
hydrazine  by  OH,  0^,  or  NO2  Is  an  order  of  magnitude  more  rapid  than 
the  decay  by  molecular  oxygen  and  these  processes  are  thought  to  be  more 
prevalent  when  the  concentration  of  hydrazine  in  the  atmosphere  is  at  or 
near  the  part-per-mil lion  (ppm)  level.  At  higher  hydrazine 
concentrations  it  Is  believed  that  the  concentration  of  these  reactive 
photooxidation  products  is  rapidly  depleted  by  hydrazine  reaction  and  the 
reaction  of  hydrazines  by  molecular  oxygen  becomes  the  major  atmospheric 
reaction  pathway. 

The  nominal  concentrations  of  these  reactive  species  in  the  clean 
lower  troposphere  is  1  x  10^  molecule  cm"^  hydroxide  radical,  40  ppb 
ozone1,  and  between  10  and  50  ppb  N0X  (NO2  and  NO).  Polluted 
industrial  atmospheres,  however,  will  exhibit  substantailly  higher 
concentrations  of  these  atmospheric  reactants.  Low-level  hydrazine 
concentrations  are  very  reactive  in  these  polluted  urban  atmospheric 
systems. 

The  reaction  between  hydrazine  and  hydroxide  radical  is  second-order 
according  to  Equation  (1). 

0H  +  N2H4  — ~h2°  +  n2h3  (D 

^1  =  6.1  +  1.0  x  10"11  cm^  molecule"1  sec"1 
(9.0  x  10^  ppm"1  min"1) 
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and  the  rate  constant  kj  i$  relatively  Independent  of  temperature 
between  300  and  425°K.3  Assuming  an  average  hydroxide  concentration 
of  1  x  10^  cm"3  in  the  lower  troposphere,  the  half-life  for  hydrazine 
in  the  ppm  concentration  range  can  vary  between  less  than  one  hour  In 
polluted  atmospheric  environments  up  to  six  hours  in  clean  environments. 

The  reaction  between  hydrazine  and  atmospheric  ozone  is  also  second 
order  and  its  reported  rate  constant  is  temperature  dependent.  Tuazon 
et.  al.  report  a  rate  constant  of  1.3  +  0.1  x  10"16  cm3  molecule”1 
sec-1  at  24°C  (0.2  ppm"1  min*1)  with  an  activation  energy  of  3.5 
to  5  kilocalories  per  mole.  Hydrogen  peroxide  (H2O2)  was  the  major 
product  In  this  oxidation  reaction,  and  nitrous  oxide  (^0)  was  a  minor 
product.  The  remaining  gas-phase  decomposition  products  were  postulated 
to  be  nitrogen,  water  vapor,  and  hydrogen  gas.  The  tropospheric 
half-life  for  ppm  concentration  of  hydrazine  vapor  were  estimated  by 
Tuazon  to  be  about  10  minutes  during  ozone  pollution  episodes  to  less 
than  two  hours  in  an  unpolluted  atmosphere. 

The  reaction  between  ppm  levels  of  hydrazine  vapors  and  atmosphere 
N0X  concentrations  also  at  ppm  levels  produced  no  observable  change  in 
the  initial  N0X  reactant  when  performed  in  the  dark,  but  proceeded 
rapidly  under  photolytic  conditions.^  In  the  later  case,  hydrazine  is 
rapidly  consumed,  NO  is  converted  to  N02  at  a  rate  approximately  equal 
to  that  of  the  hydrazine  decay,  and  03  is  formed  following  NO 
consumption.  The  only  nitrogen-containing  product  observed  was  N2O  in 
trace  amounts.  When  a  1.3  ppm  hydrazine  sample  was  sunlight  irradiated 
with  0.9  ppm  NO  and  0.12  ppm  NO2  the  hydrazine  was  completely  depleted 
in  30  minutes.  Later  studies  by  Tuazon  et.  al.  indicated  that  none  of 
the  hydrazines  studied,  which  included  hydrazine,  MMH,  and  UDMH,  reacted 
with  NO  alone  at  measureaole  rates,  but  did  react  significantly  with 
NO2.  The  rate  constant  for  the  dark  reaction  between  NO2  and 
hydrazine  reported  by  Tuazon  was  approximately  2.5  x  lO’1^  cm3 
molecule"1  sec"1  (3.7  x  10"^  ppm"1  min"1).  The  reaction  of 
NOj,  with  N2H4  resulted  in  the  formation  of  high  yields  of  nitrous 
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acid  (HONO),  hydrazinium  nitrate,  diazine  (in  excess  hydrazine  only),  and 
traces  of  N20  and  NH-j.  Results  of  these  reported  experiments 
indicate  that  hydrazine  undergoes  reactions  at  significant  rates  with 
hydroxide  radical,  ozone,  and  nitrogen  dioxide  when  present  at  ppm  levels 
in  air,  and  Pitts  et.  al.2  has  suggested  that  under  most  conditions, 
reaction  with  ozone  will  be  the  major  fate  of  ppm  levels  of  hydrazines 
released  into  the  atmosphere. 

The  reaction  between  hydrazine  and  molecular  oxygen  is  believed  to 
be  the  major  atmospheric  reaction  pathway  when  fuel  concentrations 
greatly  exceed  the  part-per-million  level.  The  reaction  between 
hydrazine  vapor  and  molecular  oxygen  is  somewhat  complex,  and  proceeds  at 
a  rate  which  is  strongly  dependent  upon  reaction  cell  geometry  and 
surface  composition. 

Stone^  has  reported  that  the  main  oxidation  reaction  of  hydrazine 
by  molecular  oxygen  is  expressed  by  equation  2: 

N2H4  +  02  ~>N?  +  2H20  (2) 

In  addition  to  the  main  reaction  which  produces  nitrogen  and  water  vapor, 
there  are  side  reactions  which  produce  ammonia  and  are  largely 
heterogeneous  in  nature.  The  oxidation  half-lives  studied  by  Stone 
depended  on  both  the  surface  area  and  composition  of  the  reaction  vessel 
and  ranged  from  25  minutes  for  a  five  liter  flask  (surface  area  =  1600 
cm^)  two  ancj  a  ha^f  hours  for  a  specially  constructed  44  x  2  -cm 
cell  (surface  area  =  300  cm^).  The  decay  rate  for  each  of  these 
reactions,  however,  was  determined  to  be  first  order  with  respect  to 
initial  hydrazine  concentration.  The  rate  of  ammonia  production,  and 
thus  the  total  amount  produced  was  also  a  function  of  the 
surface-to-volume  ratio  of  the  reaction  vessel.  The  5-liter  flask  had 
produced  approximately  0.2  torr  ammonia  (.03%  v/v)  from  an  initial 
hydrazine  concentrations  of  5.0  torr  (.66%  v/v)  in  an  80%  nitrogen  -20% 
oxygen  gas  mixture  in  which  the  total  pressure  was  760  torr.  The 
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half-life  for  ammonia  production  appeared  to  be  approximately  40 
minutes.  The  44  x  2  -  cm  cell  had  produced  approximately  0.3  torr 
ammonia  (.04%  v/v)  from  an  initial  hydrazine  concentration  of  3.5  torr 
(.46%  v/v)  in  an  80%  nitrogen  -  20%  oxygen  gas  mixture  in  which  the  total 
pressure  was  760  torr.  In  this  latter  case,  the  half-life  for  ammonia 
production  was  approximately  160  minutes.  The  surface-to-volume  ratio 
for  the  5-liter  flask  and  44  x  2  -  cm  cell  was  0.3  cm"*  and  2  cm“^, 
respectively. 

Tuazon  et.  al.  measured  the  dark  decay  of  hydrazine  at  low 
concentration  (5-12  ppm)  in  a  purified  air  environment,  i.e.  ozone  and 
concentration  were  negligible.  The  study  was  conducted  in  a  3800 
liter  chamber  fabricated  from  50  micron  (2-ml)  thick  Teflon  sheeting. 

The  calculated  surface-to-volume  ratio  of  this  chamber  was  ca.  0.06 
cm"*,  a  factor  of  five  lower  than  the  reaction  vessels  used  by  Stone. 

The  reported  half-life  for  hydrazine  decay  at  22°C  and  12%  relative 
humidity  was  10.8  hours,  and  at  55%  relative  humidity  was  4.9  hours. 
Ammonia  was  also  observed  as  a  reaction  product  in  these  tests,  but  only 
account  for  5-10%  of  the  hydrazine  lost  during  reaction. 

The  decay  of  hydrazine  in  air  therefore  was  seen  to  be  highly 
dependent  on  the  geometry  of  the  gas  reaction  vessel  as  well  as  the 
amount  of  humidity  in  the  air,  showing  an  approximately  two-fold  faster 
oxidation  rate  in  humid  air  than  in  dry  air.  Actual  atmospheric 
oxidation  rates,  in  which  a  typical  ambient  surface-to-volume  ratio  of  9 
x  10"®  cm"l  is  taken,®  is  unknown  at  present  but  will  most  likely 
approximate  those  values  performed  in  the  larger  reaction  chambers. 

2.2  Unsymmetrical  Dimethylhydrazine-Air  Reaction 

In  general,  the  same  atmospheric  components  which  were  found  to  be 
reactive  with  hydrazine  vapors  are  also  reactive  with  substituted 
hydrazine  vapors,  which  include  both  unsymmetrical  dimethyl  hydrazine 
(UDMH)  and  monomethyl  hydrazine  (MMH).  The  substituted  hydrazines  did 
exhibit  different  reaction  rates  with  the  reactive  atmospheric  species 
and  the  reaction  products  were  more  diverse  in  nature. 
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The  rate  constant  for  the  reaction  between  hydroxide  radical  and 
UDMH  was  estimated  by  Tuazon  et.  al.*  to  be  approximately  5  x  10"** 
cm-*  molecule"*  sec"*  (7.4  x  104  ppm"*  min"*).  The  rate 
constant  was  assumed  to  be  second  order  and  relatively  temperature 
independent  between  300  and  425°K.  The  half-life  for  UDMH  decay  by 
atmospheric  hydroxide  radical  was  estimated  to  be  approximately  3  hours, 
hut  could  be  two  to  five  times  larger  in  pristine  atmospheres  and  less 
than  one  hour  in  polluted  urban  atmospheres. 

When  UDMH  vapors  were  reacted  with  ozone  in  ppm  concentration 
levels,  the  fuel  reacted  at  a  rate  which  was  too  rapid  to  accurately 
measure.  The  half-life  for  substituted  hydrazines  due  to  ozone  oxidation 
was  estimated  by  Tuazon^  to  be  a  factor  of  10  shorter  than  that  of 
unsubstituted  hydrazine.  This  would  place  the  window  for  ppm  UDMH  decay 
hy  atmospheric  ozone  between  1  minute  during  ozone  pollution  episodes  to 
.2  hr  (12  minutes)  in  clean  atmospheres.  The  overall  apparent  rate 

constant  for  UDMH  +  O3  reaction  therefore  must  be  greater  than  10’*^ 
cm^  molecule-*  sec-*  (1.5  ppm“l  min-*).  In  these  reactions, 
Nitrosodimethyl amine  (NDMA,  a  known  carcinogen)  was  formed  in  high  yields 
and  the  observed  O^/UOMH  stoichiometry  was  1.5/1.  Tuazon*  also 
performed  studies  on  atmospheric  destruction  mechanisms  for  the  NDMA 
formed  from  the  UDMH/ozone  reaction.  He  determined  that  the  major 
atmospheric  sink  for  NOMA  destruction  is  the  photolysis  by  ambient 
sunlight.  The  reaction  desplayed  a  quantum  yield  of  1.0,  and  no 
appreciable  decay  of  NDMA  by  air  or  ozone  was  observed  under  atmospheric 
conditions.  The  photolysis  reaction  has  a  half-life  of  ca.  30  minutes 
and  results  in  the  formation  of  dimethylnitramlne  as  well  as  other  minor 
products. 

The  reaction  between  UDMH  and  N02  produced  tetramethy Itetrazine 
(TMT)  as  the  major  organic  reaction  product  and  displayed  an  apparent 
second-order  rate  constant  of  2  x  10-*^  cm^  molecule-*  sec-*  (.03 
ppm-*  mir'l),  with  a  half-life  of  about  5  minutes.  Nitrous  acid 
(HONO)  was  also  a  product  in  the  UDMH/NO^  reaction.  In  these  three 
cases  for  low-level  UDMH  decay  by  atmospheric  species,  the  oxidation  by 
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ozone  was  assumed  to  be  the  major  fate  of  UDMH  vapors  released  into  the 
atmosphere.  The  main  reaction  product  from  this  reaction  is  NDMA  which 
displays  unacceptable  personnel  exposure  hazards  at  any  measurable 
concentration. 

The  reaction  between  UDMH  and  molecular  oxygen  is  believed  to  be  the 
major  environmental  fate  for  elevated  fuel  concentrations  in  the 
atmosphere.  The  vapor-phase  reaction  between  UDMH  and  molecular  oxygen 
is  very  complex  and  may  involve  free-radical  processes."^  Free  radical 
reactions  are  greatly  influenced  by  the  geometry  and  surface  composition 
of  the  reaction  vessel,  as  rate  determining  intermediates  are  pressumed 
to  occur  on  the  surface  of  the  reaction  vessel.  Additionally,  different 
reaction  mechanisms  (and  thus  reaction  rates)  appear  to  be  prevalent  at 
higher  fuel  concentrations  (0. 2-1.0  percent)  than  at  lower  fuel 
concentrations  (2-10  ppm).4  Stone4  reported  an  atmospheric  half-life 
of  approximately  175  hours  for  an  initial  UDMH  concentration  of  0.7%, 
which  was  identical  for  reaction  with  dry  nitrogen,  20%  oxygen  in  dry 
nitrogen,  and  80%  oxygen  in  dry  nitrogen.  Tuazon  reported  a  half-life 
for  12-13  ppm  UDMH  in  "dry"  air  of  341jf384  hours  in  a  6400  liter 
environmental  chamber  and  341+64  hours  in  a  3800  liter  chamber.  When  the 
air  was  humidified  to  50%  R.H.,  this  half-life  was  decreased  to  70.9_^9.6 
hours.  These  results  substantiate  the  effect  of  vessel  geometry  and 
perhaps  reactant  fuel  concentration  on  the  determination  of  kinetic  rate 
constants  and  half  lives  for  UDMH  decay  by  atmospheric  air.  Perhaps  more 
importantly,  they  indicate  that  the  reaction  between  UDMH  and  dry  air  is 
relatively  insignificant  under  the  conditions  studied,  and  the  driving 
force  in  the  air  oxidation  of  UDMH  becomes  the  degree  of  humidity. 

The  reaction  between  UDMH  and  molecular  oxygen  is  first  order  (zero 
order  with  respect  to  oxygen),  and  thus  the  reaction  rate  is  directly 
proportional  to  the  UDMH  concentration.  The  main  reaction  as  proposed  by 
loper^  proceeds  according  to  Equation  3. 

3UDMH  +  202  c>2FDH  +  4Hj,0  +  H?  (3) 
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Thus  the  major  product  in  the  air  oxidation  of  UDMH  is  formaldehyde 
dimethyl hydrazone  (FDH)  and  it  is  produced  in  a  reaction  stoichiometry  of 
1  mole  FOH  formed  per  1.5  mole  UDMH  reacted.  Other  minor  products  in  the 
UDMH  air  reaction  include  ammonia  (NH^),  dimethyl  amine,  MDMA, 
di azomethane,  nitrous  oxide,  methane,  carbon  dioxide,  formaldehyde, 
formaldehyde  monomethyl  hydrazone  (FMH)  and  TMT. 

2.3  Nitrogen  Dioxide-Air  Reactions 

Nitrogen  dioxide  ( NO2 )  reacts  with  atmospheric  air  by  two  main 
pathways:  (1)  it  interacts  with  sunlight  to  produce  photochemical  smog 
which  is  representative  of  polluted  urban  atmosphere  such  as  Los 
Angeles®;  and  (2)  it  reacts  with  atmospheric  water  vapor  to  produce 
both  vapor  phase  and  condensed-phase  nitric  acid.  These  atmospheric 
reactions  of  nitrogen  dioxide  and  the  implications  for  dispersion 
modeling  of  released  NO^  vapors  are  discussed  in  the  sections  that 
follow. 

2.3,1  Photochemical  Reactions  of  Nitrogen  Dioxide. 

The  photolysis  of  nitrogen  dioxide  by  ambient  sunlight  is  the  major 
reaction  resulting  in  the  formation  of  ozone  in  the  troposphere®**®. 

The  ozone  produced  by  the  photolysis  reaction  can  react  with  nitric  oxide 
to  reform  the  nitrogen  dioxide,  with  nitrogen  dioxide  to  form  symmetrical 
nitrogen  trioxide  (NO^),  or  with  hydrocarbon  pollutants  in  the 
atmosphere  to  produce  oxidized  carbon  species. 

In  general,  photochemical  reactions  producing  tropospheric  (or 
stratospheric )  ozone  as  well  as  the  transient  nitrogen  trioxide  species 
and  dinitrogen  pentoxide  ( ^2^5 )  *ias  been  extensively  studied  by 
atmospheric  and  meteorological  scientists**’*^’*®. 

The  major  photolysis  reactions  of  N0j>  can  be  expressed  by  the 
following  series  of  equations. 
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0(3P)  +  02  4.  M  ->  O3  +  M 


(5) 


or  Oq  +  NO - *N0p  +  Op 

(6) 

^3  +  NO?  — NO3  +  0,> 

(7) 

M3  +  NO?  +  M - ->  Np05  +  M 

(8) 

In  equations  (4)  -  (6),  nitrogen  dioxide  is  regenerated  by  the 
action  of  ozone  on  the  nitric  oxide  formed  by  reaction  (4).  In  equation 
(7)  and  (8)  which  represent  an  alternate  reaction  pathway  for  the  ozone 
formed,  dinitrogen  pentoxide  (^O^)  is  formed  through  the  transitory 
and  unstable  nitrogen  trioxide  molecule  (NO3).  The  equilibrium 
composition  of  the  final  mixture  and  the  kinetics  are  complex  since  the 
reaction  mechanisms  are  competitive  for  the  reactiv.  v/;one  molecule.  The 
reaction  rates  for  each  of  the  individual  reactions  have  been  elucidated 
and  are  described  as  follows. 

The  photolysis  of  nitrogen  dioxide  (Eq.  4)  is  the  preliminary  step 
in  the  production  of  photochemical  smog.  This  reaction  provides  both 
nitric  oxide  and  photochemical ly  excited  triplet  oxygen  atom  (03P). 

This  process  was  observed  to  have  a  quantum  yield  of  1.0  (one  excited 
oxygen  atom  produced  per  photon  absorbed)  when  irradiated  with  sunlight 
of  wavelength  d95  -  398  nm^.  The  first-order  rate  constant  was 
determined  to  be  0.43  min’*  under  ambient  conditions®.  The  actual 
ohotolytic  kinetic  rate  constant,  k^t  however,  was  determined  by 
Sickles  et.  al.3^  to  be  highly  dependent  on  the  latitude,  solar  zenith 
angle,  and  cloud  cover.  Other  factors  influencing  the  transparency  of 
the  atmosphere  to  sunlight  irradiation,  including  aerosol  loading  and 
particulate  count,  were  also  found  to  preclude  the  accurate  prediction  of 

k4  during  the  sunlight  portion  of  the  day.  Sickles  did  observe  an  apex 
of  the  of  the  rate  constant  value  during  solar  noon  (.46  for  a  sunny  -  0% 
sky  cover  day)  which  was  reduced  to  ca.  .76  (clear  day)  and  .14  (cloudy 
day)  at  early  morning  hours  and  late  afternoon  hours,  respectively. 


The  highly  reactive  triplet  P-oxygen  formed  in  air  collides  with 
oxygen  molecules  forming  ozone  by  Equation  (5).  In  this  case,  M 
represents  a  3rd  molecular  species  (nitrogen,  oxygen,  water  vapor,  etc) 
that  removes  a  fraction  of  the  energy  released  during  the  interaction  of 
the  oxygen  atom  with  the  oxygen  molecule.  This  energy  removal  stabilizes 
the  ozone  product.  Typically,  reaction  (5)  is  very  rapid  and  Garvin1^ 
reports  a  third  order  rate  constant  of  5.8  x  10*^  cm^  molecule-^ 
sec"1  at  25°C  and  for  M  =  N2.  A  steady  state  concentration  of 
ozone  is  quickly  reached  by  the  reaction  pathways  (6)  and  (7).  Reaction 
(6)  is  very  rapid  as  is  reaction  (7).  Graham  and  Johnston^  report  a 
second-order  rate  constant  for  reaction  (7)  of  3.4  x  lO"1^  cm^ 
molecule’1  sec-1  at  25°C  (5  x  10"2  ppm-1  min"1).  Reaction 
(8)  which  combines  the  transient  nitrogen  trioxide  with  nitrogen  dioxide 
to  form  dinitrogen  pentoxide  is  also  very  rapid. 

\ 

Nitrogen  dioxide  destruction  by  photolysis  is  a  complex  phenomenon 
and  the  subsequent  reformation  by  reaction  wit/i  nitric  oxide  and  ozone 
can  lead  to  a  net  decrease  in  the  atmospheric  nitrogen  dioxide  equal  to 
the  steady-state  ozone  level  developed.  The  calculation  of  steady-state 
ozone  production  (and  nitrogen  dioxide  decay)  for  an  atmospheric  release 
of  hyper go  1 1c  fireball  components  is  difficult  at  best,  and  a  qualitative 
description  of  these  photochemical  reactions  of  nitrogen  dioxide  will 
suffice  for  the  present  task  effort. 

2.3.2  Gas-Phase  Reaction  Between  Nitrogen  Dioxide  and  Water  Vapor 

The  kinetics  of  the  gas-phase  reaction  between  nitrogen  dioxide  and 
water  vapor  was  studied  by  England  and  Corcoran1^.  Under  normal 
atmospheric  conditions  (in  which  N02  gas  is  present  at  ppm  levels  in 
polluted  urban  atmospheres),  the  reaction  with  water  vapor  was  very  slow 
and  formation  of  nitric  acid  aerosols  did  not  form.  This  observation  has 
been  contrasted  to  the  atmospheric  nucleation,  hydration,  and  misting  of 
sulfuric  acid  from  the  reaction  of  sulfur  dioxide  with  moist  air®.  The 
apparent  difference  in  the  acid  rain  effect  for  sulfur  dioxide  and  that 
for  nitrogen  dioxide  was  accounted  for  by  the  high  vapor  pressure  of 


10 


27 


nitric  acid  compared  to  the  low  vapor  pressure  (10~®  to  10-^  torr) 
of  sulfuric  acid.  The  vapor  pressure  of  pure  nitric  acid  was  51  torr 
(7.2%  at  one  atmosphere)  at  20°C.  The  vapor  pressure  of  a  50%  nitric 
acid-water  solution  was  ca.  470  ppm  nitric  acid  and  1.02%  water  vapor 
also  at  20°C.  The  gas-phase  reaction  between  nitrogen  dioxide  and 
water  vapor  supplies  nitric  acid  nuclei  when  partial  pressure  of  acid 
exceeds  the  vapor  pressure.  Water  vapor  readily  condenses  on  each 
nucleus  to  form  strong  nitric  acid  droplets.  England  and  Corcoran 
suggest  that  the  threshold  for  the  two-phase  region  of  nitric  acid  occurs 
at  about  50  ppm  HNO^  at  one  atmosphere  pressure  and  25°C.  Above  this 
concentration,  mists  or  condensation  of  nitric  acid  were  observed  to 
form.  The  overall  reaction 

3N°2(g)  +  H20(g) - *N0(g)  +  2HN03(aq)  (9) 

is  favored  thermodynamically  at  ambient  conditions. 

When  gas-phase  reactions  were  performed  at  concentrations  low  enough 
to  prevent  the  formation  of  a  two-phase  system,  the  rate  of 
disappearance  of  NO2  was  first-order  with  respect  to  nitrogen 
dioxide.  The  overall  reaction  in  the  presence  of  oxygen 

4no2  +  2H20  +  02 - >  4HN03  (10) 

goes  to  completion  in  the  gas  phase  and  the  third  order  rate  constant  Is 
represented  by  equation  (11) 

RN02  *  -k  (N02)2  (h20)  (11) 

k  =  5.5  x  104  1?  mole-?  sec"l  at  250C. 

This  rate  constant  indicates  that  the  initial  rate  of  reaction  increases 
as  the  square  of  the  nitrogen  dioxide  concentrations  and  therefore  would 
he  very  slow  at  low  N02  concentrations  and  very  fast  at  higher 
concentrations. 
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For  a  IX  concentration  of  water  vapor  and  an  initial  concentration 
of  10  ppm  of  nitrogen  dioxide,  the  half-life  for  reaction  would  be  about 
5  hr.  If  the  initial  concentration  of  nitrogen  dioxide  was  50  ppm  (the 
maximum  level  of  NO2  before  misting  or  condensation  occurs);  the 
reaction  half-life  would  be  12  minutes. 


Therefore,  in  cases  where  nitrogen  dioxide  vapors  are  released  into 
the  atmosphere  from  a  hypergolic  explosion,  misting  and  condensation  is 
expected  to  form  where  the  concentration  is  elevated  above  50  ppm  NO2, 
at  which  time  the  residual  NO2  vapor  will  quickly  react  to  form 
vapor-phase  nitric  acid. 
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3.0  PROPELLANT-ATMOSPHERE  REACTION  EXPERIMENTS 


Chemical  reactions  between  the  reactive  molecular  components  of  air 
(oxygen,  carbon  dioxide,  and  water  vapor)  and  vaporized  propellant  species 
(hydrazine,  UDMH,  and  nitrogen  dioxide)  were  performed  In  order  to  supplement 
the  reaction  rate  data  obtained  In  current  literature  sources  and  to  provide 
Information  on  the  gas-phase  reaction  products  between  vaporized  hypergolic 
propellants  and  air.  Results  of  these  experiments  will  be  correlated  with 
previously  reported  atmospheric  reaction  data  and  incorporated  into  an 
existing  atmospheric  plume  model  to  incorporate  time-dependent  atmospheric 
reactions  into  the  dispersion  pattern  of  an  instant  release  source  of  a 
hypergolic  fireball.  The  test  matrix  for  the  propellant-atmospheric  reaction 
experiments  is  included  in  Table  1. 

3.1  Experimental  Methods 

3.1.1  Hydrazine-Atmospheric  Reactions 

Gas  phase  mixtures  of  0.5%  hydrazine  vapor  were  prepared  by  injecting  a 
known  amount  of  anhydrous  hydrazine  liquid  (5.2  microliters)  into  a  one  liter 
gas  sampling  bulb  (Ace  Glass  Inc.)  which  was  previously  purged  with  pure 
helium,  a  21%  oxygen  helium  mixture,  or  a  21%  oxygen  .05%  carbon  dioxide  in 
helium  mixture  as  detailed  in  Table  I.  Atmospheric  gases  were  obtained  In  the 
ultra-pure  grade  from  Linde  Specialty  Gases  and  mixed  using  appropriate 
flowraters  to  give  the  desired  oxygen  and  carbon  dioxide  composition.  Water 
vapor  was  added  to  the  sampling  bulbs  Nos.  8  and  10  by  injection  of  6.0 
microliters  of  liquid  distilled  water.  This  corresponds  to  1%  water  vapor 
under  Denver  barometric  conditions  (25°C,  0.806  atmospheres).  Vaporization 
of  the  liquid  components  and  gas  mixing  was  accomplished  using  l/4-1nch  teflon 
beads.  Gas  samples  were  removed  from  the  bulbs  at  the  time  Increments 
detailed  In  Table  I  and  analyzed  for  hydrazine  by  the 
paradimethylamlnobenzaldehyde  (PDAB)  colorimetric  method  and  for  gaseous 
reaction  products  (nitrogen,  oxygen,  water  vapor,  and  ammonia)  using  gas 
chromatographic  techniques.  The  PDAB  colorimetric  analysis  method  for 
hydrazine  is  summarized  below. 
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Table  I  Test  Matrix  Atmospheric  Reaction  of  Hypergols 

Test 

Propellant 

Prop. 

Reactant  Gas* 

Sampling 

No. 

Spec  1 es 

Cone. 

Composition 

Time,  Hrs 

1 

Hydrazine 

0.5% 

Helium 

0,0. 5, 1,2, 4 

2 

Hydrazine 

0.5% 

21%  02  in  He 

0,0.25 

3 

Hydrazine 

0.5% 

21%  02  +  XX  h20  In 

He 

0,0.25 

4 

Hydrazine 

0.5% 

21%  02  +  0.05%  C02 

in  He 

0,0.25 

5 

Hydrazine 

0.5% 

21*  02  +  0.05%  C02 

+  1%  H20  in  He 

0,0.25 

6 

UOMH 

0.5% 

Hellian 

0,5,12,21,45,70 

7 

UDMH 

0.5% 

21*  02  In  He 

0,5,12,21,45.70 

8 

UOMH 

0.5% 

21%  02  +  1%  h20  in 

He 

0,5,12,21,45,70 

9 

UDMH 

0.5% 

21*  02  +  o.05<  C02 

In  He 

0,5,12,21,45,70 

10 

UDMH 

0.5% 

21*  02  +  0.05*  C02 

+  1%  H20  In  He 

0,5,12,21,45,70 

11 

no2 

0.25% 

Helium 

0,1,3,7,24 

12 

no2 

0.25% 

21*  02  In  He 

0,1,3,7,24 

13 

no2 

0.25% 

21*  02  +  1%  H20  in 

He 

0,1,3,7,24 

14 

no2 

0.25% 

21*  02  +  0.05%  C02 

In  He 

0,1,3,7,24 

*  1 %  H2o  in  Helium  Corresponds  to  50%  Relative  Humid  Air  at  25°C 


(77°F) 


Test  temperatures  were  ambient  (21-23°C)  for  amine  fuel  tests  and 
elevated  (35°C)  for  nitrogen  dioxide  tests. 


Analysis  of  Hydrazine  Vapor  Composition  by  the 
Paradimethylaminobenzaldehyde  Method  (PDAB) 


Reagents 

1.  PDAB  Stock  Solution  -  Dissolve  10.0  grams  PDAB  In  approximately  150  ml 
methyl  alcohol.  Add  10  ml  concentrated  hydrochloric  acid  and  dilute  to 
250  ml  with  methyl  alcohol. 

2.  0.1  N  Sulfuric  Acid  -  Dilute  2.78  ml  concentrated  H2SO4  to  1000  ml 
with  distilled  water. 

3.  PDAB  Final  Solution  -  Dilute  100  ml  of  the  PDAB  stock  solution  to  500  ml 
with  0.1  N  H2so4. 

i 

1 

Hydrazine  Calibration  Curve  1 

A  hydrazine  calibration  curve  was  prepared  by  injecting  krcw.  amounts  of 
hydrazine  stock  solution  Into  the  PDAB  stock  solution,  allowing  liie  >Jllow 
color  to  fully  develop  for  10  minutes,  and  reading  the  absorbance  of  the 
resulting  solution  at  a  wavelength  of  455  nm  In  a  Spectronlc  20  Colorimeter. 

The  resulting  calibration  curve  for  anhydrous  hydrazine  Is  depicted  in 
Figure  1. 

Sample  Analysis 

1.  Pipet  10.0  ml  PDAB  final  solution  into  a  stoppered  30  ml  sample  vial. 


Inject  1.0  ml  of  the  gas  mixture  contained  In  the  gas  sampling  bulb 
through  the  septum  In  the  30  ml  sample  vial  using  a  gas  tight  syringe. 
The  gas  sample  to  be  analyzed  Is  added  to  the  gas  ullage  above  the  PDAB 
solution  contained  in  the  vial. 
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Analysis  of  Hydrazine  Vapor  Composition  by  the 
Paradimethylamlnobenzaldehyde  Method  (PDAB) 
(Continued) 


3.  Shake  the  vial  vigorously  for  30  seconds  to  ensure  complete  hydrazine 
reaction  by  the  POAB  solution. 

4.  After  10  minutes  color  development  time,  read  the  absorbance  of  the 
solution  at  455  nm. 

5.  Compare  the  absorbance  of  the  analyzed  sample  with  the  absorbance  of  the 
calibration  standards  and  report  the  hydrazine  vapor  concentration  In  the 
sample. 

Interference  Studies 


Because  ammonia  was  observed  by  some  researchers  to  be  a  side-product  in 
the  hydrazine-air  reaction, 5  the  PDAB  colorimetric  reagent  was  tested  for 
Its  reaction  and  color  development  with  ammonia  gas.  The  Interference  ratio 
for  ammonia  with  the  PDAB  reagent  was  approximately  100:1  (100  ppm  NH-j 
corresponded  to  1  ppm  N2H4). 
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3.1.2  UOMH-Atmospheric  Reactions 


Gas  mixtures  containing  0.5*  UDMH  vapor  (Test  Nos.  6-10  in  Table  I)  were 
prepared  by  Injecting  12.6  microliter  of  high  purity  1,  1-dlmethylhydrazine 
liquid  (  0.2*  FOH  by  volume)  Into  a  1-liter  gas  sampling  bulb  (Ace  Glass 

Inc.)  which  was  previously  purged  with  the  desired  reactant  gases  as  described 
in  Section  3.1.1.  To  the  UDMH- atmospheric  mixtures  contained  In  the  gas 
sampling  bulb,  1.0  ml  of  pure  ethane  gas  (Scott  Calibration  Gases)  was  added 
as  an  Internal  standard. 

Water  vapor  was  added  to  the  appropriate  sample  bulbs  by  injection  of  6.0 
microliters  of  liquid  distilled  water  and  vaporization  and  mixing  of  the 
liquid  component  was  accomplished  using  the  l/4-1nch  teflon  beads  as 
previously  described.  Aliquots  of  the  gas  reaction  mixture  were  removed  from 
sample  bulb  at  the  reaction  time  intervals  indicated  In  Table  1,  and  the  bulbs 
were  stored  In  the  dark  at  ambient  temperatures  (21-23°C)  between  sampling 
intervals.  Oxygen,  water  vapor  and  UDMH  reactant  gas  composition  as  well  as 
FDH,  NDMA,  and  nitrogen  reaction  products  were  quantitated  using  the  gas 
chromatographic  techniques  to  be  subsequently  described.  The  gas 
chromatographic  analysis  of  hydrazine-air  mixture  Is  also  included  in  this 
section. 

Preparation  of  Standards  for  GC  Calibration 

Gas  calibration  standards  were  prepared  In  1-liter  gas  sampling  bulbs 
using  nitrogen  as  the  matrix  gas.  The  0.5*  UDMH  standard  was  prepared  by 
injecting  12.6  microliters  of  purified  UDMH  liquid  inJ’o  the  1-liter  bulb  which 
had  been  previously  purged  with  dry  nitrogen  gas.  To  this  mixture,  1.0  ml  of 
ethane  was  added  as  an  internal  standard.  The  0.5*  FDH  standard  was  prepared 
in  the  Identical  manner,  using  14.6  microliters  of  purified  formaldehyde 
dimethylhydrazlne  to  give  the  final  desired  composition.  The  calibration  of 
the  vapor-phase  concentration  of  nltrosodimethylamlne  (NDMA)  was  accomplished 
by  Injection  of  a  liquid  NDMA  standard  in  methanol  such  that  1.0  microliter  of 
the  liquid  standard  corresponds  to  the  same  amount  of  NDMA  which  would  be 
present  in  0.5  milliliters  of  a  0.008*  NDMA  gas  standard.  Additional  NDMA 
standards  in  methanol  were  prepared  to  give  equivalent  responses  of  0.004*, 
0.0016*,  and  0.0008*  NDMA  in  the  gas  phase. 


18 

35 


The  1*  water  vapor  standard  was  prepared  by  adding  6  microliters  of 
liquid  water  to  a  1-liter  gas  sampling  bulb  previously  purged  with  dry 
nitrogen. 

A  0,1<  ammonia  standard  was  prepared  by  injecting  1.0  milliliter  of  high 
purity  ammonia  gas  into  a  1-liter  gas  bulb  purged  with  nitrogen. 

Analysis  of  Standards  by  Gas  Chromatography 

The  method  for  analyzing  the  UDMH,  FOH,  and  NOMA  standards  employed  a  2 
meter  long  glass  column  packed  with  Tenax  GC  60/80  mesh  connected  to  a  flame 
ionization  detector.  The  GC  conditions  were:  column  initial  temperature 
100°C  for  1.0  minute,  temperature  ramp  to  250°C  at  15°C/min,  injection 
port  temperature  120°C,  detector  base  300°C,  carrier  gas  helium  at  30 
ml/min.  Due  to  the  high  concentrations  of  UDMH  and  FDH  used,  the  detector 
range  was  initially  set  at  10*^  ar  i  then  changed  to  10"^  (most 
sensitive)  during  the  elution  of  the  NDMA  peak.  The  Injection  volume  used  in 
this  analyses  was  0.50  milliliters  for  the  gas  standards  and  1.0  microliter 
for  the  NDMA  standard  in  methanol. 

The  method  used  for  analyzing  water  vapor,  oxygen,  and  nitrogen  for  the 
UDMH-atmospheric  tests  used  a  four  foot,  1/8-inch  stainless  steel  column 
packed  with  Tenax  GC  60/80  mesh  connected  in  series  with  a  9  foot  1/8-inch 
stainless  steel  column  packed  with  molecular  sieve  5A  30/60  mesh. 

The  GC  conditions  were:  column  initial  temperature,  40°C  for  1.0 
minute,  temperature  ramp  to  140°C  at  15°C/min,  injection  port  temperature 
100OC,  carrier  gas  helium  at  30  ml/min,  detector  oven  200°C,  filament 
temperature  230°C.  The  injection  volume  used  was  1.0  milliliter. 

The  method  for  analyzing  ammonia  and  water  standards  for  the  hydrazine 
atmospheric  tests  utilized  a  four  foot,  1/8-inch  stainless  steel  column  packed 
with  Tenax  GC  60/80  mesh  and  a  thermal  conductivity  detector.  The  GC 
conditions  were:  column  40°C  for  1.0  minute,  ramp  to  140°C  at  15°C/min, 
injection  port  100°C,  carrier  gas  helium  at  30  ml/min,  detector  oven  200°C 
filament  temperature  250°C.  The  injection  volume  used  was  1.0  milliliter. 
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Calibration 


Calibration  curves  for  NOMA,  water,  and  ammonia  were  prepared  from  the 
gas-chromatographic  results  of  standards  prepared  from  these  commodities. 

These  calibration  curves  are  presented  in  Figures  2  through  5,  respectively. 
Response  factors  for  oxygen  and  nitrogen  were  calculated  using  room  air  as  the 
calibration  gas  (21<  oxygen,  79X  nitrogen  nominal  composition).  Calibration 
for  UOMH  and  FOH  were  performed  using  propane  as  an  internal  standard  and  were 
calculated  on  a  Varian  Vista  6000  gas  chromatography  data  system  in  volume 
percent  of  UDMH  and  FD!!,  respectively. 

The  GC  retention  time  and  detector  responses  for  the  analyzed  gas 
commodities  using  the  Instrument  parameters  described  above  are  Included  in 
Table  II. 

The  GC  retention  time  and  detector  responses  for  the  analyzed  gas 
commodities  using  the  instrument  parameters  described  above  are  included  in 
Table  II. 
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Table  II  Gas  Chromatography  Results 
Atmospheric  Gas  Standards 


Retention  Detector  Response 

Commodity _ Concentration _ Time  (min.) _ (Counts) _ 


Propane 

1%  (int.  std.) 

0.316 

3.01  X  105 

UDMH 

0.5% 

2.995 

9.925  X  105 

FDH 

0.5 % 

4.531 

1.838  X  106 

NDMA 

0.005% 

6.609 

1.057  X  105* 

Water  Vapor* 

1.0% 

0.769 

5.82  X  104 

Oxygen 

21%  (nom.) 

10.946 

3.43  X  106 

Nitrogen 

79%  (nom.) 

12.534 

1.166  X  107 

Ammonia2 

0.1% 

0.487 

4.0  X  103 

Water  Vapor2 

1.0% 

0.737 

1.85  X  105 

1  UDMH  Tests 

2  Hydrazine 

Tests 

*  Attenuation  Change  on  Detector  X  10* 
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Figure  2.  Ammonia  Calibration  Curve  Gas  Chromatography 


VOLUME  %  AMMONIA 


Figure  3.  NOMA  Calibration  Curve  Gas  Chromatography 
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VOLUME  %  NOMA 


Figure  4.  Water  Calibration  Curve  (UDMH)  Gas  Chromatography 


c  01 * ISNOJSIU 1)0133130 
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VOLUME  %  WATER 


Figure  S.  Water  Calibration  Curve  (Hydrazine  Tests)  Gas  Chromatography 
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3.1.3  Nitrogen  Dioxide-Atmospheric  ions 

Gas-phase  mixtures  of  N02  in  helium,  21%  oxygen  in  helium,  and  21% 
oxygen  +  0.05%  carbon  dioxide  in  helium  were  prepared  by  mixing  ultra  pure  gas 
mixtures  provided  by  Linde  Specialty  Gases.  Sample  No.  13  (21%  oxygen  +  1% 
water  vapor  in  helium)  was  prepared  by  humidifying  the  21%  oxygen-helium 
mixture  by  bubbling  the  gas  through  a  gas  washing  bottle  filled  with  distilled 
water  (flow  rate  0.5  1pm)  prior  to  mixing  with  the  nitrogen  dioxide  gas.  Gas 
reaction  mixtures  were  then  fed  into  large  (15-liter)  gas  sample  bags 
(Calibration  Instruments,  Inc.)  which  were  equilibrated  at  an  elevated 
temperature  (35°C)  by  a  mechanical  convection  oven.  The  large  sample  bags 
were  required  for  the  N0x-N0  sample  analysis  to  be  described  below. 

Sample  Analysis 

Gas  aliquots  from  Test  Nos.  11  through  14  were  sampled  at  the  time 
Intervals  detailed  in  Table  I  and  analyzed  for  N0X  an(j  no  concentrations 
using  a  chemiluminescent  analyzer  (Thermoelectron  Corp.).  The  Instrument  was 
calibrated  using  a  440  ppm  NO  standard  prior  to  each  series  of  analysis. 
Chemiluminescent  readings  were  accomplished  by  attaching  a  teflon  line 
connecting  each  of  the  four  sample  bags  to  a  five-port  switching  valve,  the 
outlet  (port  5)  of  which  was  connected  to  the  N0X  analyzer.  Readings  could 
be  quickly  taken  from  all  four  sample  bags  by  sequentially  switching  the 
five-port  valve.  The  sampling  flow  rate  from  the  chemiluminescent  N0X 
analyzer  was  0.5  1pm,  which  accounted  for  the  need  of  a  large  reaction 
container. 

3.2  Results 

3.2.1  Hydrazine-Atmospheric  Reactions 

Hydrazine  decay  results  In  the  1-liter  gas  sampling  bulb  as  measured  by 
the  P0A8  colorimetric  method  are  presented  in  Table  III.  Hydrazine  did  not 
decay  in  the  pure  helium  mixture,  indicating  that  very  little 
autodecomposltion  occurred  on  the  glass  surface.  All  samples  exposed  to  21% 
oxygen  decayed  competely  within  15  minutes.  The  addition  of  0.05%  carbon 
dioxide  seemed  to  slow  or  inhibit  the  reaction.  The  very  rapid  decay  of  the 
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hydrazine  vapor  In  oxygenated  atmosphere  Is  most  certainly  heterogeneous,  and 
Is  probably  due  to  the  large  surface-to-volume  ratio  of  the  1-1  Iter  gas  bulbs 
used  In  these  tests. 
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Table  III  POAB  Colorimetric  Analysis  Hydrazine  Atmospheric  Reaction 


Percent 


Test 

No. 

Initial  Gas 

Composition 

Reaction  Time 
Minutes 

Absorbance 
455  nm 

Hydrazi 

Vapor 

1 

0.5X  N2h4  in  He 

5 

0.90 

0.52 

30 

0.90 

0.52 

60 

0.89 

0.51 

120 

0.90 

0.52 

240 

0.86 

0.49 

2 

0.5%  N2h4  ♦  21X  O2  in  He 

5 

0.015 

0.003 

15 

0.010 

0.002 

3 

0.5%  N2h4  +  21%  02  +  IX 

8 

0.005 

0.002 

H20  In  He 

15 

0.009 

0.002 

4 

0.5X  N2h4  +  21X  02  +  .05% 

2 

0.08 

0.037 

CO2  in  He 

5 

0.108 

0.049 

15 

0.001 

0.002 

8 

0.5%  N2h4  +  21X  02  +  1% 

5 

0.004 

0.002 

H20  +  0.05X  C02  in  He 

15 

0.000 

0.002 

A  graph  of  the  hydrazine  concentration  for  the  five  sample  bulbs  versus 
time  is  shown  in  Figure  6. 
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Results  of  the  water  vapor/nitrogen/oxygen  gas  analysis  for  the  hydrazine 
atmosphere  reaction  mixtures  are  presented  in  Table  IV  and  indicate  that  both 
nitrogen  and  water  vapor  are  formed  in  these  reactions  and  that  ammonia 
formation  reported  by  Stone^  was  well  below  detectable  limits  in  the  present 
study.  Water  vapor  was  formed  In  the  ratio  2  moles  H20  formed  per  1  mole 
N2H4  reacted  which  is  in  excellent  agreement  with  the  results  by  Stone. 
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Table  IV  Gas  Chromatography  Results  Hydrazine  Atmospheric  Reaction 


Test  Initial  Gas 

No.  Composition _ *  Water  *  Oxygen  X  Nitrogen  X  Ammonia 


1 

0.5*  N2h4  in  He 

/\ 

0 

• 

O 
*— * 

0.07 

0.21 

-'.05 

2 

0.5*  N2h4  +  21*  02 

in  He 

0.89 

21.9 

1.05 

'.05 

3 

0.5*  N2h4  +  21*  02 
H2c  in  He 

+  1* 

1.91 

26.0 

0.71 

in 

O 

• 

\ 

4 

0.5*  N2h4  +  21*  02 
C02  in  He 

♦  0.05* 

0.82 

18.0 

0.66 

'.05 

5 

0.5*  N2h4  +  21*  O2 

+  1* 

1.96 

21.6 

1.45 

'.05 

H20  +  0.05*  C02  in 

He 
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UPMH-Atmospherlc  Reaction 


Results  of  the  unsymmetrlcal  dimethyl  hydrazine  atmospheric  reactions  as 
analyzed  by  gas  chromatographic  techniques  are  included  in  Tables  V  through 
IX.  The  UDMH  decay  was  fitted  logarithmically  to  a  pseudo-first  order 
reaction,  and  the  correlation,  slope,  and  rate  constant  of  the  plot  of  the 
logarithm  of  the  UDMH  concentration  versus  reaction  time  are  also  included  in 
these  figures.  The  mathematical  development  and  determination  of  the 
pseudo-first  order  rate  constant  for  UDMH  decay  are  presented  in  Section  3.3. 

The  reaction  of  UDMH  with  dry  CO^-free  air  Is  very  slow  under  the 
conditions  studied,  exhibiting  a  reaction  half-life  of  approximately  257 
hours.  When  the  UDMH-oxgen  mixture  Is  humidified  with  \X  H2O,  the  decay  of 
UDMH  Is  accelerated,  with  a  reaction  half-1  if®  of  24  hours  and  a  concomitant 
increase  in  both  FDH  and  NDMA  reaction  products.  The  addition  of  0.05X  carbon 
dioxide  to  both  the  dry  and  humidified  oxygen  mixture  seems  to  Inhibit  the 
oxidation  of  UDMH  to  FDH.  The  half-life  for  these  reactions  increased  to 
approximately  990  hours  for  the  dry  reaction  and  46  hours  for  the  humid 
reaction.  This  result  was  unexpected  and  not  easily  explained  on  mechanistic 
grounds;  however,  the  hydrazine  decay  also  appeared  to  be  retarded  by  carbon 
dioxide  addition  in  less  than  stoichiometric  amounts. 


32 


49 


Table  V 


Gas  Chromatographic  Analysis  UDMH  +  He  Reaction 
(Test  No.  6) 


Reaction 

Time 

(hr) 

%  UDMH 

*  FOH 

%  NOMA 

%  02 

%  n2 

%  H20 

0 

0.50 

... 

0.02 

0.096 

4.9 

0.49 

0.0013 

< 0.0001 

-  -  - 

--- 

12.6 

0.50 

0.0008 

■'  0.0001 

... 

...» 

.  .  — 

21.8 

0.49 

0.0008 

<0.0001 

... 

... 

... 

68.5 

0.46 

0.0011 

<0.0001 

0.05 

0.149 

0.04 

first  order  kinetic  fit: 


Correlation 

slope 

k 

half-life 


.949 
-.0005 
0.001  hr-1 
693  hr 
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Table  VI 


Gas  Chromatographic  Analysis  UDMH  +  21%  Oxygen  in  Helium  Reaction 

(Test  No.  7) 


Reaction 

Time 

(hr) 

*  UDMH 

%  FDH 

%  NDMA 

X  02 

X  n2 

X  H2O 

0 

0.54 

0.002 

0.0001 

22.59 

0.152 

0.022 

5.2 

0.48 

— 

--- 

11.9 

0.44 

isis  mWf 

0.0002 

— 

-- - 

21.3 

0.45 

0.0003 

22.22 

0.170 

0.085 

45.0 

0.41 

UlWlMlMI 

0.0004 

21.10 

0.509 

0.116 

69.2 

0.43 

0.102 

0.0004 

22.05 

1.212 

0.158 

first  order  kinetic  fit:  Correlation  0.73 

slope  -.0012 

k  0.0026  hr-1 

half-life  257  hr 
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Table  VII 


Gas  Chromatographic  Analysis  UDMH  +  21%  Oxygen  +  1%  Water  Vapor 

In  Helium  Reaction 
(Test  No.  8) 


Reaction 

Tim^ 


( hr) 

%  UOMH 

%  FOH 

%  NOMA 

%  02 

%  n2 

%  H20 

0 

0.49 

0.003 

HMBI 

21.14 

0.134 

1.08 

6.1 

0.42 

0.097 

■■ 

— 

— 

— 

11.8 

0.31 

0.066 

... 

... 

21.1 

0.28 

0.122 

21.31 

1.38 

44.9 

0.13 

0.128 

\Hj$g 

20.71 

0.428 

1.50 

69.1 

0.11 

0.219 

0.0012 

22.05 

1.193 

1.35 

first  order  kinetic  fit: 


Correlation 

slope 

k 

half-life 


0.992 
-.0126 
0.029  hr-1 
24  hr 
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Table  VIII 


Gas  Chromatographic  Analysis 
UDMH  +  2\%  Oxygen  +  0.05<  Caruon  Dioxide 
In  Helium  Reaction 
(Test  No.  9) 


Reaction 

Time 

(hr) 

%  UDMH 

*  FDH 

%  NDMA 

*  o2 

X  n2 

%  h2o 

0 

0.50 

■•mm 

21.95 

0.135 

2.8 

0.50 

0.002 

<0.0001 

— 

... 

— 

10.5 

0.51 

0.004 

0.0002 

.MM 

..  . 

— 

19.7 

0.49 

0.005 

0.0001 

... 

... 

... 

67.6 

0.48 

0.006 

0.0001 

23.35 

0.158 

0.11 

first  order  kinetic  fit:  Correlation  0.83 

slope  -.0003 

k  0.0007  hr-1 

half-life  990  hr 
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Table  IX 


Gas  Chromatographic  Analysis  IJOMH  +  2\%  Oxygen  +  1%  Water  Vapor 
+  0.05?  Carbon  Oioxlde  In  Helium  Reaction 
(Test  No.  10) 


Reaction 

Time 

(hr) 

*  UDMH 

%  FOH 

%  NDMA 

%  02 

%  n2 

%  H20 

0 

0.52 

0.002 

<0.0001 

21.7 

0.21 

1.19 

5.5 

0.43 

0.022 

0.0003 

— 

- — 

— 

12.9 

0.40 

0.038 

0.0004 

... 

— 

— 

21.4 

0.39 

0.103 

0.0004 

21.3 

0.22 

1.19 

45.1 

0.26 

0.125 

0.0006 

20.4 

0.24 

1.14 

69.3 

0.17 

0.172 

0.0005 

21.1 

0.55 

1.02 

In  the  0.5%  UDMH  reaction  with  21%  oxygen,  1%  water  vapor,  and  0.05% 
carbon  dioxide  (Test  No.  10)  0.49  moles  of  FDH  were  formed  per  mole  of  UDMH 
reacted.  This  result  Is  In  conflict  with  the  stoichiometry  expressed  In 
Equation  (3).  In  this  case  0.67  moles  of  FDH  were  formed  per  mole  UDMH 
reacted.  These  results  Indicate  that  the  reaction  between  UDMH  vapor  and 
humid  air  may  react  In  a  different  stoichiometry  ratio  than  the  dry  air 
reaction.  These  results  also  Indicate  that  0.002  to  0.003  moles  of  NDMA  are 
formed  per  mole  of  UDMH  oxidized,  under  the  conditions  studied. 

3.2.3  Nitrogen  Dioxide  Atmospheric  Reactions 

The  N0X  concentrations  of  Test  Nos.  11  through  14  as  a  function  of 
reaction  time  at  35°C  are  recorded  In  Table  X.  The  decrease  in  NO2 
concentration  as  a  function  of  time  Is  pictured  In  Figure  7  and  Is  most  likely 
attributed  to  the  reaction  of  nitrogen  dioxide  with  the  walls  of  gas  sampling 
bag.  It  Is  also  Interesting  to  note  that  the  nitric  oxide  concentration 
decreases  In  the  presence  of  oxygen,  forming  nitrogen  dioxide. 
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REACTION  TIME  IHOURSI 


Table  X  Nitrogen  Dioxide  -  Atmospheric  Reaction  Analysis 


Test 

No. 

Initial  Gas 

Composition 

Reaction  Time 
(Hours) 

N0X 

Percent 

NO 

ppm 

11 

0.25X  NO £  in  He 

1 

0.22 

113 

3 

0.20 

118 

7 

0.20 

138 

24 

0.08 

100 

12 

0.25X  N02  +  21X  02  in  He 

1 

0.22 

107 

3 

0.20 

105 

7 

0.17 

42 

24 

0.04 

8 

13 

0.25%  N02  +  21*  02  +  IX 

1 

0.25 

115 

H20  In  He 

3 

0.22 

95 

7 

0.18 

45 

24 

0.02 

5 

14 

0.25X  N02  +  21X  02  +  0.05X 

1 

0.23 

107 

C02  in  He 

3 

0.20 

90 

7 

0.14 

29 

24 

•  •• 

7 
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3.3  Rate  Constant  Evaluation 


First  order  rate  constants  for  the  UDMH  atmospheric  reactions  were 
calculated  as  follows:  If  the  UDMH  oxidation  rate  with  excess  oxygen  (21X 
0 2  In  air)  Is  first  order  with  respect  to  the  UDMH  concentrations,  then 

let  |uDMHj0  =  initial  propellant  concentration  at  t*0. 

JuDMhJ  *  Propellant  concentration  at  later  time  t. 

the  rate  of  propellant  decay  with  reaction  time  Is  directly  proportional  to 
the  propellant  concentration: 


•d(UDMH]/dt  —  k[UDMH]  K  *  1st  order  reaction  rate  constant 

•d(UDMH]  =  k[UDMH]dt 

UOMH  t 

-  /  d[UDMH]/[UDMH]  =  k  /  dt 

UOMH,  0 

Upon  Integration: 

-ln([UDMH]/[UDMH]o)  =  ln([UDMH]0/(UDMH])  =  kt 
log([UDMH]0/[UDMH])  =  (k/2.303  )t 

log[UDMH]  -  [-kt/2.303]  +  iog(UDMH]0 
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Therefore  a  plot  of  the  log  of  the  UDMH  concentration  versus  time  will  give  a 
slope  of  -  k/2.303. 

The  rate  constant  can  also  be  determined  by  reaction  half-life  which  is 
the  time  required  for  the  reactant  species  (UDMH)  to  decrease  to  half  Its 
Inti al  value. 

for  a  first  order  reaction,  the  half-life  (tjyg)  1$  defined 
at  t  =  t1/a;  c  =  1/2(UOMH]o 

then: 

IOfl(lUDMH]o/0.5lUDMH]o)  =  (k/2.303)!,* 
ti/2  =  o.e«3/k 

Figure  8  shows  a  plot  of  the  log  of  UDMH  concentration  vs  reaction  time  for 
three  conditions  (pure  helium,  21 %  oxygen  ♦  1%  water  vapor,  and  21%  oxygen  ♦ 

1%  water  vapor  ♦  0.05%  carbon  dioxide).  It  Is  Interesting  to  note  that  the 
oxidation  reaction  Is  retarded  In  the  presence  of  carbon  dioxide.  In  this 
latter  case,  a  rate  constant  of  0.0152  hr"1,  with  a  reaction  half-life  of  46 
hours  Is  reported  (as  compared  to  a  half  life  of  24  hours  In  the  absence  of 
carbon  dioxide). 

Using  the  stoichiometry  expressed  In  Equation  (3),  for  every  mole  (or 
liter)  of  UDMH  decreased  during  the  oxidation  reaction,  2/3  mole  (0.67)  of  FDH 
Is  formed.  In  addition,  0.003  mole  of  HDMA  Is  formed  from  every  mole  UDMH 
oxidized,  under  the  conditions  studied  In  our  tests. 
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Figure  8.  Atmospheric  Reaction  Rates  of  UDMH 


0.5%  UOMH  IN  Ha 

0.5%  UMOH ♦  1%  H,0  *21%  0,  IN  Ha 

0.5%  UMOH  ♦  1  %  H,  0  ♦  2 1  %  0,  ♦  0.05%  C02  IN  Ha 


_ _ _ _ _ 

K  -  0.00 1  HR' 


2Q  30  40  50 


REACTION  TIME  (HOURS) 


4.0  Condensation  of  Fireball  Components 


Since  the  combustion  products  In  the  fireball  are  cooled  by  the 
entrainment  of  ambient  temperature  air  upon  fireball  liftoff  from  ground  zero, 
there  Is  a  possibility  that  certain  components  will  condense  out  of  the 
cloud.  The  mathematical  method  used  to  determine  whether  condensation  occurs 
compares  the  vapor  pressure  of  the  component  as  a  pure  liquid  to  Its  partial 
pressure  In  the  fireball.  If  the  partial  pressure  of  the  component  ever  meets 
or  exceeds  Its  vapor  pressure  as  a  pure  liquid,  condensation  will  occur. 

These  condensed  materials  will  coalesce  and  be  deposited  on  the  ground  below 
the  approximate  point  of  condensation. 

The  calculations  require  the  Input  of  the  following  data: 

1)  Atmospheric  pressure 
?.)  Atmospheric  temperature 

3)  Fireball  temperature 

4)  Volume  and  molar  composition  of  fireball 

5)  Heat  capacities  and  Enthalpies  of  fireball  components 

With  this  data  Input,  the  calculations  can  be  performed  for  either  N2o4  or 
N2H4.  These  chemicals  were  singled  out  because  of  their  potential  for 
high  concentrations  and  their  hazardous  nature.  Hydrazine  was  tested  for 
condensation  In  an  A-50  va/>or  mixture  because  it  has  a  lower  vapor  pressure  at 
ambient  temperature  than  UDMH,  and  is  therefore  more  likely  to  condense. 

The  first  calculation  necessary  Is  to  find  the  temperature  of  the  fireball  as 
a  function  of  the  amount  of  air  entrained,  (eq.  12) 


T(°K)  ^  (H«  +  H„,N.*)/(CpfNf  +  Cp  airNair)  (1?) 
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where: 


T  (°K)  ■  Temperature  of  fireball  in  degrees  Kelv.n 
Hpg  »  Total  heat  content  of  original  fireball  in  calories 
(13)  Hair  s  6.76T  *  1/2  (0.606xl(T3)T2  -2044  +  1/3  (0.13xl(T6)T3 


(14) 


(15) 

(16) 

(17) 


where  T  is  in  degrees  Kelvin 


N.„  =  (PV/RT)X1000 


Cpf  =  IN„%  C 

n»1 


pn 


P  *  atmospheric  pressure  in  atm 
V  =  volume  of  air  entrained  in  m3 
R  =  0.0821  m3atm/kg>|0K 
T  ■  atmospheric  temperature,  °K 
C=  Average  heat  capacity  of  fireball 

rr 

7 -jffi  mole  percent  of  component  n 
c pns  heat  capacity  of  component  n 

Volume  of  fireball  at  liftoff  (m3) 
Tj=  Temperature  of  fireball  at 
liftoff  °K 

Cp  =  6.76  +  (0.606x10'3)T  +  (0.13x10‘6)T2  ca1/77°K, 

sir*  \ 

where  T  is  temperature  of  air  in  degrees  kelvin 


N,  =  (PV,/RT,  )X1000 


After  the  temperature  of  the  fireball  as  a  function  of  the  amount  of  air 
entrained  is  determined,  both  the  vapor  pressure  of  the  pure  liquid  component, 
(a  function  of  temperature),  and  the  partial  pressure  of  the  component,  (a 
function  of  concentration),  can  be  calculated. 


Vapor  Pressure  Equation  for 

Pv  1  (mm  Hg)  =  lofl^l-e-SOS  -(653.88/T (°K))  +  0.0479T(°K)  -4.989X10'*T2(°K2))  (18) 


where  Py 


1 


=  vapor  pressure  of  pure  N2H4  in  mmHg 


T(°K)  =  Fireball  temperature  as  per  eq.  12 


Vapor  Pressure  Equation  for  N2O4: 


Pva  =  log  [9.824  -  2332/T(°K)  -84567/T2(°Ky)] 


where  Py^  =  Vapor  pressure  of  pure  N2O4  In  mnHg 
T(°K)  =  Fireball  temperature  as  per  eq.  12 


Partial  pressure  for  either  component: 

P,  =  [Ni/(Ni  +  N«i,)JXP  (20) 

where:  p,  =  partial  pressure  of  component  of  Interest 

Ni  =  total  moles  of  component  of  interest 
N.k  ~  total  moles  of  air  entrained 
P  =  atmospheric  pressure  (mm  Hg) 

The  equations  for  vapor  pressure  are  only  valid  for  temperatures  at  or  below 
the  normal  boiling  point  of  the  chemical  species  in  question.  This  poses  no 
problem  because  condensation  cannot  occur  above  the  boiling  point.  The 
equation  for  partial  pressure  Is  valid  for  all  situations. 


With  these  equations.  It  can  be  determined  whether  or  not  condensation 
will  occur.  This  Is  done  by  plotting  both  the  vapor  pressure  and  the  partial 
presure  of  the  component  of  Interest  as  a  function  of  the  number  of  moles  of 
air  entrained.  If  the  two  lines  ever  Intersect,  condensation  will  occur  at 
that  point.  In  the  following  section  two  example  calculations,  one  for  an 
excess  oxidizer  case  and  the  other  for  an  excess  fuel  case,  are  carried  out. 
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Excess  Oxidizer  Case:  O/F  °  5.1 


Data:  Fireball  radius  at  ground  zero 

«  35.35*m 

composition*:  CO  =  0.020  mole  % 

n2 

=  0.1387  mole  % 

C02  *  0.027  " 

OH 

*  0.022  " 

H2  =  0.025  " 

°2 

rH 

O 

• 

O 

II 

H20  =  0.144  " 

no2 

*  0.565  “ 

NO  »  0.005  " 

H 

=  0.012  " 

Cloud  radius  at  stabilization2  650  m 

Total  heat  content  (Hpg)  »  1.40  x  10*^  calories* 

Initial  temperature  of  fireball  *  810°K* 

Assume  ambient  conditions  to  be  1  atm  &  20°C  (298°K) 

1.  S.  Prince.  Atmospheric  Oispersion  of  Hypergolic  Liquid  Rocket  Fuel 
Phase  I:  Source  Characterization,  Final  Report  30  September  1982. 

2.  H.E.  Cramer  Co.  Report  on  Cloud  Growth  and  Dispersion  Model  for 
Hypergolic  Fuels.  Personal  Communication,  8  April  1983. 


Calculations: 


Volume  of  fireball  at  ground  zero  ■  4/3t r  (35.35m)3 

V  »  185.05  x  103m3 


From  eq.  16, 
From  eq.  13, 


Hr  =  (latm  •  188.05Xl03m3/(0.0821  •  810°K)]X1000  =  2.78X10* 

H  <  «  6.76(298)  +  1/2  (0.606xl0‘3) (298)2  + 

1/3  (0.13x10"®)  (298)  -2044  »  -1.466  cal/g7£ 


From  eq.  15,  *  C 


pp  *  9.8  cal /g?20K 


From  eq.  17,  C  =  6.76  +  (0.606X10"3) (298)  +  (0.13X10'6)(298)2 

"aSr 


=  6.95  cal/j?£0K 


With  these  values,  the  equation  for  vapor  pressure  as  a  function  of  the  number 
of  moles  air  entrained  is  derived  using  equations  12  and  19. 

(21)  T(°K)  =  [1.4X1010  +  (-1.486  N.(,)!/(9.8X2.78X10®  +  6.95N„1 

Pv,  =  log1  (9.824  -  2332/TCK)  -  84587/T2(°K2)] 


The  equation  for  partial  pressure  as  a  function  of  the  number  of  moles  air 
entrained,  using  equation  20,  Is: 


(22)  P*o  =  I(0.585X2.78XlO*)/(  (0.565X2.78X10*)  +  N.„)]X760  (mm  Hg) 


See  Figure  9  for  the  plots  of  equations  21  and  22.  From  the  graph.  It  Is 
obvious  that  the  two  lines  will  never  Intersect.  Note  the  scale  change  for 
the  pressure  values.  There  Is  a  difference  of  approximately  three  orders  of 
magnitude  betwee.i  the  two  pressures.  For  this  particular  situation  there  will 
be  no  N2O4  condensing  out  of  the  cloud. 


A  VAPOR  PRESSURE 
★  PARTIAL  PRESSURE 


.21 


-y-r 

i*'  k 


1.#  1.70  2.0*'*  1( 

PRESSURE  (psi) 


Excess  Fuel  Case:  0/F  =  .204 


Oata:  Fireball  radius  at  ground  zero  =  34.4m* 


CO 

--  0.037 

mole  X 

NO 

=  0.007 

mole  % 

C02 

=  0.031 

II 

n2 

=  0.193 

ll 

H2 

=  0.035 

H 

OH 

=  0.032 

II 

H 

=  0.017 

II 

°2 

=  0.057 

II 

h2o 

=  0.203 

II 

UOMH 

=  0.134 

ll 

n2h4 

=  0.254 

II 

HF8  *  2.304  x  1010cal(D 

Temp  of  Fireball  =  1046°k(*) 
Assume  ambient  temp  =  20°C  & 

Cloud  radius  at  stabilization 


1  atm  pressure. 
=  736  m 


Calcul ations: 

Volume  of  Fireball  at  ground  zero  *  4/3 Tr' (34.4m) 3 

=  1.705  x  lO^m^ 


From  eq  16 
From  eq  13 


Nf  =  [(latmx  1.705X1 05m3)/(0.082i xi 046° K))X1 000  =  1.99X106  moles 
H.„  =  -1.466  cal/gN 


From  eq  15  CpF  =  7  97  ca|/gN°K 
From  eq  17  Cp  alr  =  6.95  ca|/gN°K 


Now,  using  equations  12,  18,  20 


and 


T(°K)  =  (2.304X10’°  -  1.466XN8if)/'(7.97X1. 99X10®  +  6.95XNa„) 


68 


GRAM-MOLES  OF  AIR  ENTRAINED 


Pnch.  =  log*1  (-6.506  -  653.88/T(°K)  +  0.0479T(°K)  -  4.989X10  5T2(eK2)) 


(23) 


P»n,  =  ( (0.2S4X  1 ,49X  1 0®)/(0.254X  1 .49X10®  +  N8„))X760  (mm  Hg)  (24) 

See  figure  10  for  a  graph  of  equations  23  and  24.  Equation  23  Is  only  valid 
for  temperatures  below  700°R,  which  Is  the  boiling  point  of  N2H4.  From 
this  graph  it  Is  evident  that  condensation  will  not  occur  in  this  case 
either.  There  Is  about  three  orders  of  magnitude  difference  between  the  two 
pressures. 


5.0  CONCLUSIONS 


Atmospheric  reaction  rate  data  for  reactive  fireball  spec'es  to  be 
Incorporated  Into  a  computer  atmospheric  dispersion  model  are  presented  in 
Table  XI.  The  table  lists  the  conditions  in  which  the  rate  constants  and 
reaction  half-lives  apply. 


Table  XI  Atmospheric  Reaction  Rate  Constants 
for  Fireball  Dispersion  Modeling 


Concentration 

Major 

Fireball 

in 

Meteorological 

Rate 

Half- 

Reaction 

% 

Species 

Atmosphere 

Condition 

Constant 

Life 

Products 

Yield 

Hydrazine 

0-1  ppm 

dry  or  humid  air 

0.2  ppm-1 
min-1 

2  hr* 

H2°2 

NA 

Hydrazine 

Above  1  ppm 

dry  air 

0.064  hr-1 
0.141  hr-1 

10.8  hr 

n2,h2o 

100 

Hydrazine 

Above  1  ppm 

humid  air 

4.9  hr 

n2,h2o 

100 

UOMH 

0-1  ppm 

dry  or  humid  air 

0.15  ppm-1 
min-1 

0.2  hr* 

NDMA 

60 

UDMH 

Above  1  ppm 

dry  air 

0.0007  hr-1 

990  hr 

FDH** 

67 

UDMH 

Above  1  ppm 

humid  air 

0.015  hr-1 

46  hr 

FDH** 

67 

N0? 

10  ppm 

humid  air 

5.5X1Q4!2 

mole-2sec"l 

5  hr 

HN03(g) 

100 

no2 

50  ppm 

humid  air 

5.5xlQ4!2 

mole"2sec_1 

0.2  hr 

HN03(g) 

100 

no2 

Above  50  ppm 

humid  air 

5.5X10412 

mole-2sec"l 

<0.2  hr 

HN03(1) 

100 

no2 

All  Cone. 

dry  air 

NO  REACTION 

*  Clean, 

Unpolluted  Atmosphere 

Reacted) 

**  NDMA  Formed  in  0.2%  Yield  (.002  Moles  NDMA  Found  Per 

1  Mole  UDMH 

NO  =  Not  Available 
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USER'S  MANUAL  FOR  THE  HYPERGOLIC  ACCIDENT  RELEASE 
MODEL  (HARM)  COMPUTER  PROGRAM 
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SECTION  1 
INTRODUCTION 


1 . 1  BACKGROUND 

The  U.  S.  Air  Force ,  and  particularly  the  Air  Force  Engineering 

and  Services  Center  at  Tyndall  Air  Force  Base,  has  pursued  the  analysis  of 

existing  dispersion  models  and  the  development  of  new  approaches  for 

modeling  dispersion  in  support  of  the  Titan  II  safety  engineering  program. 

The  Hypergollc  Accidental  Release  Model  (HARM)  is  the  result  of  an  effort 
to  develop  a  computer  program  incorporating  source  characterization, 
cloud-rise  and  dispersion  model  algorithms  appropriate  for  describing  the 
transport  and  dispersion  downwind  from  a  hypergollc  reaction  involving  the 
fuel  and  oxidizer  propellants  of  the  Titan  II  missile.  The  hypergollc 
reaction  results  in  the  formation  of  a  thermally  buoyant  cloud  which  rises 
vertically  to  a  height  related  to  the  amount  of  heat  generated  by  the 
reaction,  which  is  a  function  of  the  amount  of  propellant  components 
involved  in  the  accidental  release.  The  existing  model  (Ocean  Breeze/Dry 
Gulch)  used  to  assess  propellant  hazards  (U.  S.  Air  Force,  1983)  is  not 
suited  for  modeling  dispersion  downwind  from  elevated  volume  sources.  The 
HARM  computer  program  Is  based  on  source  characterization  algorithms  for 
hypergollc  reactions  developed  by  Martin  Marietta  Aerospace  (Prince,  1982, 
1983),  buoyant  cloud  models  developed  by  Briggs  (1970)  and  dispersion 
models  for  rocket  launches  (Bjorklund,  ex_  al^  ,  1982)  developed  by  the  H.  E. 
Cramer  Company. 


1.2  PURPOSE 

The  technical  objective  of  the  work  described  in  this  report  is 
to  provide  a  computer  program  that  predicts  ground-level  concentrations  of 
toxic  materials  resulting  from  a  catastrophic  reaction  of  hypergollc 
propellants.  The  purpose  of  this  report  is  to  provide  documentation  of  the 
following: 
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•  A  description  of  the  algorithms  used  to  define  the  source 
characteristics  of  selected  catastrophic  hypergollc  reactions 
Involving  Titan  II  propellants 

•  A  mathematical  description  of  the  cloud-rise,  dispersion 
models  and  other  mathematical  formulas  used  in  the  HARM  code 

e  User's  instructions  for  the  HARM  computer  program 


1.3  ORGANIZATION  OF  THE  MANUAL 

The  main  body  of  the  manual  contains  five  sections.  Section  2 
contains  a  brief  overview  of  the  major  components  of  the  HARM  code,  the 
accident  scenarios  the  code  is  designed  to  accommodate,  and  the  meteoro¬ 
logical  and  source  model  input  parameters  required  to  operate  the  program. 
The  algorithms  used  in  the  program  to  characterize  the  source,  to  describe 
the  reaction  products  and  decay /production  rates,  and  to  predict  the 
cloud-rise  of  the  thermally  buoyant  cloud  are  described  in  Section  3. 
Section  4  contains  a  mathematical  description  of  the  atmospheric  dispersion 
models  Implemented  by  the  computer  code.  Section  5  provides  user  instruc¬ 
tions  for  the  HARM  program. 

Appendix  A  of  this  manual  contains  example  executions  of  the  HARM 
program.  A  program  listing  and  magnetic  tape  containing  the  program  have 
been  provided  to  Martin  Marietta  Aerospace  under  separate  cover. 
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SECTION  2 

OVERVIEW  OF  THE  HARM  COMPUTER  PROGRAM 


The  catastrophic  hypergolic  reaction  of  Titan  II  rocket  propellants 
results  in  the  nearly~instantaneous  formation  of  a  relatively  large  cloud 
of  hot,  buoyant  reaction  products  near  ground  level  which  subsequently 
rises  and  entrains  ambient  air  until  the  temperature  and  density  of  the 
cloud  reach  an  approximate  equilibrium  with  ambient  conditions.  Depending 
on  the  amount  of  fuel  (A-50)  and  oxidizer  (NTO)  involved  in  an  accident  and 
the  atmospheric  stability,  the  stabilization  height,  defined  as  the  point 
the  cloud  ceases  buoyant  rise,  can  vary  from  10' s  of  meters  to  heights 
exceeding  i  km.  The  cloud  may  contain  toxic  components  of  unreacted 
propellants  (hydrazine,  unsymmetrical  dlmethythydrazine,  or  UDMH,  and 
nitrogen  tetroxide,  or  NTO)  in  addition  to  the  products  of  the  reaction. 

Toxic  products  of  the  reaction  Include  dimethylnltrosamine  (NDMA) ,  methyl 
amine,  formaldehyde,  hydrogen  cyanide,  ammonia  and  formaldehyde  dime thy 1- 
hydrazone  (FDH) .  The  HARM  computer  program  is  designed  to  calculate  peak 
concentrations  and  dosages,  peak  deposition  due  to  precipitation  scavenging 
and  isopleths  of  dosage,  concentration  and  deposition  due  to  precipitation 
scavenging  downwind  from  catastrophic  reactions.  The  program  calculations 
are  automated  for  considering  the  five  major  toxic  components  of  hydrazine, 
UDMH,  NTO,  NDMA  and  FDH.  The  HARM  program  is  designed  to  assist  the 
Missile  Potential  Hazard  Team  (MPHT)  in: 

•  Planning  hazard-safety  programs  for  catastrophic  hypergolic 
accidents  involving  Titan  II  propellants 

•  Real-time  support  of  hazard-safety  operations  involving 
Titan  II  propellants 

•  Post-accident  environmental  analyses  of  hypergolic  accidents 
involving  Titan  II  propellants 
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While  the  program  is  designed  for  considering  hypergolic  accidents  involving 
Titan  II  propellants,  it  can  easily  be  modified  to  consider  other  types  of 
hypergolic  reactions. 


2.1  COMPONENTS  OF  THE  HARM  COMPUTER  PROGRAM 

Figure  2-1  is  a  schematic  diagram  showing  the  major  components  of 
the  HARM  computer  program.  Requisite  meteorological  Inputs  to  the  computer 
program  are  obtained  from  the  vertical  profiles  of  wind  direction,  wind 
speed,  air  temperature,  atmospheric  pressure  and  dew  point  or  relative 
humidity  between  the  earth's  surface  and  3000  m.  This  information  is 
obtained  primarily  from  rawinsonde  measurements  made  at  selected  upper-air 
stations  by  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  and 
by  meteorological  stations  at  selected  USAF  air  bases.  Rawinsonde  measure¬ 
ments  are  routinely  made  twice  daily  at  0000  Greenwich  Mean  Time  (GMT)  and 
1200  GMT  at  these  stations.  The  HARM  program  can  be  configured  to  accept 
the  rawinsonde  data  from  tape  and  disc  files  or  from  TTY  Input  devices  such 
as  a  modem.  As  shown  in  Figure  2-1,  the  rawinsonde  data  can  be  manually 
edited  to  provide  for  any  changes  in  the  vertical  profiles  that  weather 
forecasters  assigned  to  the  MHPT  consider  to  have  occurred  between  the  time 
the  rawinsonde  measurements  were  made  and  the  time  of  the  accident. 
Similarly,  meteorological  observations  from  the  accident  site  can  manually 
be  entered  to  reflect  the  observed  near  ground-level  conditions  at  the 
scene. 


The  HARM  program  is  controlled  by  operator  input  and  internal 
management  routines  based  on  the  operator’s  response  to  plain-language 
queries  displayed  on  a  CRT  terminal.  In  Figure  2-1,  this  complex 
interactive  function  ii  simply  designated  by  the  box  labeled  CRT  Program 
Control.  Provision  nas  been  made  in  the  computer  program  for  selection  of 
calculations  for  the  following  5  major  accident  scenarios: 

•  Above-ground  slow  leak 
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•  SCENARIO  TYPE 

•  PRINCE  (1982)  FIREBALL 
MATHEMATICS 


DOSAGE/ CONCENTRATION  MODELS 
AND  ALGORITHMS 

•  NASA/REEDM  MULTI-LAYER 
DOSACE/CONCENT RATION  MODELS 
AND  ALGORITHMS 

•  PEAK  CENTERLINE  DOSACE/ 
CONCENTRATION 

•  DOSACE /CONCENTRATION 
ISOPLETKS 


CLOUD  RISE,  MATERIAL 
DISTRIBUTION  AND  CHEMISTRY 

•  BRIGCS  (1970)  INSTANTANEOUS 
CLOUD  RISE  MODEL 

•  PRINCE  (1983)  CLOUD 
CHEMISTRY 

•  NASA/REEDM  MATERIAL 
DISTRIBUTION  ALGORITHMS 


PRECIPITATION  SCAVENGING 
MODEL  AND  ALGORITHMS 

•  NASA/REEDM  PRECIPITATION 
SCAVENGING  MODEL/ ALGORITHMS 

•  PEAK  CENTERLINE  DEPOSITION 

•  DEPOSITION  ISOPLETHS 


CRAPHICS  ROUTINE 


PRINT  DISPLAYS 


PLOT  DISPLAYS 


FIGURE  2-1.  Schematic  diagram  illustrating  major  components  of  the  HARM 
Computer  Program. 


•  Above-ground  major  spill  (nearly-instantaneous  detonation) 

•  In-silo  accident  with  the  silo  door  closed  and  either  a  slow 
leak  ultimately  resulting  in  a  hypergolic  reaction  or  an 
instantaneous  hypergolic  reaction 

•  In-9llo  accident  with  the  silo  door  open  and  either  a  slow 
leak  ultimately  resulting  in  a  hypergolic  reaction  or  an 
instantaneous  hypergolic  reaction 

•  Stage  ejection  with  a  subsequent  t  losion  at  a  user  specified 
height 

The  operator  is  also  asked  to  supply  the  amount  of  fuel  (A-50)  and  oxidizer 
(NTO)  involved  in  the  accident  and  the  specific  coordinates  of  the  accident 
location.  After  the  operator  has  responded  to  the  these  queries,  the 
program  automatically  selects  a  proper  set  of  source  inputs  to  be  used  in 
the  calculation  of  the  dimensions  of  the  fireball  created  by  the  hypergolic 
reaction  and  the  heat  of  the  reaction  for  use  in  the  cloud-rioe  calculation. 
The  cloud-rise  model  then  calculates  the  spatial  position  and  dimensions  of 
the  cloud  as  a  function  of  time  and  distance  from  the  accident  until  the 
stabilization  height  is  reached.  The  vertical  distribution  of  material  in 
the  various  layers  of  the  atmosphere,  where  the  layer  boundaries  are  given 
by  the  heights  at  which  rawinsonde  measurements  are  available,  is  also 
automatically  calculated.  The  concentration  of  material  in  the  cloud  is 
used  in  an  algorithm  based  on  the  work  of  Prince  (1983)  to  determine  the 
half-lives  and  decay/production  constants  to  be  used  in  the  subsequent 
dispersion  calculations.  At  this  point,  the  rawinsonde  meteorological  data 
and  the  results  of  the  cloud-rise,  cloud  dimension  and  material  distribution 
algorithms  are  output  to  a  nrinter.  The  operator  also  has  the  option  to 
output  a  plot  display  of  the  vertical  profiles  of  wind  direction,  wind 
speed,  temperature  and  virtual  potential  temperature  as  well  as  a  simple 
projection  of  the  stabilized  cloud.  The  operator  and/or  weather  forecaster 
then  have  the  option  o',  modifying  the  default  values  selected  and  calculated 
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by  the  program  to  represent  the  major  meteorological  layer  structure 
parameters  (the  height  of  the  base  and  top  of  an  elevated  inversion  layer, 
for  example)  and  the  turbulence  parameters  that  will  be  used  in  the  disper¬ 
sion  calculations.  After  the  final  selection  of  model  input  parameters  has 
been  made  by  the  operator,  the  program  performs  either  the  selected 
dosage/concentration  or  precipitation  scavenging  (deposition)  calculations 
using  the  multi-layer  dispersion  models  described  in  Section  4.  When  these 
calculations  are  complete,  the  results  are  printed  and,  at  the  option  of 
the  operator,  plotted.  If  the  dosage/concentration  option  was  selected, 
the  print  output  includes  peak,  concentration  at  1-km  intervals  downwind 
from  site  of  the  accident,  and  the  total  dosage  and  time-mean  concentration 
for  the  period  of  interest  at  these  points.  The  operator  can  request  plots 
of  these  results  for  each  pollutant  versus  distance  and/or  plots  of  the 
isopleths  of  concentration,  dosage,  and  time-mean  concentration  on  maps  or 
map  overlays.  If  the  option  to  calculate  deposition  due  to  precipitation 
scavenging  is  selected,  peak  centerline  maximum  possible  or  time-dependent 
peak  centerline  deposition  is  printed  and,  at  the  user's  option,  plotted. 

The  user  can  also  request  isopleth  plots  of  deposition.  Example  of  the 
print  and  plot  output  produced  by  the  HARM  program  are  included  with  the 
example  problems  described  in  Appendix  A. 

Although  not  shown  in  Figure  2-1,  there  are  three  major  modes 
that  the  program  user  can  choose  for  making  calculations  using  the  HARM 
code  (operational,  research  and  production).  The  operational  mode  is 
designed  for  assisting  the  Missile  Potential  Hazard  Team  in  evaluating 
hazards  and  automatically  selects  many  of  the  required  input  parameters, 
although  the  program  user  is  provided  an  opportunity  to  modify  the  selected 
default  values.  When  the  research  mode  is  selected,  more  information  is 
usually  input  by  the  operator.  For  example,  the  operator  can  specify 
values  of  the  turbulence  parameters  at  each  height  where  rawinsonde  data  is 
available  rather  than  using  the  default  turbulence  profiles.  The  production 
mode  of  the  HARM  code  can  he  used  to  process  multiple  rawinsonde  soundings 
which  are  read  from  a  tape  or  disc  file.  While-  the  production  mode  can  be 
run  interactively  I rom  a  CRT  terminal,  the  primary  putpose  of  the  production 


mode  Is  to  facilitate  batch  processing  of  multiple  cases  without  operator 
attention. 


It  should  be  noted  that  the  HARM  computer  program  is  intended  for 
use  in  predicting  the  environmental  impact  resulting  from  the  nearly 
Instantaneous  hypergolic  reaction.  In  some  accidents  there  may  be  residue 
material  in  the  silo  or  on  the  open  ground  not  consur  d  in  the  reaction. 
Material  from  these  ground-level  sources  may  be  emitted  continuously  over  a 
period  of  hours.  The  HARM  computer  program  does  not  consider  the  environ¬ 
mental  impact  of  these  sources. 


2.2  ACCIDENT  SCENARIOS  AND  SOURCE  PARAMETERS 

The  five  major  scenarios  of  hypergolic  accidents  that  the  hARM 
computer  program  can  process  are  listed  In  Section  2.1  above.  The  major 
source  parameter  affecting  the  outcome  of  the  calculations  is  the  ratio  of 
fuel  to  oxidizer  involved  in  the  accident.  The  algorithm  describing  the 
resulting  products  from  the  user  input  -f  the  total  weight  of  fuel  and 
oxidizer  involved  in  the  accident  is  described  In  Section  3.1.  In  those 
accidents  that  involve  stoichiometric  quantities  of  fuel  and  oxidizer  (no 
excess  fuel  or  oxidizer) ,  vhere  is  of  course  no  excise  fuel  or  oxidizer 
available  for  downwind  transport  in  the  chemical  cloud  produced  by  tha 
accident.  In  those  accidents  where  the  hypergolic  reaction  is  nearly- 
instantaneous  (regardless  of  whether  the  accident  is  an  above-ground  major 
spill,  in-silo  accident  with  the  silo-door  closed  or  open,  or  stage  ejection), 
there  is  never  any  excess  hydrazi-.e  available  f t  r  downwind  transport. 

Finally,  all  the  accidents  are  assumed  to  occur  at  ambient  atmospheric 
pressure  except  those  that  occur  in  the  silo  with  the  door  closed,  where 
the  pressure  is  assumed  to  b._  12. o  atmospheres.  Based  on  these  assumptions, 
the  HARM  model  computer  program  treats  accidents  involving  an  above-ground 
slow  leak  and  an  in-silo  slow  leak  with  the  silo  door  open  in  the  9ame 
manner,  since  both  ate  assumed  to  occur  at  a  pleasure  of  1  atmosphere. 
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Similarly,  the  program  treats  above  ground  major  spills  and  the  instan¬ 
taneous  hypergolic  reaction  occurring  in  the  silo  with  the  silo  door  open 
in  the  same  manner. 

The  stage  ejection  accident  scenario  used  in  the  computer  program 
assumes  that  an  instantaneous  explosion  occurs  in  the  silo  and  an  explosion 
of  the  second  6tage  of  the  Titan  II  occurs  at  a  user  defined  height.  The 
user  is  also  required  to  input  the  total  amount  of  fuel  and  oxidizer 
involved  in  the  accident.  The  program  assumes  that  approximately  80 
percent  of  the  total  propellant  weight  is  involved  in  the  silo  explosion 
and  that  approximately  20  percent  is  involved  in  the  explosion  at  the  user 
defined  height.  In  this  scenario,  the  user  is  also  asked  if  the  6ilo  door 
is  open  or  closed.  If  the  door  is  closed,  the  in-silo  reaction  is  assumed 
to  occur  at  a  pressure  of  12.6  atmospheres  and  the  reaction  aloft  at  a 
pressure  of  1  atmosphere.  Both  reactions  occur  at  a  pressure  of  1  atmosphere 
if  the  silo  door  is  open.  The  cloud-rise  is  independently  calculated  for 
each  source.  In  the  case  of  the  ln-sllo  explosion,  cloud-rise  begins  from 
a  point  at  ground-level  while  the  cloud-rise  from  the  in-air  detonation 
begins  at  the  user  designated  detonation  height.  The  amount  of  material  in 
each  atmospheric  layer  from  both  sources  is  calculated  using  an  algorithm 
defined  in  Section  3  and  then  summed  during  the  dispersion  calculations. 


2.3  METEOROLOGICAL  INPUT  PARAMETER  REQUIREMENTS 

As  noted  above,  the  primary  meteorological  input  to  the  HARM 
computer  program  u  in  the  form  of  rawinsonde  observations.  Each  level  of 
lc formation  (standard,  mandatory  a..<.  significant  levels)  in  the  rawinsonde 
data  stream  (,vth  observation  level  in  th<-  model  descriptions  in  Section  3 
and  is  used  to  obtain  wind  and  temperature  profiles.  The  HARM  code  is 
eurrentlv  di-sign-d  to  perform  dispersion  calculations  in  two  major,  meteoro¬ 
logies  1  ly-c!.  tlned,  layers.  The  base  of  the  lower  layer  Is  assumed  to 

hr  at  tne  earth's  surface  and  t'ae  top  of  the  layer  Js  assumed  to  be  located 


at  the  base  of  an  elevated  inversion  (top  of  the  surface  nixing  layer). 

The  boundaries  of  the  upper  layer  (L*»2)  are  set  by  the  operator.  If  the 
concentration/dosage  of  the  material  assumed  to  be  trapped  in  the  elevated 
inversion  is  ot  interest,  for  example,  the  operator  would  set  the  lower 
boundary  of  the  second  layer  equal  to  the  top  of  the  first  layer  and  the 
top  of  the  second  layer  to  the  height  of  the  top  of  the  elevated  Inversion. 
The  boundaries  of  the  two  major  layers  must  coincide  with  a  height  ot  level 
reported  in  the  rawlnsonde  message.  If  the  user  designates  a  boundary  of  a 
major  layer  other  than  a  Kth  height,  the  program  automatically  selects  the 
Kth  height  closest  the  user  input  value. 

The  selection  of  the  boundaries  of  the  two  major  layers  is 
critical  to  the  outcome  of  the  concentration  and  dosage  dispersion  calcula¬ 
tions.  The  gaseous  products  in  the  cloud  produced  by  the  hypergolic 
accident  are  assumed  to  be  reflected  towards  the  earth's  surface  at  the 
tops  of  the  major  boundaries  and  upward  at  the  bases  of  the  major  boundaries 
(including  the  earth's  surface)  in  the  concentration  and  dosage  calculations. 
The  boundaries  of  these  two  major  layers  are  also  used  in  the  specification 
of  the  vertical  profiles  of  the  turbulence  parameters  used  in  the  concentra¬ 
tion,  dosage  and  deposition  calculations  as  indicated  in  Section  3.  In  the 
absence  of  user  input  of  the  layer  boundaries,  the  program  automatically 
sets  the  top  of  the  lower  layer  equal  to  the  Kth  rawlnsonde  observation 
level  closest  to  twice  the  cloud  stabilization  height.  Also,  the  height  of 
the  base  of  the  upper  layer  is  defaulted  to  equal  the  top  of  the  lower 
layer  and  the  program  sets  the  top  of  the  upper  layer  to  the  Kth  observation 
level  closest  to  3000  m.  Rawlnsonde  measurements  at  heights  greater  than 
3000  m  are  not  used  by  the  HARM  code.  This  default  depth  of  the  surface 
mixing  layer  is  used  because  the  use  of  a  Kth  observation  level  between  the 
cloud  stabilization  height  and  twice  the  stabilization  height  to  represent 
the  depth  of  the  mixing  layer  tends  to  maximize  the  calculated  ground-level 
vapor  concentrations. 

In  general,  program  users  should  select  tha  height  of  the  base  of 
an  elevated  inversion  to  represent  the  top  oe  t*>e  lower  (L“ll  layer. 
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Because  the  top  of  the  lower  layer  is  defined  in  the  program  as  an  upper 
bound  to  turbulent  mixing,  the  user  should  select  a  height  above  which  the 
temperature  either  increases  or  remains  constant  for  at  least  50  to  100  m. 
Bec-use  the  primary  purpose  of  the  HARM  code  is  to  assist  in  estimating 
ground-level  environmental  effects  produced  by  hypergolic  accidents,  the 
user  should  generally  ignore  surface-based  inversions  or  elevated  inversions 
with  tops  less  a  few  hundred  meters  above  the  surface  in  selecting  the  top 
of  the  lowest  layer. 


2.4  BASIS  OF  THE  HARM  COMPUTER  PROGRAM 

The  computer  program  described  in  this  manual  incorporates  many 
of  the  dispersion  model  algorithms  and  concepts  employed  in  the  REEDM 
(Rocket  Exhaust  Effluent  Diffusion  Model)  Computer  Program  developed  for 
the  National  Aeronautics  and  Space  Administration  (NASA)  for  assessing  the 
environmental  impact  of  exhaust  products  from  the  normal  and  abnormal 
launches  of  launch  vehicles  (Bjorklund,  ££  £l .  ,  1982).  The  REEDM  model  is 
the  product  of  over  a  decade  of  development  in  preparation  for  real-time 
support  of  Space  Shuttle  launches  and  is  currently  used  to  support  launches 
of  the  Space  Shuttle  at  Kennedy  Space  Center.  The  task  of  predicting  the 
dispersion  downwind  from  a  hypergolic  accident  is  in  many  respects  similar 
to  the  task  of  predicting  dispersion  of  the  ground  cloud  produced  by  rocket 
launches  because  both  scenarios  result  in  the  buoyant  rise  and  formation  of 
an  elevated  source  with  relatively  large  horizontal  and  vertical  dimensions. 
As  noted  in  Section  1.1,  models  developed  for  predicting  the  downwind 
dispersion  from  ground-level  continuous  point  sources  are  not  generally 
applicable  because  of  the  large  initial  dimensions  of  the  cloud,  because 
the  source  is  formed  nearly  instantaneously,  and  because  the  dispersion 
takes  place  in  a  relatively  deep  layer  of  the  atmosphere.  These  factors 
were  considered  in  the  development  of  the  REEDM  computer  ptugram  and  have 
been  incorporated  in  the  HARM  models.  We  also  note  that  NASA,  and  par¬ 
ti  'ul.irly  Langley  Research  Center,  has  conducted  sampling  measurement 
P'  ograni.*;  using  ground-based  and  aircraft  .sampling  equipment  in  conjunction 
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with  normal  launches  of  Titan  and  Space  Shuttle  launches  over  the  past 
decade.  These  measurements  (Gregory,  e_t  al.,  1974a,  1974b;  1976:  Woods, 
e£  al_. ,  1979;  Wornom,  et  al. ,  1979;  Sebacher,  et  al.  ,  1980  and  others)  have 
been  used  to  verify  and  improve  the  performance  of  the  REEDM  computer 
program  predictions  of  dosage  and  concentration.  In  general  the  predicted 
concentrations  have  been  within  a  factor  of  two  of  the  measured  concentra¬ 
tions  and  the  models  tend  to  overpredict  ground-level  concentrations, 
yielding  conservative  estimates  of  the  ground-level  environmental  impact. 
Thus,  the  decision  was  made  to  adapt  many  of  the  concepts  Incorporated  In 
the  NASA  REEDM  computer  Program  for  use  in  the  HARM  Computer  Program  for 
the  following  reasons: 

•  REEDM  has  proven  successful  in  predicting  concentrations 
downwind  from  the  release  of  large  thermally  buoyant  clouds 

•  REEDM  contains  the  requisite  algorithms  for  calculating 
cloud-rise,  material  distributions  and  other  source  charac¬ 
teristics  of  large  thermally  buoyant  clouds 

•  REEDM  contains  the  requisite  dispersion  algorithms  for 
treating  the  downwind  transport  and  dispersion  of  clouds 
with  initially  large  horizontal  and  vertical  dimensions 

•  REEDM  uses  available  meteorological  Input  information 
(ravinsonde  data)  describing  meteorological  conditions  in 
and  above  the  surface  mixing  layer 

e  REEDM  has  moderate  computer  memory  and  run-time  requirements 
suited  for  real-time  support  of  hazard-safety  operations. 

Those  concepts  and  the  algorithms  employed  in  the  HAIM  Computer  Program  to 
describe  the  source  characteristics  of  hypergolic  reactions  are  described 
in  detail  in  Section  3  and  4  below. 
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SECTION  3 

HARM  SOURCE,  CLOUD  RISE  AND  METEOROLOGICAL  ALGORITHMS 

3.1  SOURCE  CHARACTERISTIC  ALGORITHMS 

As  mentioned  above,  the  algorithms  describing  the  source  charac- 

i 

teristics  included  in  the  HARM  model  are  bused  on  the  work  of  Prince  (1982, 
1983).  The  following  discussion  is  principally  based  on  these  reports. 

The  mixing  of  the  hypergolic  rocket  propellant  A-50  and  NTO  during  an 
accidental  spill  or  missile  tank  rupture  converts  chemical  energy  into 
thermal  energy.  The  thermal  energy  is  used  to  heat  the  hypergolic  combustion 
products,  vaporize  excess  unreacted  propellant  and  to  heat  the  environment 
in  the  vicinity  of  the  accident.  The  thermal  energy  of  the  fireball  formed 
by  the  hypergolic  reaction  and  the  size  of  the  fireball  are  required  in  the 
HARM  model  for  calculating  the  stabilization  height  of  the  cloud  of  reaction 
products  and  the  chemical  composition  of  the  cloud  at  the  stabilization 
height. 


The  thermochemical  analysis  performed  by  Martin  Marietta  Aerospace 
(Prince,  1982)  proceeded  as  follows: 

•  Fireball  combustion  products  were  identified  and  adiabatic 
flame  temperatures  calculated  using  theoretical  thermo¬ 
dynamic  combustion  properties  of  the  propellant  and  the 
gaseous  reaction  products  for  stoichiometric  conditions 

•  In  cases  where  the  oxidizer  to  fuel  ratio  was  non- 
stoichiometric ,  the  chemical  reaction  was  created  as  a 
non-equilibrium  condition  where  the  thermal  energy  of  the 
fireball  was  reduced  by  that  energy  required  to  heat  and 
vaporize  the  excess  propellant 
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•  Tne  thermal  eaergy  of  the  fireball  was  further  reduced  by 
accounting  for  radiative  heat  transfer  to  the  surrounding 
environment 

A  description  of  the  algorithms  used  in  the  HARM  program  for  calculating 
the  chemical  composition  of  the  material  in  the  fireball ,  the  thermal 
energy  available  for  rise  of  the  fireball  and  the  dimensions  of  the  fireball 
are  given  below. 


3.1.1  Chemical  Composition  and  Fireball  Temperatures 

Determination  of  the  fireball  combustion  products  and  the  calcula¬ 
tion  of  fireball  adiabatic  temperatures  resulting  from  stoichiometric 
reactions  of  A-50  and  NTO  was  accomplished  by  Martin  Marietta  Aerospace 
using  a  computer  program  for  the  calculation  of  complex  chemical  equilibrium 
composition  (NASA  SP-273).  The  program  reiteratively  solves  equations  that 
minimire  the  Gibbs  free  energy  of  the  chemical  reaction  products  and 
maintains  a  mass  balance  between  the  chemical  reactants  (bypergoiic 
propellants)  and  chemical  products  (combustion  products,  oxidation  products 
and  unreacted  propellants).  The  output  for  stoichiometric  combinations  of 
A-50  and  NTO  reacting  at  1  and  12.6  atmospheres  is  reproduced  in  Table  3-1. 

Flame  temperatures  of  the  fireball  provided  to  the  H.  E.  Cramer 
Company  were  calculated  under  an  adiabatic  assumption  where  conductive, 
convective  and  radiative  heat  losses  to  the  environment  were  not  considered. 
Thus  the  heat  of  reaction  caused  by  combining  A-50  with  NTO  at  25#C  was 
used  to  heat  the  reaction  products  in  the  fireball  from  25°C  to  the  temper¬ 
ature  at  time  t,  or  T{t}.  The  reaction  equation  is: 


(a l  +  a2^  C. 696  H5.39  N- 


*  a3  N2°4 


P  +AH 

J  j  reaction 

J 


16 

8:> 


(3-1) 


TABLE  3-1 


THERMODYNAMIC  PROPERTIES  AND  MOLE  FRACTIONS  FOR  STOICHIOMETRIC  MIXING 
OF  A-50  AND  NTO  (AFTER  PRINCE,  1982) 


where 


a^  *  moles  of  reacted  liquid  hydrazine 


a^  *  moles  of  reacted  liquid  UDMH 
*»  moles  of  reacted  liquid  0^ 


3^  *  moles  of  the  product 


AH 


reaction 


-  heat  evolving  from  the  chemical  reaction 


1.54  x  10^  cal  mole  *  of  reacting  A-50 


For  an  adiabatic  process 


/T{t) 

c’  d  T'  (products) 
p 


(3-2) 


where 


low  pressure  heat  capacity  of  the  products  (cal  mole  *  #K  S 


A  «•  BT{ t }  +  C (T { t ))2  +  D(T(t})3 


A,  B,  Ci  D  ■  coefficients  defined  below 


T{ t )  ■  temperature  in  ®K  at  time  t 


Integrating  Equation  (3-2)  yields  the  expression 


AH j (reactant  s) 


-  AHj (products)  +  ATft)  ♦  (B/2)(T(t})‘ 
♦  (C/3) (Tf t } ) 3  +  (D/4) (T { t } ) 4  -  E 


(3-3) 
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Scffiffiff'WAVj .-.m  a*  *• «.  *-w  *  „  * 


where 

E  -  A(298)  +  (B/2X298)2  +  (C/3)(298)3  +  (D/4)(298)A  (3-4) 

Values  of  the  coefficients  A,  B,  C,  D  and  E  required  to  solve  Equation 
(3-3)  are  given  in  Table  3-2  and  are  based  on  the  values  recommended  by 
Prince  (1982)  in  his  Table  IV.  The  4th  degree  polynomial  in  T{ t )  is  solved 
for  T{t}  using  Newton's  method.  The  value  of  T{t)  of  the  fireball  at  the 
time  of  lift-off  from  the  ground  T{t-tg}  is  used  in  calculating  the  fireball 
size  and  heat  content  of  the  cloud  at  liftoff  from  the  ground  as  outlined 
below. 

3.1.2  Fireball  Size  and  Heat  Content 

Prince  (1982)  based  the  calculations  of  the  fireball  size  and 
heat  content  of  the  cloud  at  lift-off  from  the  ground  on  a  mathematical 
model,  developed  by  Sandia  Laboratories,  that  expresses  the  reaction  time 
tD  in  seconds  for  the  mass  of  hypergolic  propellant  W.  to  mix  and  react  to 

D  D 

completion  as 

tg  -  0.6845  Wb1/6  (3-5) 

where  the  units  of  are  in  kilograms.  The  radius  of  the  fireball  varies 
as  a  function  of  the  fireball  temperature  T{t)  according  to  the  expression 

/  3  W  R'  T{ t }  \1/3 

r  { t )  **  I  - I  (3-6) 

V  p  *w  / 


17 


92 


F.V' 

1 


m 

I 

co 


c 

o 


<0 
p 
c r 
U4 


c 

<v 

tH 

u 

•»H 

Uh 

o 

o 


(U 

cO  *—t 

'©  I 


W  €0 

a 


on  i 
o  V 


co 

LO  I 
O  4) 

rH  rH 

o 

•  fl 


eg 
*2 
CM  | 
O  4» 

rH  rH 

o 

•  a 


r-s 

* 

* 


a.  a>  ^ 

e  oo  ^ 

Ci  Co 
H  <0  w 


CO 

>3 

O 

© 

IM 

>3 

© 

U1 

•H 

© 

rH 

© 

CO 

On 

ON 

CM 

© 

N© 

CN 

PM 

© 

PH 

CO 

CM 

PH 

m 

O 

1 

1 

i 

o. 

CO 

ON 

CM 

© 

© 

00 

© 

fM 

ON 

rH 

CM 

•“* 

1 

1 

i 

♦ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

rH 

CM 

H 

~ 

CM 

CM 

CM 

CN| 

CO 

CO 

CM 

CM 

•o 

CO 

<o 

fv 

O' 

ao 

1 

1 

1 

o 

o 

o 

© 

© 

00 

o 

CM 

o 

© 

© 

On 

© 

1 

1 

1 

o 

© 

00 

m 

• 

• 

• 

» 

• 

U 

o 

rH 

rH 

1 

O' 

U0 

00 

MT 

o 

P^» 

NO 

On 

© 

NO 

CO 

o 

PN. 

ao 

Ph 

CM 

00 

ph» 

00 

CM 

© 

1 

1 

i 

CM 

*3 

O 

<r 

Ph 

CM 

cO 

rH 

n© 

oo 

CO 

H 

1 

1 

o 

CO 

o 

o 

© 

o 

o 

o 

O 

vO 

«H 

© 

CM 

fH. 

• 

* 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

r 

© 

o 

© 

o 

o 

© 

a 

o 

© 

CM 

o 

o 

CM 

t 

i 

1 

i 

1 

i 

i 

i 

i 

i 

i 

i 

n© 

o 

O' 

•o 

00 

oo 

sT> 

m 

nO 

© 

© 

CO 

© 

•A 

O) 

© 

© 

NO 

CM 

1 

1 

1 

UO 

rH 

© 

© 

CO 

-3 

© 

m 

in 

NO 

rH 

N© 

1 

1 

i 

i-h 

© 

rH 

CO 

CM 

H 

rH 

rn 

© 

.o 

CO 

© 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

rH 

O 

© 

© 

o 

© 

o 

rH 

NO 

rH 

o 

i 

o 

CO 

© 

^3 

© 

CM 

O' 

m 

CM 

rH 

V© 

in 

1 

1 

00 

© 

\© 

CM 

r-» 

© 

CM 

00 

CO 

CM 

vO 

© 

CO 

^3 

1 

1 

i 

^3 

CO 

O' 

^3 

O' 

<r 

tn 

c- 

rH 

© 

m 

ST' 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

« 

• 

© 

© 

© 

© 

© 

© 

m 

© 

o 

*n 

NO 

rn 

rH 

O 

© 

o 

© 

o 

o 

© 

© 

© 

o 

© 

© 

o 

o 

O 

© 

g 

© 

© 

o 

o 

g 

© 

o 

g 

o 

o 

o 

r- 

ph 

r- 

pH. 

fH 

o 

m 

m 

m 

i 

( 

i 

<o 

CO 

i/o 

cO 

CO 

CO 

CO 

4A 

CO 

m 

CM 

H 

»H 

i 

I 

i 

i 

1 

1 

1 

1 

i 

1 

i 

1 

i 

i 

1 

1 

1 

CO 

CO 

© 

CO 

CO 

CO 

r> 

© 

CO 

© 

© 

co 

CO 

CO 

r* 

pH 

o 

CM 

pH 

pH 

ph 

o 

fH 

© 

fH 

PH 

fH* 

CM 

CM 

o 

CM 

CM 

CM 

CM 

© 

CM 

© 

CM 

CM 

CM 

rH 

rH 

rH 

ON 

© 

n© 

CM 

-3 

CO 

© 

in 

•o 

«n 

© 

o 

•3 

wO 

co 

uo 

O' 

ON 

© 

CM 

CO 

o 

NO 

^H 

© 

CO 

\o 

o 

o 

© 

ph 

© 

© 

NO 

O 

«3 

r» 

ON 

o 

CM 

CM 

<r 

•© 

*3 

CM 

fs. 

rH 

ON 

CM 

O 

PH 

CM 

H 

t 

CM 

1 

o\ 

1 

UO 

m 

1 

CM 

CM 

CM 

1 

/~s 

O'. 

rH 

rH 

— M 

? 

<r 

*3 

X 

a; 

3C 

o 

CM 

o 

X 

X 

CM 

CO 

CM 

p 

CM 

o 

O 

CM 

CM 

© 

CM 

X 

CM 

CM 

© 

o 

X 

z 

5 

z 

u 

o 

32 

X 

X 

is 

z 

© 

o 

Z 

5 

z 

z 

a0 

f*'- 

rH 

H 

<0 

r 

• 

© 

•H 

CM 

CO 

U-N 

NO 

ON 

*  « 

O*) 

CO 

UO 

%© 

r*- 

00 

O' 

rH 

^4 

•H 

«r4 

M 

TO 

<0 

<D 

CO 

CD 

<0 

CO 

<0 

CD 

<T> 

a 

«0 

<0 

<0 

CD 

<9 

Species  are  gaseous  unless  otherwise  noted. 


where 


P  ■  ambient  atmospheric  pressure  (Pascals) 
R'  -  8.31432  J°K~l  mole'1 


•*  average  molecular  weight  of  the  cloud  (kg  mole  weight) 


?  %  ^ 

‘  ^ 

Mwj  =  molecular  weight  of  species  j 

and  Tit)  is  the  fireball  temperature  at  time  t.  The  final  radius  of  the 
cloud  at  cloud  lift-off  r  is  obtained  when  t  is  equal  to  t  .  That  is, 

K  R 


(3-7) 


rR  '  rU'tB) 


(3-8) 


when  T{t}  ■  T{t*t  }.  The  value  of  T{t"t„}  is  calculated  from  the  time 

O  D 

history  of  the  fireball  temperature  and  is  obtained  numerically  by  using  a 
fourth-order  Runge-Kutta  method  to  solve  the  following  nonlinear  differ¬ 
ential  equation: 


dT 

dt’ 


AH  (reactants)  -  AH  (products) 
r  r 


-4tt  c  o 


B 


1.667  W, 


1/6 


1  t'  R' 


4  ti  P  M 


w — 


2/3 


T(t’)14/3 


(3-9) 


19 


94 


where 


t'  =  nondimensional  time 
'  t/CB 

e  •  emissivity  (set  equal  to  1  for  black  body  radiation) 
o_  *  Boltzmann's  constant 

D 

-  5.67  x  10-8  W.  J  s~ 1  m"3  #K“4 

o 


(3-1C) 


The  total  effective  heat  available  for  cloud  rise.  H,  is  calculated  from 
the  expression 


H  -  c'  W.  T{  t“t„ } 
p  b  B 


(3-11) 


3.1.3  Propellant  Decay  and  Decay  By-Products 


Non-stolchiometric  ratios  of  oxidizer  to  fuel  allow  excess 
propellants  to  remain  in  the  cloud  of  reactants.  Both  components  of  the 
fuel  as  well  as  the  oxidizer  also  react  with  chemical  species  found  in  the 
atmosphere  to  produce  species  not  originally  present  in  Che  fireball.  The 
HAPJ1  model  assumes  that  the  decay  of  material  In  the  cloud  after  the  time 
of  cloud  stabilization  is  of  the  exponential  form 


exp  { -k±c ]  (3-12) 
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where  is  the  rate  coefficient  of  decay  for  the  ith  species.  Rate 
coefficients,  half-lives,  major  products  and  yields  of  the  fuel  and  oxidizer 
are  given  in  Table  3-3  based  on  the  values  reported  by  Prince  (1983).  At 
concentrations  below  1  ppm,  both  hydrazine  and  UDMH  have  second  order  rate 
coefficients  which  are  much  larger  than  at  higher  concentrations  and  are 
dependent  on  concentration.  They  both  disappear  rapidly  as  the  concentration 
decreases  below  1  ppm.  Since  the  reaction  products  of  hydrazine  are  the 
relatively  harmless  species  and  1^,  changing  to  when  the  concentra¬ 

tion  is  below  1  ppm,  the  HARM  model  uses  the  rate  coefficient  0.141  hr  *  to 
describe  the  decay  of  hydrazine  at  all  concentration  levels.  The  principal 
exposure  limits  for  hydrazine  are  10  to  30  ppm  and  the  overestimation  of 
hydrazine  at  concentrations  less  than  1  ppm  will  not  affect  hazard  calcula¬ 
tions.  The  exposure  limits  for  UDMH  (30  to  100  ppm)  are  also  much  larger 
than  the  concentration  level  where  the  reaction  rate  becomes  second  order 
and  the  use  of  the  linear  rate  of  0.015  hr  ^  shown  in  Table  3-3  will  not 
affect  the  hazard  calculations  for  UDMH.  The  major  products  of  the  UDMH 
reaction  with  air  and  water  vapor  are  MDMA  and  FDH  at  concentrations  above 
1  ppm.  The  production  of  these  species  is  modeled  using  the  expression 

PK(i}  -  EkU}(1  -  exp(-kit)  ]  (3 

where  E  { i }  is  the  efficiency  (yield)  of  the  conversion  from  UDMH  to  NDMA 

K  -1 
or  FDH  and  the  value  of  equals  0.015  hr 

The  decay  of  the  oxidizer  is  of  third  order  as  shown  by  the  units 
2  -2-1 

(l  mole  s  )  given  in  Table  3-3.  The  half-life  varies  from  less  than 
12  minutes  for  concentrations  greater  than  50  ppm  to  greater  than  5  hours 
for  concentrations  less  than  10  ppm.  Principle  exposure  limits  for  NO2  and 
^0^  are  in  the  range  from  1  to  5  ppm.  For  use  in  the  HARM  model  calcula¬ 
tions,  a  first  order  decay  coefficient  appropriate  at  a  concentration  of  l 
ppr,  was  estimated  from  considering  the  rate  of  conversion  of  NO^  to  HNC^  in 
the  chemical  reaction 


TABLE  3-3 


ATMOSPHERIC  REACTION  RATE  CONSTANTS  (PRINCE,  1983)  USED  IN  HARM 


Species 

Concentration 

(ppm) 

Rate 

Constant 

k 

Half-Life 

(Hours) 

Major 

Reaction 

Products 

Yield  EL 
(%)^ 

Hydrazine 

0-1 

0.2  (a) 

2  (b) 

H2°2 

NA  (e) 

Hydrazine 

>  1 

0.141  (c) 

4.9 

n2,  h20 

100 

UDMH 

0-1 

0.15  (a) 

0.2  (b) 

NDMA 

60 

UDMH 

>  1 

0.015  (c) 

46 

FDH 

67 

NDMA 

0.2 

no2 

10 

5.5  x  10A  (d) 

5 

HNO- 

3 

100 

no2 

50 

5.5  x  104  (d) 

0.2 

HNO- 

100 

(a) 

Units 

are  ppm  *  min  * 

(b) 

For  a 

clean  unpolluted  atmosphere 
are  hr  *. 

(c) 

Units 

(d) 

Units 

are  i  mole  s 

(e) 

Not  available. 
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(3-14) 


4N0,  +  2H2<3  +  02  -  4HN03 


The  differential  equation  expressing  the  rate  of  conversion  is: 


d(N02] 

dt 


k3;H2o][No2]2 


(3-15) 


where 


■  third-order  decay  -lte 

c  c  »2  i  -2  -1 

■  5.5  x  10  £  mole  s 

and  the  brackets  indicate  the  concentration  of  the  enclosed  chemical 
compound.  Separating  variables,  integrating  and  solving  for  the  time  tj 
required  for  the  initial  concentration  of  nitrogen  dioxide  [N02]0  to  decay 
to  1  ppm  yields  the  expression 


*N02*O  "  1  ppm 
Cl  "  k3(H20][N02y  1  ppm] 


(3-16) 


The  equivalent  first  order  coefficient  k^  required  to  produce  the  same  time 
t^  for  the  decay  of  N02  to  a  concentration  of  1  ppm  is 


ln|N02]0 


(3-17) 
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where  tN02^o  i8  «*pr«asad  in  ppm.  Equation  (3-17)  ie  used  in  the  HARM 
model  to  obtain  the  decay  coefficient  for  the  oxidizer. 

3.2  CLOUD  RISE  AND  ASSOCIATED  ALGORITHMS 


3.2.1  HARM  Cloud  Rise  Model 


The  determination  of  the  stabilized  height  of  the  buoyant  cloud 
produced  by  hypergolic  reactions  ia  an  important  factor  in  the  dosage/ 
concentration  calculations  because  the  suximum  dosage /concent rat ion  calcu¬ 
lated  at  the  earth's  surface  is  approximately  inversely  proportional  to  the 
cube  of  the  stabilized  height.  According  to  the  algorithm  described  in 
Section  3.1  above,  the  cloud  of  buoyant  gas  or  "fireball"  from  the  hypergolic 
reaction  forms  and  is  ready  to  lift  from  the  ground  over  a  time  period  of 
about  3  seconds.  Experience  in  predicting  the  cloud  rise  from  the  launch' 
of  solid  rocket  launches  (Bjorklund,  at  al. ,  1982)  with  similar  formation 
times  indicates  that  the  Instantaneous  cloud  rise  model  based  on  the  work 
of  Briggs  (1970)  yields  stabilization  heights  that  are  in  good  agreement 
with  observed  heights. 

The  time  for  the  cloud  produced  by  the  hypergolic  reaction  to 
reach  a  height  z^  in  a  stable  atmosphere  is  given  by  the  expression 


where  t.  is  constrained  to  be  less  than  the  cloud  stabilization  time  t*  and 


-0.5 


(3-10) 


t*  *  ns 


a  •  stability  parameter 


T  Lt 

«2 

g  »  gravitational  acceleration  (9.8  os  ) 


(3-20) 


—  »  vertical  gradient  of  virtual  potential  teoperacure  (®K  rr  *) 

z  between  the  surface  and  height 

y  *  entrainment  coefficient  (0.64) 


F 


1 


H 


initial  buoyancy  tern 


it 


4  ttc  Tc 
P 


(3-21) 


total  effective  heat  release  (cfll)  obtained  ln>m  the  nwrre 
algorithm 


Cp  ■  specific  heat  of  air  at  constant  pressure  (0,?4  cal  g  1  °K  ') 
T  •  ambient  air  temperature  (  K)  near  the  surface 
p  *  ambient  air  density  (g  m  near  the  surface 


r 


R 


initial  radius  of  the  fireball  (m)  obtained  from  the  source 
algorithm 


The  cloud  stabilizut Ion  height  is  defined  as 


C  I  -  .  ’ I 


2  3 

100 


According  to  the  above  formulas,  the  final  stabilization  height  is 
related  to  the  vertical  gradient  of  virtual  potential  temperature  measured 
over  the  same  height.  The  value  of  z^  in  Equation  (3-22)  must  therefore  be 
determined  through  iteration.  The  iteration  process  is  initiated  by 
assuming  that  z^  lies  in  the  first  height  interval  (K*l)  above  the  surface 
bounded  by  the  rawlnsonde  observations  at  (k"l)  and  (k*2)  and  solving 
Equation  (3-22)  for  z^  with 


A$ 

Az 


(3-23) 


If  z^  exceeds  z^,  the  iteration  continues  using  the  virtual  potential 
temperature  from  the  next  height  Kth  observation  level  with  the  vertical 
gradient  A<t/Az  estimated  from  the  least-squares  approximation 


(3-24) 


Providing  that  <  3000  m,  Che  ptogram  finds  a  value  of  z^  within  an 
interval  (z^  ^  S  z^  <  z^).  At  this  point  the  program  assumes  that  the 
gradient  of  virtual  potential  temperature  in  this  height  interval  is  linear 
and  linearly  interpolates  to  determine,  within  ±10  m,  the  value  of  z^. 
Momentum  flux  is  not  considered  in  this  application  of  the  cloud-rise  model 
because  inclusion  of  the  momentum  term  does  not  affect  the  final  cloud 
stabilization  height. 


The  cloud-rise  model  described  above  was  theoretically  derived 
for  use  in  thermally  stable  atmospheric  layers  (A$/Az  >0).  A  model  based 
on  similar  considerations  is  easily  derived  for  an  adiabatic  atmosphere. 
However,  the  adiabatic  model  predicts  the  cloud  to  rise  over  all  time,  and 
the  rate  of  rise  approaches  zero  asymtot ically  at  longer  times.  Experience 
shows  that  the  height  at  which  the  rate  of  rise  determined  from  the  adiabatic 
model  becomes  negligible  for  practical  purposes  can  be  predicted  using  a 


small  value  of  A$/Az  in  the  stable  model  described  above.  For  this  reason, 

-U 

the  program  sets  A#/Az  equal  to  3.322  x  10  when  the  value  developed  from 

-4 

the  rawinsonde  measurements  indicate  A$/Az  i9  less  than  3.322  x  10 


3.2.2  Dimensions  and  Material  Distribution  in  the  Stabilized 
Cloud 


The  dispersion  models  described  in  Section  4  are  derived  under 
the  assumptions  that  a  vertical  finite  line  source  can  be  used  to  represent 
the  source  of  material  in  each  of  the  K  layers  defined  by  the  rawinsonde 
observation  heights  and  that  the  radius  of  the  cloud  at  the  stabilization 
height  Zj  is  consistent  with  that  used  in  deriving  the  cloud-rise  model, 
e.g. : 


r 


I 


Y  Z!  *  "R 


(3-25) 


Under  these  assumptions,  the  alongwind  and  crosswind  dimensions  of  the 
cloud  in  the  Kth  layer  are 


(zTK  +  ZbK)/2 


(3-26) 


(3-27) 

(3-28) 


TK 


height  of  the  top  of  the  Kth  layer 
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height  of  the  base  of  the  Kth  layer 


'BK 


The  HARM  model  assumes  that  the  source  material  in  each  Kth  layer  is 

uniformly  distributed  in  the  vertical  and  Gaussian  distributed  alongwind 

and  crosswind  in  the  horizontal  plane  consistent  with  the  Gaussian  model 

approach  used  in  the  dispersion  models.  The  initial  standard  deviations  of 

the  alongwind  a  and  crosswind  o  source  dimensions  used  in  the  cloud 
xo  yo 

expansion  formulas  described  in  Section  4  are  defined  as  follows: 


°*o{K)  ■  VK)  ■  V2-15 


(3-29) 


under  the  assumption  that  the  concentration  at  one  radius  from  the  cloud 
centerline  is  one-tenth  of  the  maximum  concentration  at  the  cloud  centerline. 
While  the  vertical  distribution  in  each  Kth  layer  is  assumed  uniform  with 
height  in  the  layer,  the  fraction  of  the  total  material  in  the  cloud 
residing  in  the  Kth  layer  at  the  time  of  cloud  stabilization  is  derived 
under  the  assumption  that  the  material  is  Gaussian  distributed  over  the 
entire  vertical  extent  of  the  cloud.  Thus  the  fraction  by  weight  of  the 
material  of  each  species  F { K }  determine!  from  the  algorithm  described  in 
Section  3.1  above  is  given  by  the  expression 


F(K)  -  -  pUBK» 


(3-30) 


where 


P  { z  } 


(3-31) 
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3.2.3  Spatial  Position  of  the  Cloud  at  Stabilization  Time 


The  spatial  position  in  the  plane  of  the  horizon  of  the  portion 
of  the  cloud  produced  by  the  reaction  and  residing  in  the  Kth  layer  at  the 
stabilization  time  t*  is  determined  in  the  HARM  computer  program  from  the 
wind  directions  and  wind  speeds  obtained  from  the  rawinsonde  and  the  time 
from  Equation  (3-18)  required  for  the  cloud  to  rise  through  the  tops  of  the 
various  K  layers.  The  distance  R^  and  bearing  0^  at  time  t*  of  that 
portion  of  the  cloud  in  the  Kth  layer  is  defined  by  the  following 
expressions: 


R 


CK 


[xK  ♦  uK(t*-t  )  sin  eK]2  ♦  [yK  ♦  ^<e*~tp >  cos  fy2 


1/2 


(3-32) 


eCK  “  "  Can 


-1 


yK  *  uK(t*~tp)  °K 


xK  +  uJC(t*-t  )  cos  0R 


(3-33) 


where 


V  ~  *K-1  ~  sin(<!)sK) 


yK-i  -  «k  cos(nK) 


(3-34) 

(3-33) 


u  *  mean  wind  speed  in  the  Kth  layer 
K 


(UTK  +  UBK)/2 


(3-36) 


u  =  wind  speed  measured  at  the  top  of  the  Kth  layer 
1 K 


u  =  wind  speed  measured  at  the  base  of  the  Kth  layer 
1»K 


t  =  time  the  cloud  passes  through  the  top  z  of  the  Kth  layer, 

P  obtained  by  setting  z  in  Fquation  (3-18)  to  z 

K  IK 
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0  ■  mean  wind  direction  in  the  Kth  layer 

l\ 


^0TK  +  6BK^/2 


6  "  winc*  direction  measured  at  the  top  of  the  Kth  layer 

8  ■  wind  direction  measured  at  the  base  of  the  Kth  layer 

EK 

For  the  layers  below  the  Kth  layer  containing  the  stabilization  height  z_. 


«K  *  Cp  UK 


A  m  6 

VSK  K 


For  the  layer  containing  z^,  and  $sK  are  defined  as 


-  (  2i~zbk\ 

kVw-ibJ  “bk 


sK  2 


(- 


+  2  'bk 


Finally,  for  all  K  layers  above  z7(z  >z_), 

X  Dl\  X 


rck  *  rck(z=7-i> 


9CK  ~  ^CK  ^ 2_  z  i  ■* 


(3-37) 


(3-38) 


(3-39) 


(3-40) 


(3-41) 


(3-42) 


(3-43) 


30 

105 


3.3 


TURBULENCE  PROFILE  ALGORITHM 


The  HARM  dispersion  models  use  profiles  of  the  standard  deviations 
of  the  azimuth  wind  angle  and  elevation  angle  o^.  as  prime  predictors  of 
lateral  and  vertical  cloud  growth  downwind  from  the  point  of  cloud  stabiliza¬ 
tion.  The  computer  program  calculates  default  vertical  profiles  of  turbulence 
which  can  be  adjusted  by  the  user.  The  algorithm  used  to  calculate  the 
profiles  begins  by  determining  a  reference  standard  deviation  of  the  wind 
azimuth  angle  <j^{to=600s},  assumed  representative  of  a  measurement  made 
over  a  ten-minute  period  at  the  lowest  height  reported  on  the  rawinsonde 
observation  (usually  from  2  to  10  m  above  the  surface). 


3.3.1  Calculation  of  the  Default  Value  for:  o!„{t  =600s) 
-  AR  o 

The  calculation  method  (Goldford,  et^  a_l,  1977),  based  on  the 
application  of  similarity  relationships  outlined  by  Golder  (1972),  assumes 
that 


JAR‘ 


0'.n{  T  “600s) 


ov{600s) 


f  { B } 


In  ( — ^  - 
\  z  I 
\  °/ 


Ri  *  0 


^{Ri } 


(3-44) 


where 


o 


V 


T 

O 


f  (B> 


k 


A 


standard  deviation  of  the  crosswind  component  of  the  wind 

mean  wind  speed  at  the  measurement  height  of  o^r 

reference  time  for  the  measurement  of  o'  and  o., 

AR  V 

function  of  the  bulk  Richardson  number  B 

Von  Kantian' s  constant  =  0.4 
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z  »  geometric  mean  height  of  the  layer  of  interest 
zq  *  roughness  length  (default  value  set  to  0.2  o) 
iHRi)  “  function  cf  the  Richardson  number  Ri 
In  the  program,  values  of  f { B }  and  ip{Ri)  are  obtained  from  the  expressions 


For  Ri  <  0,  the  following  equation  is  solved  by  Newton's  method  to  obtain 
a  value  of  £  for  use  in  Equation  (3-46) 
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1  -  c 


16S2  Jin(z/2o)  +  0.30864  -  2(in(l+OI  -  *n(l+S2)  ♦  2tan“l£ 
where  C  is  defined  by  Equation  (3-47). 

Finally,  for  RiB0,  0^(1^=6003}  is  calculated  from  the  relationship 


-B 


(3-50) 


o , _ { t  *600s} 
AK  O 


48.816 


Ri  *  0 


(3-51) 


The  program  does  not  permit  o'  (t  =600s)  to  be  greater  than  0.349  radians 

AK.  C 

(20  degrees). 


3.3.3  Vertical  Profiles  of  a!  and  o' 

_ A _ E 

The  dispersion  models  described  in  Section  4  use  mean  values  of 

and  a^,  in  the  two  major  meteorological  layers  (L-1,2).  In  the  case 

where  the  user  enters  values  of  o'  and  o'  at  each  k  rawinsonde  observation 

A  E 

level,  the  program  computes  height-weighted  mean  values  f  >r  the  Lth  layer 
from  the  expression 


-  zR)  [<o'{k+l}  +  o’  { k } ) / 2I 


ZTT.  '  ZBL 


(3-52) 


where 


o'  =  mean  value  of  the  standard  deviation  (o'  or  o')  in  the  Lth 
L  «  n  h 

Laver 


=  top  of  the  Lth  layer 


7.  =  base  of  the  Lth  layer 

nl. 
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When  values  of  o’  and  o'  are  not  entered  at  each  k  level,  the  program 
calculates  the  mean  value  of  o'  under  the  assumption  that  o’  decreases  with 
height  in  the  surface  mixing  layer  according  to  a  power-law  relationship 
(Swanson  and  Cramer,  1965;  Lumley  and  Panofsky,  1964;  Dumbauld,  1982)  given 
by 


c !  { z }  »  o'  (t  "  600s} 
A  AR  o 


(4 


(3-53) 


where  m  takes  on  negative  values  and  is  the  reference  height  (lowest 
measurement  level  reported  on  the  rawinsonde  message).  The  mean  value  of 
in  the  L*1  layer  is  defined  as 


[V600a)  * 

etl*zr  J 


dz  ■ 


(3-54) 


It  can  be  shown,  for  reasonable  combinations  of  the  values  of  m  and  z^, 

that  the  value  of  o'{L*>1,t  "6008 }  obtained  from  Equation  (3-54)  can  be 

A  o 

approximated  by  the  simple  expression 


c! (L-l ,t  -600s} 
A  O 


a  *  { t  -600s} 
AR  o 


(3-55) 


The  value  of  o'  at  the  top  of  the  L*I  layer  can  be  approximated  by  the 
A 

expression 


o’ {z  (L"l ) , t  "600s  } 
A  TL  o 


o’  {x  "600s} 
AR  o 


(3-56) 
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The  assumption  is  also  made  (Cramer,  et_  £l . ,  1964;  Osipov,  1972)  that  the 
value  of  can  be  adjusted  for  the  time  t*  required  to  form  the  stabilized 
cloud  from  the  relationship 

5;{L-l;r-t*}  -  o;{L«l,v600s}  ^ 

and 

5A{2TL(L’1)’Vt*}  *  5AtzTL(W>*TO"600s}  (m)  (3'58) 


j°'2 


(3-57) 


As  noted  in  Section  2.3,  the  HARM  code  in  the  default  mode  assumes  the 
surface  mixing  layer  is  capped  by  an  elevated  inversion  where  atmospheric 
turbulence  levels  are  expected  to  be  minimal.  To  account  for  the  expected 
reduction  in  turbulence  levels  above  the  surface  mixing  layer,  the  program 
reduces  the  level  of  turbulence  at  the  top  of  the  next  Kth  layer  above  the 
surface  nixing  layer  to  0.01745  radians  (1  degree).  The  average  value  in 
this  Kth  layer  is  then  assumed  to  be  given  by 

S’{zTI  { L* 1 } , t  -t*}  +  .01745 

o^{K,r«tM  -  ---- - - -  (3-59) 


The  top  of  the  second  major  meteorological  layer  (L=2)  is  defaulted  to  the 
top  of  the  data  considered  by  the  program  (-3000  m) .  The  height-weighted 
mean  value  of  c*  in  the  second  layer  is  defined  by  the  expression 


; '  { I.«2  ,  x  =  t*} 
A 


S^{K,T»t*}(zTK-zT|!.=  l})  ♦  .01745  (zT[  {1-2}-ztk) 


ZTL(La2}  *  ZTI.,L“1 


(3-60) 
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If  the  user  chooses  to  assign  the  base  of  the  major  upper  layer  at  some  Kth 
layer  above  the  surface  mixing  layer  not  coincident  with  the  top  of  the 
surface  mixing  layer  and  uses  the  default  turbulence  profile,  the  program 
assumes  the  average  value  of  o ^  in  the  layer  is  0.01745  radians.  Finally, 
the  program  assumes  that  the  mean  turbulence  over  the  layer  depths  of 
interest  is  approximately  isotropic  and  thus  that  the  mean  effective  value 
of  0^  is  equal  to  the  mean  value  of  calculated  for  the  layer. 
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SECTION  4 

HARM  DISPERSION  MODELS 


The  dispersion  models  used  in  the  HARM  computer  program  are  based 
on  Gaussian  model  concepts  which  experience  has  shown  to  be  well  suited  for 
most  practical  applications.  Pasqulll  (1975)  and  Gifford  (1975)  discuss 
Gaussian  dispersion  modeling  concepts  and  alternative  approaches.  As 
pointed  out  by  Gifford,  the  Gaussian  approach,  when  properly  used,  "is 
peerless  as  a  practical  diffusion  modeling  tool.  It  is  mathematically 
simple  and  flexible,  it  is  in  accord  with  much  though  not  all  of  working 
diffusion  theory,  and  it  provides  a  reliable  framework  for  the  correlation 
of  field  diffusion  trials  as  well  as  the  results  of  both  mathematical  and 
physical  diffusion  modeling  studies."  In  the  HARM  diffusion  model  applica¬ 
tions,  the  material  in  the  cloud  formed  by  the  hypergolic  reaction  is 
assumed  to  be  uniformly  distributed  in  the  vertical  within  each  Kth  layer 
bounded  by  rawlnsonde  observations  and  to  have  a  bivariate  Gaussian  distri¬ 
bution  in  the  plane  of  the  horizon  at  the  time  of  cloud  stabilization.  It 
follows  from  these  assumptions  that  the  models  are  of  the  general  form 
identified  with  Gaussian  models  for  line  sources  (in  this  case,  a  vertical 
line  source)  of  finite  extent. 


4.1  DOSAGE  AND  CONCENTRATION  MODELS 


For  convenience,  the  dosage  and  concentration  formulas  described 
below  are  written  with  reference  to  a  rectangular  coordinate  system  with  an 
origin  at  the  ground  beneath  the  spatial  position  of  the  cloud  residing  in 
the  Kth  layer  at  the  time  of  cloud  stabilization.  The  x  axis  is  directed 
along  the  axis  of  the  mean  wind  direction  in  the  Lth  layer  and  the  y  axis 
is  directed  crosswind  or  perpendicular  to  the  mean  wind  direction.  In  the 
computer  program,  the  origin  of  the  polar  coordinate  system  is  at  the  point 
of  the  accident  and  appropriate  coordinate  transformation  formulas  are  used 
to  adjust  for  the  displacement  of  the  source  cloud  In  the  Kth  layer  from 
this  point. 
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B 
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The  dosage,  or  time-integrated  concentration,  at  any  point  (x, 
z)  in  the  Lth  layer  due  to  a  source  in  an  internal  Kth  layer  is  given  by 
the  expression 

F{K)  1  /  v  \2 

v  ^  -  ■«>=,.  ‘TT(v) 


i-o  L  L 


2i(zTL  "  ZBL*  "  ZBK  *  2 


+  erf 


r-2i(2TL  “  ZBL 


^  *  ZTK  ~  Z\ 

.  '  - 


*  Y1*1[erf( 


21(V  -  W  -  2*BL  *  *TK  *  *' 


/-2i (z_,  -  z„,  )  *  2z„.  -  z  -  z\  1 


(4-1) 

(Continued) 


+ 


ZBL)  +  2zBL 


/? 


zL 


+  erf 


*2i(2TL  '  ZBL)  "  2zBL  +  ZTK 


/2 


zL 


-) 


where,  for  convenience,  0°  is  set  equal  to  unity  and 


F{K)  *  fraction  of  the  material  in  the  Kth  layer  (see  Section  3.3) 

Y  ■  reflection  coefficients  permitting  the  user  to  assume,  if 
Justified,  partial  reflection  of  material  at  the  base  of 
major  layer  boundaries  (default  *  1  for  complete  reflection) 

o  ^  •  standard  deviation  of  the  vertical  distribution  of  material 

in  the  Lth  layer  due  to  the  source  in  the  Kth  layer 


x  *  distance  downwind  from  a  vertical  point  source  over  which 
the  vertical  cloud  expansion  is  linear  (default  value 
equals  100  m  in  the  HARM  code). 


B  -  coefficient  of  vertical  cloud  expansion  (default  value 
equals  1.0  in  the  HARM  code). 


o'  ■  effective  value  o£  o'  in  the  Lth  layer  (see  Section  3.3) 
F.L  fc 


o 


yL 


standard  deviation  of  the  crosswind  distribution  of  material 
in  the  Lth  layer  due  to  the  source  in  the  Kth  layer 


'  r- 

O 

d. 

2 

/  x+x  -x  ( 1 -u)  \  U 

AG'  x 

X  v  ry  1 

4. 

L 

Al.  rv  \  a  x  ) 

4.3 

\  ry  / 
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distance  downwind  from  a  virtual  point  source  over  which 
the  crosswind  cloud  expansion  is  linear  (default  equals 
100  tn  in  the  HARM  code) 


virtual  distance 


V{K} 


ry  \  o’  x 
3  \  AL  ry 


+  x  (1-a) 
ry 


(4-4) 


a  *  coefficient  of  crosswind  cloud  expansion  (default  equals  1 
in  the  HARM  code) 


o’  ■  effective  value  of  o'  in  the  Lth  layer  (see  Section  3.3) 

AL  A 


48L  '  (9TL  -  V  ('/,80) 


(4-5) 


^  zT u 


k*l  "  Zv)  Uk 


(4-6) 


The  total  dosage  at  the  receptor  position  (x,  y,  z)  is  calculated  by 
summing  the  contributions  from  all  sources,  i.e., 

The  peak  concentration,  or  highest  concentration  which  occurs  as 
the  exhaust  cloud  passes  the  point  (x,y,z),  is  given  by  the  expression 


(4-7) 


where 


standard  deviation  of  the  alongwlnd  distribution  of  material 
in  the  Lth  layer  due  to  the  source  in  the  Ktn  layer 
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L{x}  «  alongwind  cloud  length  at  the  distance  x 


i/2 


(4-8) 


L{x) 


0.28  Au^  x 


;  Au^  £  0 


0.28  |  AUgJx  . 

- = - !  A?  '  °-4V  0 


A*  - 

;  --  >  0,Aul<  0 


(4-9) 


'TL 


Aut  ■ 


L 

k-1 


(z.  ,  -  z.  )  (u,  ,  -  u,  ) 

k+1  k  k+1  k 


ZTL  “  ZBL 


(4-10) 


The  peak  time-mean  concentration,  or  highest  time-mean  concentra¬ 
tion  to  occur  as  the  exhaust  cloud  passes  the  point  (x,y,z),  is 


Xp  ^  } 


erf 


UL  ta 


2/2 


xL, 


(4-11) 


where 


*  time  in  seconds  over  which  the  concentration  is  averaged 
(default  is  10  minutes  in  the  HARM  code) 

A. 2  PRECIPITATION  SCAVENGING  MODEL 


The  weight  of  material  from  the  Kth  layer  deposited  on  the  ground 
as  a  result  of  washout  of  material  by  rain  is  given  by  the  expression 
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(4-12) 


-  " 

A  F (K } 

-  i.  /-U 

exp 

_  A  /*_  _  t  \ 

m  o  ,  ut 
yL  L 

2  WJ 

W  V 

i 

where 

A  *  fraction  of  material  removed  per  unit  time 
t  -  time  precipitation  begins 
The  principal  assumptions  made  in  deriving  Equation  (A— 12)  are: 

(1)  The  rate  of  precipitation  is  steady  over  an  area  that  is 
large  compared  to  the  horizontal  dimension  of  the  cloud  of 
material 

(2)  The  precipitation  originates  at  a  level  above  the  top  of  the 
cloud  so  that  hydrometeors  pass  vertically  through  the 
entire  cloud 

(3)  The  time  duration  of  the  precipitation  is  sufficiently  long 
so  that  the  entire  alongwind  length  of  the  cloud  passes  over 
the  point  (x,y) 

Most  laboratory,  theoretical  and  field  studies  to  determine  the 
scavenging  coefficient  A  for  gases  have  focused  on  sulfur  dioxide  (SO^) 
because  of  concerns  about  the  environmental  impacts  of  this  industrial  air 
pollutant.  As  referenced  by  McMahon  and  Dennison  (1979),  Chamberlain 
(1953)  derived  a  theoretical  expression  relating  A  to  the  rainfall  rate  J  for 
SO^  given  by 


A  { s  ~ 1 ) 


icf4  J0'53 


(4-13) 
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where  J  has  units  of  min  per  hour.  Also  as  referenced  by  McMahon  and 
Dennison  (1979),  the  laboratory  experiments  of  Beilke  (1970)  indicate  that 
A  for  SO^  is 


AU'1} 


,  4  0.6 

1.7x10  J 


Assuming  a  linear  relationship  between  A  and  J,  Maul  (1978)  inferred  from 
sequential  hourly  SO^  concentrations  in  rural  areas  that 


A{ s-1 ) 


3  x  10'5  J 


On  the  basis  of  a  literature  review,  McMahon,  et^  a^l.  (1976)  selected  for 
use  in  their  long-range  transport  model  the  similar  expression 


A{s"1} 


6  x  10"5  J 


Table  4-1  lists  values  of  A  calculated  using  Equations  (4-13)  through 
(4-16)  for  a  range  of  rainfall  rates.  The  theoretical  values  of  A  given  by 
Equation  (4-13)  are  approximately  double  the  values  obtained  from  laboratory 
experiments  given  by  Equation  (4-14)  over  the  range  of  precipitation  rates 
shown  in  the  table.  Although  Equation  (4-15)  should  tend  to  overestimate 
scavenging  coefficients  because  it  implicitly  includes  other  removal 
mechanisms  such  as  chemical  conversions,  the  values  of  A  obtained  are  the 
lowest  in  the  table.  The  A  value  obtained  using  Equations  (4-13)  and 
(4-16)  tend  to  become  equivalent  as  the  intensity  of  prccipi tation  increases 
For  a  precipitation  rate  of  5  mm  hr  Thompson  and  Ticcronc  (1982)  give  A 
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for  nitric  acid  vapor  (HNO^)  as  2  x  10  ,  which  is  consistent  with  the 

predictions  calculated  from  all  four  equations.  It  is  important  to  note 
that  most  field  estimates  of  A  are  lower  than  the  values  given  in  Table 
4-1.  For  example,  excluding  a  case  in  which  desorption  of  SO^  from  rain 
drops  is  believed  to  have  occurred,  the  field  measurements  of  A  for  SO 

_5  *■ 

summarized  by  McMahon  and  Dennison  (1979)  range  from  1.3  x  10  to  6.0  x 
10  The  decision  was  made  to  use  Equation  (4-16)  as  the  default  expres¬ 
sion  for  A  in  the  HARM  model  because  it  is  based  on  field  measurements  and 
because  it  predicts  values  of  A  that  are  intermediate  between  Equations 
(4-13)  and  (4-15).  Because  precipitation  rates  J  are  often  expressed  in 
inches  per  hour,  the  expression  for  A  in  the  HARM  model  is 


where  R  is  the  precipitation  rate  in  units  of  inches  per  hour. 
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SECTION  5 

HARM  COMPUTER  PROCRAM  OPERATION 


This  section  provides  the  complete  execution  and  data-preparation 
instructions  for  the  HARM  computer  program.  The  HARM  computer  program  is 
designed  for  use  on  most  computers  supporting  FORTRAN  77.  The  HARM  program 
consists  of  a  main  program  element  HARM  and  the  following  primary  subrou¬ 
tines,  all  scheduled  by  the  main  program:  HINOM,  HDATM,  GDATM , HCLDM , 

FIREM,  HCONM,  HCDRM,  HPDPM,  HPDRM,  HCIMM,  HMMRM  and  HISOM.  The  HARM 
operating  instructions  in  Section  5.1  assume  the  user  has  assembled  and 
loaded  the  HARM  program  and  has  prepared  a  meteorological  upper-air 
(rawinsonde)  data  file.  Section  5.2  describes  the  method  made  available  in 
the  HARM  program  for  automated  generation  of  forms  for  use  in  subsequent 
plotting  of  meteorological  data  and  peak  centerline  concentrations,  peak 
centerline  time-average  concentrations  and  peak  centerline  dosages. 

Formats  for  rawinsonde  input  data  that  the  HARM  program  is  designed  to 
accept  are  given  in  Section  5.3. 


5.1  GENERAL  HARM  OPERATING  INSTRUCTIONS 

The  HARM  computer  program  can  be  executed  in  a  completely  inter¬ 
active  environment,  in  a  completely  batch  environment  or  in  a  combination 
of  the  two.  Schedule  the  HARM  computer  program  according  to  the  normal 
method  on  your  computer.  Section  5.1.1  below  specifies  the  operating 
instructions  for  an  interactive  environment  and  Section  5.1.2  contains 
instructions  for  a  batch  environment. 


5.1.1  Interactive  Processing 


The  HARM  computer  program  is  executed  in  an  interactive  environ¬ 
ment  only  when  the  logical  unit  number  specified  by  DATA  UNIT  is  an 
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Interactive  device  (CRT, TTY, etc. ) .  After  the  user  schedules  HARM,  the  HARM 
program  begins  execution  by  displaying  the  following  message: 

**************************************************************** 

******  HYPERGOLIC  ACCIDENT  RELEASE  MODEL  ****** 

******  (HARM)  ****** 

******  UPDATE  LOCATION  KSC  ****** 

******  ENTER  BA  TO  CHANGE  PREVIOUS  INPUT  VALUE.  ****** 

******  ENTER  RS  TO  START  AT  BEGINNING  OF  PROGRAM.  ****** 

******  ENTER  EX  TO  ABORT  PROGRAM.  ****** 

******  THE  FIRST  INPUT  OPTION  SHOWN  IS  THE  DEFAULT.  ****** 

**************************************************************** 

where  the  value  of  f 'Hfftlf  is  the  update  level  of  the  copy  of  HARM  the  user  is 
executing.  Tne  above  display  also  reminds  the  uqer  of  certain  conventions 
that  will  apply  to  most  input  requests  throughput  tke  operation  of  the 
program.  Specifically,  the  following  conventions  apply  to  interactive 
processing: 

•  The  first  option  shown  is  the  default  option  unless  stated 
otherwise  in  the  explanation  of  the  option.  The  default 
option  is  obtained  by  typing  the  SPACE  key  followed  by  the 
RETURN  key.  In  the  following  text,  the  SPACE  key  will  be 
denoted  by  V  and  the  RETURN  key  by  CR. 

•  Moat  options  can  be  entered  by  typing  the  first  character  of 

the  option  name  followed  by  C  .  However,  for  those  options 

K 

that  have  the  same  first  characters,  you  must  enter  sufficient 
characters  to  make  the  entry  unique  or  type  the  entire  name 
of  the  option.  Underlining  has  been  used  in  this  section  to 
indicate  the  minimum  set  of  characters  to  be  entered  for 
each  option. 

•  For  most  input  options,  the  entry  of  an  RS  followed  by  C 

K 

will  cause  the  HARM  orogram  to  restart. 
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•  For  most  input  options,  the  entry  of  a  BA  followed  by  C_ 

K 

will  cause  the  pregram  to  go  back  to  the  previous  input 
•option  display. 

•  For  input  options  where  multiple  values  are  required,  the 
values  are  positional  and  each  must  be  separated  by  a  comma. 
Two  consecutive  commas  indicate  an  empty  field  with  the 
respective  variable  retaining  its  previous  (default)  value. 
Do  not  type  until  the  entire  string  of  values  has  been 
entered. 

The  remainder  of  this  section  will  show  the  HARM  prompt  as  a 
display  asking  for  either  data  or  for  a  response  that  will  direct  the 
program  along  optional  paths.  Each  display  is  printed  in  uppercase  type 
with  a  leading  display  number  (record  number).  The  display  number  is  used 
for  order  and  cross-referencing  against  batch  data  record  numbers  and  is 
accompanied  by  a  code  in  parenthesis  (ORP;CW).  The  letters  0,  R  and/or  P 
that  appear  prior  to  the  semi-colon  indicate  that  the  display  will  appear 
under  the  Operational,  Research  or  Production  modes  of  operation.  The 
letters  C  and/or  W  indicate  that  the  display  will  appear  for  Concentration 
(dosage)  or  Washout  deposition.  If  the  particular  letter  does  not  occur  in 
parenthesis,  the  display  does  not  appear  for  that  mode  of  operation  and/or 
that  calculation.  The  program  will  prompt  with  only  those  displays  that 
are  required  for  the  mode  of  operation  and  calculation  quantity  you  select. 
The  HARM  interactive  display  prompts  are  explained  below  in  the  order  they 
appear  to  the  user. 

1.  -  (0RP;CW) 

ENTER  DATA  UNIT,  MET.  DATA  UNIT,  OUTPUT  UNIT,  PLOT  UNITS  (UP  TO  THREE) 
DEFAULTS  ARE:  l,  8,  6,  12,  12,  12: 


where 


DATA 

UNIT 


Logical  input  number  from  which  rite  HARM  program  13  to  read 
the  input  data  and  control  information.  The  number  specified 
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DATA 

UNIT 

(Cone . ) 


MET. 

DATA 

UNIT 


OUTPUT 

UNIT 


PLOT 
UNIT  1 


PLOT 
UNIT  2 


determines  whether  the  program  will  execute  in  an  interactive 
or  batch  environment.  If  the  number  entered  (or  default) 
is  the  logical  unit  number  of  an  interactive  device 
(CRT, TTY, etc. ) ,  the  program  assumes  an  interactive  environ¬ 
ment  and  solicits  all  input  data  and  commands  via  displays 
to  the  device  specified  by  DATA  UNIT.  The  default  is 
obtained  by  typing  a  comma  Indicating  an  empty  field.  The 
default  logical  unit  will  then  be  the  unit  from  which  the 
HARM  program  was  scheduled.  If  the  number  entered  for  DATA 
UNIT  is  the  logical  unit  number  of  a  non-interactive  device 
(magnetic  tape,  card  reader,  paper  tape,  etc.)  or  a  number 
equated  with  a  disc  file,  the  program  assumes  a  batch 
environment.  The  program  will  be  completely  batch  or  will 
require  some  interactive  user  input,  depending  on  the  input 
data.  All  program  interactive  messages  and  displays  are 
written  to  the  unit  from  which  the  batch  run  was  initiated. 
See  Section  5.2  for  special  input  values  for  DATA  UNIT  used 
for  generating  meteorological  and  maximum  centerline  plot 
forms. 

Logical  unit  number  previously  equated  to  the  file  containing 
the  meteorological  upper-air  (rawinsonde)  data. 


Logical  unit  number  for  print  output.  The  default  printer 
unit  (6)  is  obtained  by  typing  a  comma  indicating  an  empty 
field. 

Logical  unit  number  for  meteorological  profile  plot  output. 
The  default  unit  (12/  is  obtained  by  typing  a  comma 
Indicating  an  empty  field. 

Logical  unit  number  for  maximum  centerline  (concentration , 
dosage,  deposition)  profile  plot  output.  The  default  unit 
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PLOT 
UNIT  2 
(Cont . ) 


(12)  Is  obtained  by  typing  a  comma  Indicating  an  empty 
field. 


PLOT  -  Logical  unit  number  for  isopleth  (concentration,  dosage, 
UNIT  3 

deposition)  plot  output.  The  default  init  (12)  is  obtained 
by  typing  a  comma  indicating  an  empty  field. 

2.  -  (ORP;CW) 


ENTER  RUN  TITLE  * 

A  title  of  up  to  48  characters  in  length  may  be  entered,  the 
asterisk  indicates  the  48th  character  position.  This  title  will  appear  in 
all  headers  printed  to  the  output  unit.  The  default  is  a  blank  line. 

3.  -  (ORP; CW) 

ENTER  RUN  TYPE  (OPERATIONAL,  RESEARCH,  PRODUCTION): 

As  explained  in  Section  2.1,  the  Operational  (0)  mode  is  designed 
for  use  during  real-time  accident  support  operations  and  automatically 
calculates  many  of  the  user  inputs.  The  Research  (R)  mode  allows  the  user 
to  specify  a  more  complete  array  of  inputs  and  provides  more  detailed 
output  which  is  useful  for  case  studies.  The  Production  (P)  mode  is  used 
to  process  multiple  cases  and  is  primarily  executed  in  a  batch  environment 
as  described  in  Section  5.1.2. 

4.  -  (ORP; CW) 

METEOROLOGICAL  DATA  FILE  FORMAT  IS  (WHO,  FREF.-FIELD) : 

The  default  is  the  world-wide  standard  WMO  format  described  in 
Section  5.3.  For  Production  mode  runs,  multiple  cases  (data  sets)  may  be 
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stacked  In  the  file  provided  each  data  set  iB  separated  by  a  logical  end  of 
file  (zero  length  record).  Because  Production  mode  runs  require  a  response 
to  prompt  4  for  each  case,  both  UMO  and  free-field  formatted  cases  may  be 
contained  in  the  same  file. 

5.  -  (P;CW) 

ENTER  NUMBER  OF  RUNS  TO  BE  MADE  (1): 

This  display  only  occurs  for  production  mode  runs  and  asks  for 
the  number  of  cases  you  wish  to  process.  The  default  value  is  given  in 
parenthesis.  There  must  be  this  number  of  data  cases  stacked  in  the  file 
referenced  in  Display  4. 

6.  -  (ORP; CW) 

ENTER  MOD El.  TYPE  (CONCENTRATION/DOS. .WASHOUT  DEP.): 

The  HARM  program  contains  two  different  dispersion  model  options. 
The  user  has  the  option  of  either  calculating  concentration,  dosage  and 
time  average  concentration  (C)  or  deposition  due  to  washout  by  precipitation 
(W). 

7.  -  (OkP;CW) 

ENTER  ACCIDENT  TIME  AND  DATE  (1000  EST  22  MAR  1982): 

The  program  uses  the  current  time  and  date  for  the  default 
displayed  in  parenthesis.  If  you  desire  to  enter  another  time  and  date  or 
change  any  part  of  the  default,  the  values  are  entered  immediately  below 
the  corresponding  time,  day,  month,  etc.  of  the  default  display.  Abbre- 
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viations  for  the  month  are  -  JAN,  FEB,  MAR,  APR,  MAY,  JUN,  JUL,  AUG,  SEP, 
OCT,  NOV,  DEC. 

8.  -  (ORP; CW) 


ACCIDENT  OCCURRED  IN  A  SILO /OPEN  FIELD: 

The  following  entries  through  Display  13  define  the  conditions 
existing  at  the  moment  of  the  hypergollc  reaction  and  therefore  are  the 
determining  factors  In  the  calculation  of  cloud  composition,  size  and  heat 
content.  Enter  the  answer  most  closely  describing  the  accident  description. 
If  the  accident  occurred  above  ground  the  program  skips  to  Display  12. 

/  < 

9.  -  (ORP;CW)  .  '  '  v 

.  i 

SECOND  STAGE  EJECTION  (NO/HEIGi^'  Wl)  •  ■  .  !  J 

The  HARM  model  allows  the  user  to  enter  the  height  in  meters 
above  ground  level  (AGL)  at  which  the  ejected  second  stage  detonated.  An 
entry  of  (N)  for  no  ejection  or  a  detonation  height  £0.0  will  cause  the 
program  to  include  the  propellant  in  the  second  stage  with  that  of  the 
first  stage  in  the  computation  of  the  fireball  parameters  (see  Section 
2.2). 


10.  -  (ORP; CW) 


REACTION  WAS  SLOW  LEAK /INSTANTANEOUS: 

In  reactions  where  the  oxidizer  to  fuel  ratio  is  less  than  the 
stoichiometric  ratio,  the  initial  rate  of  reaction  determines  the  fate  of 
the  hydrazine  component  of  the  fuel.  If  the  accident  was  not  initiated  by 
the  rapid  mixing  of  large  quantities  of  fuel  and  oxidizer,  the  slow  leak 
option  should  be  chosen. 
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11.  -  (ORP;CW) 


SILO  DOOR  WAS  CLOSED/OPEN: 

A  reaction  In  a  silo  with  the  door  closed  takes  place  at  a  higher 
pressure  than  one  with  the  door  open,  resulting  in  different  reaction 
products  and  different  final  proportions  of  chemical  species. 

12.  -  (ORP; CW) 

POUNDS  FUEL  IN  REACTION;  (DEFAULT- 106 500 ) : 

The  default  is  the  total  amount  of  A-50  fuel  contained  in  the 
first  and  second  stage  of  a  Titan  II  missile  ready  for  launch,  enter  only 
the  amount  of  fuel  consumed  in  the  reaction.  If  the  second  scage  was 
ejected  and  detonated  above  the  ground,  then  the  fuel  on  board  both  stages 
should  be  entered  here.  The  program  automatically  allots  fuel  to  the  first 
and  second  stages  in  approximately  an  80-20  ratio. 

13.  -  (0RP;CW) 

POUNDS  OXIDIZER  IN  REACTION:  (DEFAULT-201015): 

The  default  Is  the  total  amount  of  NTO  contained  in  the  first  and 
second  stages  of  a  Titan  II  missile  ready  for  launch.  Enter  only  the 
amount  of  oxidizer  consumed  in  the  reaction.  The  second  stage  ejection  is 
treated  in  a  similar  way  as  explained  under  Display  12  above. 

14.  -  (RjCVJ) 

ENTER  CLOUD  SHAPE  (SPHERICAL,  ELLIPTICAL)  : 

This  parameter  defines  the  shape  of  the  cloud  and  the  distri¬ 
bution  of  source  material.  The  option  allows  an  Elliptical  (£)  shape 


54 

129 


and  a  uniform  vertical  distribution  of  source  material  in  each  Rth  meteorolo¬ 
gical  layer  (see  Section  3.2).  The  default  is  a  spherical  (S)  shape  with  a 
Gaussian  vertical  distribution  of  material  within  the  cloud.  We  suggest 
the  (JS)  option  be  used  for  most  accidents  and  the  ^E)  option  be  used  only 
when  there  is  evidence  that  no  residual  material  remained  near  the  ground. 

15.  -  (R;CV) 


ENTER  GAMMAX,  GAMMAY,  GAMMAZ  (0.64,  0.64,  0.64): 

These  entrainment  parameters  control  the  growth  of  the  cloud 
between  formation  and  stabilization  in  the  alongwind,  crosswind  and  vertical 
directions  respectively.  The  default  values  are  (0.64,  0.64,  0.64). 

16.  -  (0RP;CW) 

ENTER  ACCIDENT  SITE  LOCATION  (UTM  COORDINATES  TN  KILOMETERS): 

Enter  the  Universal  Transverse  Mercator  (UTM)  coordinates  of  the 
accident  site  as  a  pair  of  numbers  separated  by  a  comma  with  the  UTM  east 
coordinate  first  followed  by  the  UTM  north  coordinate.  The  UTM  coordinates 
in  meters  may  readily  be  determined  from  the  standard  USGS  series  of 
topographical  maps.  Divide  the  values  found  on  the  maps  by  1000  (conversion 
from  meters  to  kilometers)  and  enter  the  resulting  numbers.  If  the  users 
facility  has  a  plotting  capability  and  a  standard  hazard  map  exists,  the 
accident  site  coordinates  will  be  checked  to  sec  if  they  are  within  the  map 
boundaries.  If  not,  then  the  following  message  is  displayed: 

ACCIDENT  SITE  LOCATION  IS  OUTSIDE  OF  LOCAI.  MAP 
CONTINUE?  (Y/N): 

The  default  response  proceeds  with  data  entry.  The  negative  (N) 
response  causes  a  return  to  Display  16  to  correct  the  accident  site  location. 
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17.  -  (ORP;C) 


DO  CALCULATIONS  AT?  (SURFACE,  STABILIZATION,  DEFAULT-0.0  METERS): 

This  display  only  appears  for  concentration/dosage  calculations. 
Washout  deposition  is  calculated  only  at  the  surface.  An  entry  of  (j>) 
results  in  calculations  of  concentration/dosage  at  the  surface.  An  entry 
of  (ST)  results  in  calculations  performed  at  the  cloud  stabilization 
height.  If  a  numerical  quantity  is  entered,  the  calculations  are  performed 
at  the  entered  height. 

18.  -  (R;CW) 


ENTER  REFLECTION  COEFFICIENT  (DEF. -TOTAL- 1. 0) : 

Enter  the  fraction  of  material  reflected  at  the  surface.  The 
value  entered  here  refers  to  the  fraction  of  vapor  material  reflected  at 
the  ground  surface  and/or  the  base  of  an  upper  layer.  The  default  value  is 
total  reflection  or  unity. 

19.  -  (R; CW) 


ENTER  DIFFUSION  COEFFICIENTS  (ALPHA-1.0,  BETA-1.0): 

This  option  permits  the  user  to  change  the  coefficients  of  the 
crosswind  a  and  vertical  8  cloud  expansion  coefficients  from  the  default 
values  to  values  other  than  unity. 

20.  -  (R; CW) 

ENTER  DOWNWIND  EXPANSION  DISTANCES  (XRY-100.0  XRZ-I00.0): 

This  option  permits  t he  user  to  change  the  alongwind  distance  in 
meters  over  which  rectilinear  cloud  expansion  in  the  crosswind  and/or 
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vertical  directions  occurs  (see  Equations  4-2  and  4-4).  If  a  change  is 
desired,  enter  a  pair  of  numbers  for  XRY,  XRZ  respectively. 

21.  -  (R;C) 

ENTER  CONCENTRATION  AVERAGING  TIME  (TIMAV-600.0  SEC.): 

This  option  allows  the  user  to  change  the  time  over  which  peak 
time-mean  concentration  is  calculated  (see  Equation  4-1 1).  The  default 
value  is  shown  in  parenthesis.  If  a  change  is  desired,  enter  the  new  time. 

22.  -  (ORP; W) 

CALCULATE  (MAXIMUM  POSSIBLE .TIME-DEPENDENT)  WASHOUT  DEPOSITION?: 

This  display  appears  only  for  calculations  of  washout  deposition. 
Maximum  possible  (M)  produces  the  maximum  possible  deposition  on  the  ground 
due  to  precipitation  scavenging  independent  of  the  time  when  precipitation 
begins  (t^-x/u^  in  Equation  4-12).  Time  dependent  (T)  washout  deposition 
produces  the  deposition  on  the  ground  due  to  precipitation  scavenging  that 
begins  at  the  time  t^  specified  in  Display  25. 

23.  -  (ORP; W) 

ENTER  RAINFALL  RATE  (HEAVY,  MODERATE,  LIGHT,  ANOTHER  INCHES  PER  HOUR): 

This  display  appears  only  for  washout  deposition  and  asks  for  the 
rainfall  rate  in  Inches  per  hour.  The  standard  rates  of  heavy,  moderate 
and  light  correspond  to  0.3,  0.2  and  0.1  inches  per  hour  respectively.  If 
the  user  wishes  to  specify  a  differer.c  rainfall  rate,  enter  the  rate  in 
inches  per  hour. 
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24.  -  (R; w) 


ENTER  RAINFALL  SCAVENGING  COEFFICIENT  (LAMBDA-4 . 5 72E-04) : 

The  rainfall  scavenging  coefficient  A  appearing  in  Equation  4-12 
has  units  of  s  *  And  is  determined  from  Equation  4-17  depending  on  the 
value  of  R  in  Display  23.  The  value  of  A  in  parenthesis  above 
(LAMBDA-4. 572E-04)  is  for  heavy  rain  (R-0.3  in/hr  in  Display  23).  If  the 
user  desires  to  use  a  value  of  A  other  than  the  value  calculated  using 
Equation  4-17,  the  desired  value  must  be  entered  here. 

25.  -  (ORP;W) 

ENTER  TIME  RAIN  STARTS  AFTER  ACCIDENT  (TIME1-0.00  MINUTES): 

This  display  appears  only  for  time-dependent  washout  deposition. 
The  default  time  of  start  of  precipitation  is  displayed  in  parenthesis.  If 
the  user  wishes  to  specify  the  start  time,  enter  the  time  desired. 

26.  -  (R;CW) 


PRINT  DETAIL  MODEL  PARAMETERS?  (NO.YES): 

This  option  allows  the  user  to  select  a  summary  or  detailed 
print-output  and  is  available  only  in  the  research  mode. 

27.  -  (ORP;CW) 

At  this  point  the  input  parameters  have  been  selected  and  HARM 
displays  all  of  the  parameters  and  program  options  selected  for  review. 

DO  YOU  WISH  TO  CONTINUE  WITH  THE  MODEL  CALCULATIONS?  (YES/NO): 

If  a  review  of  the  displayed  program  options  and  parameters  discloses  an 
error,  the  user  has  two  options:  1)  either  use  "BA"  repeatedly  to  back  up 
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to  the  prompt  for  the  parameter  or  option  needing  correction  or  2)  enter 
(N)  to  restart  HARM. 

28.  -  (ORP;CW) 

DO  YOU  WISH  TO  PLOT  THE  METEOROLOGICAL  PROFILE  (YES, HO): 

The  meteorological  data  profile  is  used  to  assess  the  quality  of 
the  rawinsonde  data  and  to  determine  the  height  of  the  layer  boundaries  to 
be  used  in  the  dispersion  calculations.  The  pio^rcm  plots  vertical  profi?.es 
of  wind  direction,  wind  speed,  temperature  and  vtrtuil  ot ^ntial  temperature 
as  well  as  the  dimensions  of  the  stabilized  clou,.  pr  iec.,'G  on  a  vertical 
plane  extending  from  the  accident  site  in  the  direction  of  the  mean  wind 
direction  in  the  lowest  layer.  (Examples  are  given  in  Appendix  A). 

29.  -  (ORP; CW) 

MOUNT  A  METEOROLOGICAL  PROFILE  FORM  ON  PLOTTER  LU if# 

SPACE  -  RETURN  WHEN  READY 
ENTER  F  TO  PLOT  THE  FORM: 


This  display  appears  only  when  (Y)  is  given  in  response  to 
Display  28  above.  The  logical  unit  to  which  the  program  sends  the  plot  is 
given  by  If  If  in  Display  29.  A  meteorological  profile  form  or  blank  paper 
must  be  mounted  on  the  plotter  and  the  plotter  left  in  a  ready  state  before 
a  response  is  entered  to  this  display.  Forms  generation  for  meteorological 
profiles  is  discussed  further  in  Section  5.2. 


30.  -  (ORP;CW) 


HEIGHT  AT  THE  TOP  OF  THE  UPPER  LAYER  (METERS):  3048.00 
HEIGHT  AT  THE  BASE  OF  THE  UPPER  LAYER  (METERS):  2438.40 
HEICHT  AT  THE  TOP  OF  THE  LOWER  LAYER  (METERS):  2438.40 
+  +  +  •*•  +  +  ■►  STABILIZATION  HEIGHT  +++++++  1246.39 
++++++++  CALCULATION  HEIGHT  +++++++  0.00 
HEICHT  AT  THE  BASE  OF  THE  LOWER  LAYER  (METFP.S):  0.00 


DO  YOU  WISH  TO  CHANGE  (NEITHER .UPPER, LOWER)  TRANSITION  i.AYER: 
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The  HARM  program  U9es  two  major  meteorologically  defined  layers 
as  discussed  in  Section  2.3.  The  boundaries  of  these  layers  must  coincide 
with  heights  reported  in  the  rawinsonde  data  file.  If  the  user  selects  any 
height  other  than  those  shown  in  Display  30  above,  the  program  will  auto¬ 
matically  select  the  height  from  the  rawinsonde  data  file  closest  to  the 

value  entered  by  the  user  in  the  following  displays.  The  program  sets  the 

base  of  the  lower  layer  to  0  (surface)  and  defaults  to  twice  the  cloud 
stabilization  height  for  the  top  of  the  lower  layer.  The  program  also  uses 

the  top  of  the  lower  layer  as  a  default  value  for  the  base  of  the  upper 

layer  and  sets  the  top  of  the  upper  layer  to  the  highest  rawinsonde 
measureraent  level  read  by  the  program  (nominally  3000  m) .  These  layer 
heights  are  used  as  defaults  because  a  value  for  the  height  of  the  top  of 
the  lower  layer  between  the  stabilization  height  and  twice  the  stabilization 
height  maximizes  ground-level  vapor  concentrations.  However,  the  user  is 
expected  to  alter  the  default  values  of  the  base  and  top  of  the  major 
layers  based  on  rbe  print  output  of  the  rawinsonde  data,  the  plot  of  the 
meteorological  profile,  or  forecasts  of  meteorological  conditions  at  the 
time  of  the  accident.  The  user  should  generally  select  the  height  of  the 
base  of  an  elevated  inversion  to  represent  the  top  of  the  lower  layer  and 
base  of  che  upper  layer.  It  should  be  noted  that  the  primary  function  of 
selecting  the  top  of  the  lower  layer  is  to  represent  a  boundary  to  turbulent 
mixing.  For  this  reason,  the  Temperature  should  begin  to  increase  or 
remain  constant  above  the  selected  height  for  at  least  50  to  100  m.  The 
primary  purpose  of  the  HARM  program  is  to  assist  users  in  estimating 
environmental  effects.  The  user  must,  therefore,  reflect  on  the  expected 
effects  of  the  selected  layer  heights  on  the  calculated  concentration/dosage 
levels.  For  this  reason,  the  user  should  likely  ignore  surface-based 
inversions  or  elevated  inversions  with  tops  less  than  a  few  hundred  meters 
above  the  surface  in  the  selection  of  the  top  of  the  lower  layer  unless 
there  is  very  strong  evidence  that  the  cloud  above  this  layer  will  not 
penetrate  to  the  surface.  Also,  if  the  model  is  being  used  to  estimate 
concentrations  at  the  flight  level  rf  sampling  aircraft  and  this  flight 
level  is  in  an  elevated  inversion,  the  base  and  top  of  the  upper  layer 
should  correspond  to  the  base  and  top  of  the  elevated  inversion.  If  there 
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is  no  elevated  inversion,  the  user  can  specify  a  single  layer  by  entering  0 
for  the  top  or  bottom  of  the  upper  layer  or  by  changing  the  top  of  the 

lower  layer  to  a  height  greater  than  the  calculation  height.  Also  the  user 

is  cautioned  to  Unit  the  top  of  the  upper  layer  to  the  height  precipitation 
originates  when  calculating  washout  deposition. 

31.  -  <ORP;CW) 

ENTER  THE  HEIGHT  AT  THE  TOP  OF  THE  LOWER  LAYER  (DEFAULT=2438. AO  METERS): 

The  program  will  use  the  height  from  the  rawinsonde  data  file 
closest  to  the  value  entered  here.  If  the  value  entered  here  equals  the 
top  of  the  sounding,  the  program  assumes  that  only  a  single  layer  is  being 
considered  and  will  set  both  the  top  and  base  of  the  upper  layer  to  zero. 

If  the  value  entered  is  greater  than  the  base  of  the  upper  layer,  producing 

a  region  of  overlap  between  the  layers,  the  program  changes  the  height  of 
the  base  of  the  upper  layer  to  the  value  entered,  removes  the  overlap  and 
goes  to  Display  30.  If  the  value  entered  produces  a  gap  between  layers  the 
program  goes  to  Display  33  if  concentration/dosage  is  being  calculated.  If 
washout  deposition  is  being  calculated,  the  program  automatically  removes 
the  gap  by  changing  the  base  of  the  upper  layer  to  the  value  entered  and 
then  goes  to  Display  30. 

32.  -  (ORP;CW) 

DO  YOU  WISH  TO  CHANGE  (TOP.MSE)  HEIGHT  OF  THE  UPPER  LAYER: 

Options  to  modify  either  the  height  of  the  top  (T)  or  base  (J5)  of 
the  upper  layer. 

33.  -  (ORPjCW) 

ENTER  THE  HEIGHT  AT  THE  EASE  OF  THF.  UPPER  I. AYER  { DEFAULT* 2 4 38 . 40  METERS): 

The  program  uill  use  the  height  from  the  rawinsonde  data  field 
closest  to  the  value  entered  here.  If  the  value  entered  equals  zero. 
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Che  program  assumes  rhac  only  a  single  layer  Is  being  considered  and  will 
set  the  top  of  the  upper  layer  to  zero  and,  if  the  top  of  the  lover  layer 
is  less  than  the  top  of  the  sounding,  will  go  to  Display  31.  Otherwise  the 
program  will  go  to  Display  30.  If  the  value  entered  is  less  than  the  top 
of  the  lower  layer,  producing  a  region  of  overlap  between  the  layers,  the 
program  changes  the  height  of  the  top  of  the  lower  layer  to  the  value 
entered,  removes  the  overlap  and  goes  to  Display  30.  If  the  value  entered 
produces  a  gap  between  layers,  the  program  goes  to  Display  31  if  concen¬ 
tration/dosage  is  being  calculated.  If  washout  deposition  is  being  calcu- 
lated,  the  program  automatically  removes  the  gap  by  changing  the  top  of  the 
lower  layer  to  the  value  entered  and  then  goes  to  Display  30. 

34.  -  (ORG;CW) 

ENTER  THE  HEIGHT  AT  THE  TOP  OF  THE  UPPER  LAYER  (DEFAULT- 3048. 00  METERS): 

The  program  will  use  the  height  from  the  rawlnsonde  data  file 
closest  to  the  value  entered  here.  If  0  Is  entered  the  program  assumes 
that  only  a  single  layer  is  being  considered.  If  the  height  entered  is 
less  than  the  top  of  the  sounding  or  3048  m,  whichever  is  less,  the  program 
asks  for  confirmation  by  repeating  Display  34.  Other  conditions  which  are 
checked  are:  Calculation  height  must  be  below  the  top  of  the  highest  layer 
defined;  the  top  of  the  upper  la;  f  used,  must  be  higher  than  the  top 

of  the  lower  layer. 

35.  -  (0RP;CW) 

ENTER  SIGMA  AZ,  SIGMA  EL  (13.82  DEG): 

The  default  turbulent  intensity  SICMA  AZ  and  SICMA  EL  value  shown 

in  parenthesis  are  the  program's  approximation  to  the  standard  deviation  of 

the  azimuth  and  elevation  wind  directions  at  a  height  of  approximately  5  m. 

The  SICMA  AZ  value  corresponds  to  o'  (i  =600}  in  Equation  3-44.  The 

A  K  O 

program  follows  the  rules  outlined  in  Section  3.3  for  calculating  and 
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assigning  SIGMA  AZ  values  to  the  meteorological  sounding  levels.  If  the 
user  is  in  the  Research  mode  and  does  n^t  intend  to  Input  SIGMA  AZ  at  each 
sounding  level  (see  Display  36  below),  he  program  calculates  the  SIGMA  AZ 
at  each  sounding  level  using  the  same  ru'—s  specified  in  Section  3.3.3 
(Equations  3-55  and  3-56)  as  if  the  program  had  selected  o'  (t  ®600  s). 
However,  if  the  user  Is  in  the  Research  mode  and  intends  to  input  SIGMA  AZ 
at  each  sounding  level  in  Display  36,  the  value  entered  here  is  used 
unmodified  as  the  value  at  the  base  of  the  first  meteorological  sounding 
layer  and  the  remaining  values  are  entered  in  Display  36  below. 

The  default  value  of  SICMA  EL  shown  in  parenthesis  is  identical 
to  the  defaulted  value  of  SIGMA  AZ.  This  value  of  SIGMA  EL  is  automatically 
treated  similarly  to  SIGMA  AZ  in  the  program.  That  is,  the  layer  values 
used  in  the  program  are  derived  from  this  input  using  the  same  rules  given 
by  Equations  3-52  through  3-56  under  the  assumption  that,  after  the  rules 
have  been  followed,  the  turbulence  over  the  layer  depths  of  interest  is 
isotropic.  However,  if  the  user  is  running  in  the  Research  mode  and 
chooses  to  input  SIGMA  EL  at  each  sounding  level  in  Display  36  below,  the 
value  entered  here  is  used  unmodified  as  the  value  at  the  base  of  the  first 
meteorological  layer  and  the  remaining  values  are  entered  in  Display  36. 

36.  -  (R;CW) 

DO  YOU  WISH  TO  INPUT  SICMA  A  &  SIGMA  E  FOR  EACH  LEVEL?  (NO, YES): 

When  SIGMA  A  AND  SICMA  E  values  are  entered  for  each  meteoro¬ 
logical  sounding  level,  the  program  uses  a  height-weighted  mean  (similar  to 
Equation  3-54)  to  calculate  the  value  used  within  each  meteorological 
layer.  If  the  user  desires  to  enter  the  SIGMA  A  AND  SICMA  E  values,  enter 
a  O')  here  and  enter  the  values  for  sounding  level  2  through  the  top 
sounding  level  in  Display  37.  The  program  assumes  the  values  for  level  1 
(surface)  v/ere  entered  in  Display  35  above. 
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37.  -  (R;CW) 


ENTER  SIGMA  A,  SIGMA  E  (IN  DEG)  FOR  LEVEL  ##  (20.000,20.000): 

Enter  the  values  for  SIGMA  A  and  SIGMA  E  respectively,  in  degrees 
and  separated  by  a  comma.  The  display  gives  the  level  number  as  fit  and  the 
default  values  shown  in  parenthesis.  The  program  will  repeat  this  display 
for  each  meteorological  sounding  level. 

At  this  point,  values  for  all  of  the  program  variables  and 
options  have  been  entered  and  the  calculations  are  being  performed.  After 
a  short  period  of  time  during  which  the  line 

NAME  MODEL  PROCESSING  RANGE  AT  N000  METERS 

is  displayed  once  for  each  range  at  which  calculations  are  made,  the 
program  prints  the  centerline  values  of  concentration,  dosage,  time-mean 
concentration  or  washout  deposition.  Here  "NAME*'  will  be  either 
"CONCENTRATION/DOSAGE"  or  "WASHOUT  DEPOSITION"  and  "N 000"  will  be  the  range 
from  the  accident  site  where  "N"  is  the  distance  in  kilometers. 

38.  -  (0R;CW) 

DO  YOU  WISH  TO  PLOT  MAXIMUM  CENTERLINES?  (YES, NO): 

This  option  enables  the  user  to  plot  maximum  centerline  profiles 
of  concentration,  dosage,  time-mean  concentrations  or  washout  deposition. 

39.  -  (OR; CW) 


MOUNT  A  CENTERLINE  PROFILE  FORM  ON  PLOTTER  LU  HP 

SPACE  -  RETURN  WHEN  READY 
ENTER  F  TO  PLOT  THE  FORM: 


Mount  a  maximum  centerline  profile  form  or  blunk  paper  (see 
Section  5.2)  on  the  plotter  with  the  logical  unit  number  specified  by  HP  in 
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the  display.  Make  sure  the  plotter  is  ready  to  receive  commands  prior  to 
entering  VC^  if  a  form  is  mounted  cr  prior  to  entering  (F)  if  a  blank  sheet 
of  paper  is  mounted. 

40.  -  (0R;CW) 

PLOT  MAXIMUM  CENTERLINE  VALUES  FOR:  (N02  ,N2H4  ,UDMH  .NDMA.FDH)  : 

This  display  will  appear  only  if  there  is  more  than  one  species 
present.  Only  species  present  will  be  displayed. 

41.  -  (OR; CW) 

DO  YOU  WISH  TO  PLOT  CENTERLINE  PROFILES  FOR  ANOTHER  SPECIES?  (YES,  NO): 

This  option  enables  the  user  to  replot  the  just  completed  profile 
or  to  plot  a  maximum  centerline  concentration,  dosage,  time-mean  concentration 
or  washout  deposition  profile  for  another  species. 

42.  -  (OR;CW) 

DO  YOU  WISH  TO  PLOT  ISOPLETHS?  (YES, NO): 

This  option  allows  the  user  to  plot  isopleth  patterns  of  concentra¬ 
tion,  dosage,  time-mean  concentration  or  washout  deposition. 

43.  -  (OR;C) 

PLOT  ISOPLETHS  OF:  (CONCENTRATION .DOSAGE, TIME  MEAN  CONCENTRATION) 

This  option  is  displayed  only  if  concentration  has  been  calculated. 

44.  -  (OR; CW) 

PLOT  ISOPLETHS  FOR:  (NQ2  ,  F2H4  ,  UDMII  .MDMA  .  FDH) 
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This  option  is  displayed  only  if  there  is  more  than  one  species 
present.  Only  species  present  will  be  displayed. 

45.  -  (OR;W) 

PLOT  ISOPLETHS  FOR  SUM  OF  LAYERS  OR  LOWER  LAYER  OR  UPPER  LAYER?: 

Washout  deposition  on  the  ground  is  calculated  as  the  contribution 
from  the  lower  layer,  the  upper  layer  or  the  total  given  as  the  sum  of 
layers.  This  option  allows  the  user  to  plot  the  contribution  of  each  layer 
separately  or  the  total  deposition  and  the  default  is  the  total  deposition 
(sum  of  layers) . 

46.  -  (OR;CW) 

MAXIMUM  CONCENTRATION  OF  N2H4  -  .9723  PPM 
DEFAULT  ISOPLETHS  ARE:  .9723E-01  PPM 

.2917  PPM 

.4862  PPM 

. 6806  PPM 

.8751  PPM 

ENTER  FIRST  ISOPLETH  VALUE  (SPACE  -  RETURN  FOR  DEFAULTS) 

This  display  shows  the  maximum  value  of  the  quantity  (concentration) 
calculated  by  the  program  for  the  species  (N2H4)  to  be  plotted.  Also, 
default  laopleth  levels  chosen  by  the  program  are  shown.  If  the  user 
chooses  to  use  the  defaults,  VC^  is  entered.  However,  if  the  user  wishes 
to  specify  the  isopleth  levels,  the  first  specified  value  (in  the  same 
units  as  shown  in  the  display)  is  entered  at  this  point.  The  second 
through  tenth  values  are  then  entered  in  Display  47  below. 

47.  (OR; CW) 

ENTER  SECOND  ISOPLETH  VALUE  (SPACE  -  RETURN  TO  TERMINATE  ISOPLETH  INPUTS) 
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Upon  entering  the  second  isopleth  value,  the  program  will  return 

to  this  display  replacing  the  word  SECOND  by  THIRD,  THIRD  by  FOURTH,  etc. 

until  VC  is  entered  or  until  ten  values  have  been  entered. 

R 

48.  -  (OR;CW) 

DO  YOU  WISH  TO  PLOT  ISOFLETHS  FOR  ANOTHER  VARIABLE 

OR  SPECIES?  (VARIABLE,  SPECIES,  NEITHi  : 

This  option  allows  the  user  to  plot  isopleths  of  any  combination 
of  variables  or  species.  The  word  variable  applies  only  tc  concentration, 
dosage  and  time-mean  concentration.  When  calculating  washout  deposition 
entering  either  a  (V)  or  an  (S)  will  allow  plotting  of  an  alternate  species. 

49.  -  (ORP;CW) 

DO  YOU  WISH  DISCRETE  RECEPTOR  CALCULATIONS?  (YES,NO,LU0  OF  DATA  FILE): 

This  option  allows  the  user  to  calculate  values  of  concentration, 
dosage,  time-mean  concentration  cr  washout  deposition  at  a  specified 
location  relative  to  the  accident  site.  If  the  user  wishes  to  enter  the 
discrete  receptors  interactively,  a  (Y)  is  entered  and  the  location  is 
entered  in  Display  50.  If  the  user  desires  to  enter  the  receptors  from  a 
data  file,  the  logical  unit  number  previously  equated  with  the  receptor 
data  file  is  entered  at  this  point.  The  discrete  receptor  data  file 
contains  a  maximum  of  50  discrete  receptors.  Each  discrete  receptor  is 
entered  on  a  new  line  (record)  in  the  following  format  - 

RANGE  (METERS) ,  AZIMUTH  (DEGKEES)  ,HEIGHT(METF,RS)  .COMMENTS 

Type  the  range  to  the  receptor  in  meters  relative  to  the  accident 
site  followed  by  a  comma.  Type  t lie  azimuth  bearing  In  degrees  measured  in 
the  clockwise  direction  from  0  degrees  north  relative  to  the  accident  site. 
Also,  if  concentration/dosage  is  being  calculated  by  the  program,  the  user 
can  specify  the  calculation  height  by  entering  a  comma  after  the  azimuth 
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entry  followed  by  the  calculation  height.  Comment  information  may  be 
entered  in  columns  31  through  50  and  is  printed  along  with  the  results  of 
each  receptor  point.  If  you  do  not  enter  50  points  in  the  discrete  receptor 
data  file,  the  last  image  in  the  file  must  be  a  -1  beginning  in  column  1. 

50.  (0RP;CW) 

ENTER  DISCRETE  RECEPTOR  LOCATION  RELATIVE  TO  ACCIDENT 

SITE.  A  20  CHAR.  COMMENT  MAY  BE  ENTERED  STARTING 

UNDER  THE  ASTERISK.  * 

RANGE(M) ,BEARING(DEG) .HEIGHT(M) : 

The  word  HEIGHT  does  not  appear  for  washout  deposition  calcula¬ 
tions  because  this  quantity  can  only  be  calculated  at  the  surface.  The 
range  is  entered  in  meters  relative  to  the  accident  site  and  is  followed  by 
a  comma.  The  bearing  is  entered  in  degrees  and  is  measured  relative  to  the 
accident  site  in  a  clockwise  direction  from  0  degrees  north  and  is  followed 
by  a  comma  if  the  height  is  to  be  entered.  Enter  the  calculation  height  in 
meters  only  for  concentration  calculations.  If  comment  information  (town 
name,  etc.)  is  desired,  begin  the  comment  immediately  under  the  asterisk 
that  appears  on  y  jur  CPT  screen  or  TTY  page  by  entering  spaces  following 
the  height  or  bca.ing.  A  maximum  of  20  consnent  characters  per  discrete 
point  are  saved  fo  -  print  output.  As  each  discrete  receptor  location  is 
entered,  there  will  be  a  short  pause  to  calculate  the  quantity.  The 
program  will  then  display  the  calculation  results  on  the  CRT  or  TTY  and 
print  the  results  or  save  for  later  printing  depending  on  previous  options. 

51.  -  (ORP; CW) 

DO  YOU  WISH  TO  ENTER  ANOTHER  DISCRETE  RECEPTOR  LOCATION? 

(YFS.NO): 


This  display  appears  only  if  the  answer  to  display  50  was  yes  00 
and  allows  you  to  return  to  enter  up  to  50  discrete  receptors  interactively. 
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52.  -  (OR;C) 


DO  YOU  WISH  RESULTS  FOR  ANOTHER  CALCULATION  HEICHT  (YES, NO): 

Thi.9  display  asks  whether  or  not  the  user  wishes  to  calculate 
concentration  at  another  height.  If  (Y)  is  entered,  the  program  returns  to 
Display  17,  then  skips  to  Display  26  to  continue  processing. 

53.  -  (OR;W) 


DO  YOU  WISH  TO  CHANGE  WASHOUT  DEPOSITION 
CALCULATION  TYPE?  (YES, NO): 

This  display  asks  whether  the  user  wishes  to  change  the  washout 
deposition  calculation  from  maximum  to  time-dependent  or  from 
time-dependent  to  maximum.  If  you  enter  a  00,  the  program  will  return  to 
Display  22. 

54.  -  (ORP;CW) 

DO  YOU  WISH  TO  PROCESS  ANOTHER  METEOROLOGICAL  CASE  (YES, NO): 

This  is  the  last  HARM  display  and  asks  whether  the  user  wishes  to 
continue  with  another  meteorological  data  case.  If  (N)  is  entered,  the 
program  terminates.  If  00  is  entered,  the  program  restarts  at  Display  1. 


5.1.2  Batch  Processing 

The  HARM  computer  executes  in  batch  environment  when  the  input 
logical  unit  number  specified  by  DATA  UNIT  is  a  non-interactive  device 
(card  reader,  magnetic  tape,  paper  tape,  disc  file,  etc.).  The  program  can 
be  executed  completely  as  a  batch  run  or  may  require  some  interactive  input 
depending  on  the  user  selected  options.  All  messages  and  interactive 
displays  are  written  to  the  logical  unit  from  which  the  hatch  run  was 
initiated . 
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After  being  scheduled,  the  HARM  program  prints  the  following 
message  with  system  log  device  or  scheduling  logical  unit: 

**************************************************************** 

******  HYPERGOLIC  ACCIDENTAL  RELEASE  MODEL  ****** 

******  HARM  ****** 

******  UPDATE####  LOCATION  HEC  ****** 

******  ENTER  BA  TO  CHANGE  PREVIOUS  INPUT  VALUE.  ****** 

******  ENTER  RS  TO  START  AT  BEGINNING  OF  PROGRAM.  ****** 

******  ENTER  EX  TO  ABORT  PROGRAM.  ****** 

******  THE  FIRST  INPUT  OPTION  SHOWN  IS  THE  DEFAULT  ****** 

**************************************************************** 

Where  the  value  of  ####  is  the  update  level  of  the  copy  of  HARM  you  are 
executing.  The  reminders  in  the  above  display  (a  "BA"  returns  to  a  previous 
display  option  and  "RS"  restarts  the  program)  are  applicable  in  the  batch 
environment  only  when  the  program  solicits  an  interactive  response  and 
only"BA"  should  be  used.  In  the  following  text,  each  input  option  is  shown 
with  the  respective  record  or  reference  number.  These  record  numbers  are 
the  same  as  those  given  in  Section  5.1.1  under  Interactive  Processing.  In 
addition,  each  record  number  is  accompanied  by  a  code  in  parenthesis 
(OP;CW).  The  letters  0  and  P  that  appear  prior  to  the  semicolon  indicate 
if  the  input  option  is  required  for  the  Operational  and/or  Production  modes 
of  operation.  The  letters  C  and/or  W  indicate  if  the  option  is  required 
for  concentration  or  washout  deposition.  An  optional  parameter  must  be 
present  in  the  input  data  file  if  the  record  code  for  that  parameter 
contair.6  the  letter  (0,  P,  C,  or  W)  of  an  option  chosen  by  the  user.  In 
setting  up  the  input  batch  data  file,  the  following  conventions  must  be 
observed: 

•  The  first  option  shown  in  the  record  description  below  is 
the  default  option,  unless  stated  otherwise  in  the  explan¬ 
ation  of  the  option.  The  default  option  is  obtained  by 

entering  VC  (space  rerurn). 

K 

•  For  input  options  where  multiple  values  are  required,  the 
values  are  positional  and  each  must  be  separated  by  a  comma. 
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Two  consecutive  commas  indicate  an  empty  field  with  the 
respective  variable  being  defaulted. 

•  Input  options  that  require  an  interactive  response  are 
treated  in  the  same  manner  as  outlined  in  Section  5.1.1 
under  Interactive  Processing. 

•  All  input  options  and  values  typed  into  the  input  data  file 
must  begin  in  column  1  of  each  new  record. 

Enter  the  HARM  input  control  data  in  the  following  order: 

1.  -  (ORP;CW) 

ENTER  DATA  UNIT,  MET.  DATA  UNIT,  OUTPUT  UNIT,  PLOT  UNITS  (UP  TO  THREE) 
DEFAULTS  ARE:  1,  8,  6,  12,  12,  12. 

where 

DATA  -  Logical  input  number  from  which  the  HARM  program  is  to  read 
UNIT 

the  input  data  and  control  information.  The  number  speci¬ 
fied  determines  whether  t  ,e  program  will  execute  in  an 
interactive  or  batch  environment.  If  the  number  entered  (or 
default)  is  the  logical  unit  number  of  an  interactive  device 
(CRT,  TTY,  etc.),  the  program  assumes  an  interactive  environ 
ment  and  solicits  all  input  data  and  commands  via  displays 
to  the  device  specified  by  DATA  UNIT.  The  default  is 
obtained  by  typing  a  comma  indicating  an  empty  field.  The 
default  logical  unit  will  then  be  the  unit  from  which  the 
HARM  program  was  scheduled.  If  the  number  entered  for  DATA 
UNIT  is  the  logical  unit  number  of  a  non-interactive  device 
(magnetic  tape,  card  reader,  paper  tape,  etc.)  or  a  number 
equated  with  a  disc  file,  the  program  assumes  a  batch 
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(Cont. ) 


environment.  The  program  will  be  completely  batch  or  will 
require  some  interactive  u'er  input,  depending  on  the  input 
data.  All  program  interactive  messages  and  displays  are 
written  to  the  unit  from  which  the  batch  run  was  initiated. 
See  Section  5.2  for  special  input  values  for  DATA  UNIT  used 
for  generating  meteorological  and  maximum  centerline  plot 
forms. 


MET. 

DATA 

UNIT 


Logical  unit  number  previously  equated  to  the  file  containing 
the  meteorological  upper-air  (rawinsonde)  data. 


OUTPUT  -  Logical  unit  number  for  print  output.  The  default  printer 
UNIT 

unit  (6)  is  obtained  by  typing  a  comma  indicating  an  empty 
field. 


PLOT 
UNIT  1 


Logical  unit  number  for  meteorological  profile  plot  output. 
The  default  unit  (12)  is  obtained  by  typing  a  comma 
indicating  an  empty  field. 


PLOT  -  Logical  unit  number  for  maximum  centerline  (concentration, 
UNIT  2 

dosage,  deposition)  profile  plot  output.  The  default  unit 
(12)  is  obtained  by  typing  a  comma  indicating  an  empty 
field. 

PLOT  -  Logical  unit  number  for  isopleth  (concentration,  dosage, 
UNIT  3 

deposition)  plot  output.  The  default  unit  (12)  is  obtained 
by  typing  a  comma  indicating  an  empty  field. 

-  (OP; CW) 


RUN  -  A  title  of  up  to  48  characters  in  length  may  be  entered. 
TIT!  E 

This  title  will  appear  in  all  headers  printed  to  the  output 
unit.  The  default  is  a  blank  line. 
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3.  -  (OP;CW) 

Run  Type  -  Enter  (0)  (default)  for  operational  or  (J?)  for  pro¬ 
duction  mode.  Note  that  hatch  processing  does  not 
allow  a  research  mode. 

4.  -  (0P;CW) 


Meteorological  - 
Date  File  Format 


Enter  (W)  for  WMO  or  (£)  for  free-field  to  designate 
the  format  of  the  next  data  set  to  be  read.  For 
Production  mode  runs,  multiple  cases  (data  sets)  may  be 
stacked  in  the  file  equated  with  the  MET.  DATA  UNIT 
provided  each  data  set  is  separated  by  a  logical  end  of 
file  (zero  length  record).  The  Operational  mode 
processes  only  one  case  per  program  execution. 


5.  -  (P;CW) 


Number  of  -  Under  the  Production  mode  only,  enter  the  number  of 

data  cases  to  be  processed  from  the  file  specified  in 
Record  4.  The  default  value  is  1. 


6.  -  (0P;CW) 


Model  -  Enter  a  (C)  (default)  for  concentration,  dosage  and 
time-mean  concentration.  Enter  a  (W)  for  washout 
deposition. 

7.  -  (0P;CW) 


Accident  Date  -  Enter  the  accident  time  and  date  or  leave  a  blank 

record  for  the  default.  The  HARM  program  assumes  the 
current  tine  is  the  default  time  and  date.  The  format 
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7.  -  (OP;CW) 
(Cont. ) 


of  the  default  time  and  date  is  the  same  as  that  used 
to  enter  the  date  here  and  is 


HHHHVESTVDDVMMMVYYYY 

where: 

HHHH  -  the  hour  (1000) 

EST  -  jLastern  standard  Time,  but  may  be  entered 
as  LST,  PST,  etc. 

DD  -  the  day  (22) 

MMM  -  the  month,  selected  from  -  JAN,  FEB,  MAR, 

APR,  MAY,  JUN,  JUL,  AUG,  SEP,  OCT,  NOV,  DEC 

YYYY  -  the  year  (1982) 

V  -  space 

Entries  made  in  Records  8  through  13  define  the  conditions 
existing  at  the  moment  of  the  hypergolic  reaction  and  therefore  are  the 
determining  factors  in  the  calculation  of  cloud  composition,  size  and  heat 
content.  Enter  the  answer  most  closely  describing  the  accident  description. 
If  the  accident  occurred  above  ground,  the  program  skips  to  Record  12. 

8.  -  (OP; CW) 

Accident  Location-  Enter  (S)  if  the  accident  occurred  in  a  silo  and  (0)  if 

the  accident  occurred  above  ground. 

9.  -  (OP; CW) 

Stage  Ejection  -  Enter  the  height  in  meters  ACT,  at  which  the  ejected 

second  stage  detonated.  An  entry  of  (N)  for  no  ejection 
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-  (OP;CW) 
(Cont . ) 


-  (OP; CW) 

Reaction  Rate  * 

-  (OP; CW) 

Silo  Door  - 
Position 

-  (OP; CW) 

Fuel  Weight  - 


or  a  detonation  height  £0.0  will  cause  the  program  to 
Include  the  propellant  in  the  second  stage  with  that  of 
the  first  stage  in  the  computation  of  the  fireball 
parameters . 


Enter  (j5)  if  the  reaction  was  Initiated  by  a  slow  leak; 
enter  (.1)  if  the  reaction  was  initiated  by  the  rapid 
mixing  of  fuel  and  oxidizer.  If  the  accident  was  not 
initiated  by  the  rapid  mixing  of  large  quantities  of 
fuel  and  oxidizer,  the  slow  leak  option  should  be 
chosen. 


Enter  (C)  if  the  silo  door  was  closed;  enter  (0)  if  the 
door  wa6  open.  A  reaction  in  a  silo  with  the  door 
closed  takes  place  at  higher  pressure  than  one  with  the 
door  open,  resulting  in  different  reaction  products  and 
different  final  proportions  of  chemical  species. 


The  default  is  the  total  amount  of  A-50  fuel  in  pounds 
contained  in  the  first  and  second  stage  of  a  Titan  II 
missile  ready  for  launch,  enter  only  the  amount  of  fuel 
consumed  In  the  reaction.  If  the  second  stage  was 
ejected  ~nd  detonated  above  the  ground,  then  the  fuel 
on  board  both  stages  should  be  entered  here.  The 
program  automatically  allots  fuel  to  the  first  and 
second  stages  in  approximately  an  80-20  ratio. 


13.  -  (OP; CW) 


Oxidizer  Weight  -  The  default  is  the  total  amount  of  NTO  in  pounds 

contained  in  the  first  and  second  stages  of  a  Titan  II 
missile  ready  for  launch,  enter  only  the  amount  of 
oxidizer  consumed  in  the  reaction.  The  second  stage 
ejection  is  treated  in  a  similar  way  as  explained  under 
Record  12. 

(Records  14  and  15  are  not  used  in  batch  operations) 

16.  -  (OF;CW) 


Accident  Site 
Coordinates 


Enter  the  Universal  Transverse  Mercator  (UTM)  coordinates 
of  the  accident  site  as  a  pair  of  numbers  separated  by 
a  comma  with  the  UTM  east  coordinate  first  followed  by 
the  UTM  north  coordinate.  The  UTM  coordinates  in 
meters  may  readily  be  determined  from  the  standard  USGS 
series  of  topographical  maps.  Divide  the  values  found 
on  the  maps  by  1000  (conversion  from  meters  to  kilometers) 
and  enter  the  resulting  numbers.  If  the  usere  facility 
has  a  plotting  capability  and  a  standard  hazard  map 
exists,  the  accident  site  coordinates  will  be  checked 
to  see  whether  they  are  within  the  map  boundaries.  If 
not,  then  the  following  message  is  printed: 


17.  -  (0P;C) 


ACCIDENT  SITE  LOCATION  IS  OUTSIDE  OF  LOCAL  MAP 


Calculation 

Height 


This  record  is  required  only  for  concentration/dosage 
calculations.  Washout  deposition  is  calculated  only  at 
the  surface.  An  entry  of  (S)  results  in  calculations 
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at  the  surface,  (ST)  in  calculations  at  the  cloud 
stabilization  height  and  any  numerical  entry  will 
result  in  calculations  at  the  entered  height. 

(Records  18  through  21  are  not  used  in  batch  operations.) 


17.  -  (OP;C) 
(Cont. ) 


22.  -  (OP;W) 


Washout  - 
Deposition 
Type 


This  record  is  required  only  for  calculations  of 
washout  deposition.  Maximum  possible  (M)  produces  the 
maximum  possible  deposition  on  the  ground  due  to 
precipitation  scavenging  independent  of  the  time  when 
precipitation  begins  (t^-x/u^  in  Equation  4-12).  Time 
dependent  (T)  washout  deposition  produces  the  deposi¬ 
tion  on  the  ground  due  to  precipitation  scavenging  that 
begins  at  the  time  t^  specified  in  Record  25. 


23.  -  (OP;W) 


Rainfall  Rate  -  This  record  is  used  only  for  washout  deposition.  Enter 
the  rainfall  rate  in  Inches  per  hour.  The  standard 
rates  of  heavy  (Ip,  moderate  (M)  and  light  (L)  corre¬ 
spond  to  0.3,  0.2  and  0. I  inches  per  hour  respectively. 
These  rates  may  be  selected  by  entering  the  appropriate 
letter. 


(Record  24  is  not  used  in  batch  operations) 
25.  -  (0P;V) 


Time  Rain  Starts  -  This  record  Is  used  only  for  time-dependent  washout 

deposition.  To  specify  the  start  time,  enter  the  tine 
desired , 

(Records  26  and  27  are  not  used  In  hatch  operations) 
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28.  -  (OP; CW) 


Plot 

Meteorological 

Profile 


Enter  (Y)  (default)  for  yes,  (H)  for  no  or  (P), 
which  specifies  the  profile  form  also  be  plotted. 
Meteorological  profiles  are  plotted  on  the  logical  unit 
specified  by  PLOT  UNIT  1  in  the  run  command. 


(Record  29  is  not  used  in  batch  operations) 


30.  -  (OP; CW) 


Boundary 

Layering 


Because  the  default  boundary  layer  values  are  not  known 
aprlori,  enter  (Y)  to  display  the  default  boundary 
layer  values.  Any  other  entry  causes  the  program  to 
use  the  default  boundary  layer  values  and  go  to  Record 
35.  When  a  (Y)  Is  entered,  the  program  temporarily 
enters  interactive  processing.  The  program  will  then 
Interactively  prompt  the  user  for  the  layer  values  in 
the  same  manner  Indicated  under  Display  30  through  34 
under  Interactive  Processing,  Section  5.1.1.  At  the 
end  of  interactive  processing,  the  program  continues 
with  Record  35. 


35.  -  (OP;CW) 


Wind  Direction  - 
Standard 
Deviation 


Because  the  default  SICMA(A)  and  SIGMA(E)  values  are 
not  known  apriori,  enter  (A)  to  display  the  defau1.. 
values  and  interactively  modify  them  if  desired. 


Any  other  entry  will  cause  the  program  to  use  the 
default  values  and  go  to  Record  38. 


(For  n  Production  run.  Records  38  and  42  are  not  entered) 
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(Records  36  and  37  are  not  used  in  batch  operations. ) 


38.  -  (0;CW) 


Plot  Maximum 
Centerline 


42.  -  (0; CW) 
Plot  Isopleths 


49.  -  (0P;CW) 


Discrete 

Receptors 


Enter  (Y)  (default)  for  yes,  (fp  for  no  or  (JF) , 
which  specifies  the  maximum  centerline  form  also  be 
plotted.  If  (Y)  or  (F)  is  entered,  the  program  enters 
interactive  processing  and  will  remain  there  during  the 
plotting  of  all  maximum  centerline  displays.  The 
interactive  prompt  displays  for  maximum  centerline 
plotting  are  shown  in  Display  39  through  41  under 
Interactive  Processing,  Section  3.1.1.  All  maximum 
centerline  plots  are  plotted  on  the  logical  unit 
specified  by  PLOT  UNIT  2  in  the  run  command. 


Enter  Op  (default)  for  yes  or  (N)  for  no.  If  Op  is 
entered,  the  program  enters  interactive  processing  and 
will  remain  there  during  all  isopleth  plotting.  The 
interactive  prompt  displays  for  isopleths  are  shown  in 
Display  43  through  48  under  Interactive  Processing, 
Section  5.1.1.  All  isopleth  plots  are  plotted  on  the 
logical  unit  specified  by  PLOT  UNIT  3  in  the  run 
command . 


Enter  (Y)  (default)  for  yes,  (fl)  for  no  or  the  logical 
unit  number  of  the  data  file  containing  the  discrete 
receptors.  If  (Y)  is  entered,  the  program  enters 
interactive  processing  and  prompts  the  user  for  the 
discrete  receptors  as  shown  in  Displays  30  and  31  under 
1  nte  i  .ic  t  i  ve  Processing,  Seel  ion  3.1.1.  If  <i  logical 
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49.  -  (0?;CW) 
(Cont.) 


(Records 

52.  -  (0;C) 

Another  Calcu-  - 
lation  Height 

53.  -  (0;  W) 

Another  Washout 
Deposition  Type 


unit  number  is  entered,  the  program  will  read  all 
discrete  receptor  points  from  successive  records  from 
that  unit  in  the  following  format 

RANGE (M) .AZIMUTH (DEG) , HEIGHT (M) .COMMENTS 

Enter  the  range  in  meters  to  the  receptor  relative  to 
the  accident  site.  Enter  the  azimuth  bearing  in 
degrees  to  the  receptor,  measured  clockwise  from  0 
degrees  north  relative  to  the  accident  site.  Enter  the 
height  in  meters  only  if  this  is  a  concentration 
calculation.  These  three  values  are  separated  by 
commas  and  begin  in  column  1  of  the  record.  A  maximum 
of  20  characters  of  comment  information  may  be  entered 
beginning  in  column  31.  There  are  a  maximum  of  50 
discrete  points  possible  and  if  you  have  less  than  50 
points,  the  last  discrete  point  must  be  followed  by  a 
record  containing  -1. 

50  and  51  are  not  used  in  batch  operations.) 


Enter  (Y)  if  calculations  for  another  height  are 
desired  or  (N)  if  not.  If  (Y)  is  entered,  the  new 
height  is  entered  in  a  record  immediately  following 
this  record  and  a  continued  data  file  must  be  built 
starting  with  Record  18  above. 


Enter  (Y)  if  nn  alternate  wushout  deposition 
calculation  type  is  desired  (maximum  possible  or 
time-dependent)  or  (^0  if  not.  If  (Y)  is  entered,  this 
record  must  be  followed  with  a  continued  data  file 
starting  with  Record  22. 
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54.  -  (OP;CW) 


Another  Case  -  Enter  (Y)  (default)  if  the  program  is  to  continue  with 
another  data  case  or  file  containing  multiple  cases. 

An  (N)  terminates  the  program.  If  ( X )  ia  entered,  the 
input  data  file  must  be  continued  from  Record  1. 


5.2  PLOT  FORMS  GENERATION 

The  HARM  computer  program  has  the  option  of  generating  plot  forms 
for  the  meteorological  profiles  and  for  the  maximum  centerline  (concentration, 
dosage  and  deposition)  profiles  Independent  from  the  normal  program  opera¬ 
tions.  Prior  to  executing  the  program,  prepare  the  plotter  to  receive  the 
plot  because  the  HARM  program  starts  to  plot  immediately  upon  execution. 

To  generate  meteorological  profile  forms,  enter  99  as  the  output 
unit  number  in  response  to  Display  1 

At  the  completion  of  each  plot  form,  the  program  will  display  the 
following  message  on  your  interactive  unit: 

DO  YOU  WANT  TO  PLOT  ANOTHER  METEOROLOGICAL  PROFILE  FORM? 

CHANGE  PLOT  PAPER  BEFORE  A  YES  (YES  OR  NO) 

The  default  reply  is  (Y)  and  will  continue  to  plot  forms  as  long  as  you 
desire.  An  (N)  response  terminates  the  HARM  program. 

To  generate  maximum  centerline  profile  forms  enter  98  as  the 
output  unit  number  in  response  to  Display  1. 

At  the  completion  of  each  plot  form,  the  program  will  display  the 
ft.  'lowing  message  on  your  interactive  unit. 

DO  YOU  WANT  TO  PLOT  ANOTHER  CENTERLINE  PROFILE  FORM? 

CHANCF  PLOT  PAPER  BEFORE  A  YES  (YES  OR  NO) 
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The  default  reply  is  (Y)  and  will  continue  to  plot  forms  as  long  as  you 
desire.  An  (N)  response  terminates  the  HARM  program. 


5.3  DESCRIPTION  OF  METEOROLOGICAL  DATA  FORMATS 

The  HARM  program  accepts  rawinsonde  data  in  two  formats:  1)  The 
world-wide  standard  as  defined  by  the  World  Meteorological  Organization 
(WMO)  and  used  by  the  National  Weather  Service  for  dissemination  within  the 
United  States  and  2)  a  plain  language  unencoded  format  useful  for  entering 
data  from  the  archived  data  available  from  the  National  Climatic  Center. 

The  WMO  data  format  is  designed  for  transmitting  meteorological 
information  rapidly  and  accurately  and  therefore  the  Information  is  encoded. 
The  procedures  for  encoding  the  data  may  be  found  in  the  Manual  on  Codes » 
Volumes  1  and  2  available  from  the  World  Meteorological  Organization  and 
are  used  by  the  National  Weather  Service  and  the  Air  Weather  Service. 
Algorithms  to  decode  FM35-V  TEMP  (upper-level  pressure,  temperature, 
humidity  and  wind  report  from  a  land  station)  and  FM32-V  PILOT  (upper-wind 
report  from  a  land  station)  are  Included  in  HARM.  These  two  code  forms 
each  are  divided  into  four  parts,  two  parts  contain  data  for  levels  at  or 
below  100  hPa  (about  50,000  ft)  and  two  parts  for  levels  above  100  hPa 
(1  hPa  -  1  mb).  The  algorithms  are  designed  to  ignore  those  parts  for 
levels  above  100  hPa  and  decode  only  the  standard  and  significant  data 
levels  at  or  below  100  hPa.  The  algorithms  then  discard  all  data  above  the 
first  reported  level  greater  than  3000  m  (about  10,000  ft)  above  ground 
level.  An  example  of  a  typical  transmission  is  shown  in  Figure  5-1.  The 
ability  to  decode  these  standard  forms  permits  HARM  to  access  the  required 
rawinsonde  data  with  a  minimum  of  operator  involvement. 

The  other  format  represents  all  data  in  plain  language  and  uses 
positional  information  in  assigning  values  to  variables.  This  format  is 
designed  to  permit  rapid  transcription  of  data  from  the  hard  copies  of 
archived  rawinsonde  data  as  supplied  by  the  National  Climatic  Center.  The 
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44778 

06867 

55762 

06064 

66754 

05663 

77734 

06472 

88730 

06475 

99707 

05673 

11700 

05273- 

PFBB 

62150 

74794 

90012 

33515 

35020 

00519 

90345 

02018 

01516 

00514 

90678 

35010 

36010 

00512 

909// 

36015 

910// 

35515- 


FIGURE  5-1.  Typical  coded  WMO  upper-air  transmission.  The  coded  meteor¬ 
ological  data  In  this  figure  corresponds  to  the  meteorological 
data  used  in  the  example  problems  described  in  Appendix  A. 
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archived  daca  sheets  are  copies  of  the  output  printed  by  the  computer 
system  at  the  upper  air  reporting  station  and  therefore  may  also  be  avail¬ 
able  for  case  studies  and  post-accident  review  from  individual  station 
records.  The  data  is  organized  into  two  types  (pressure/temperature  data 
and  height/wind  data).  The  format  for  both  types  is  free  field  where  the 
individual  elements  are  separated  by  commas.  The  first  type  (the 
pressure/temperature  data)  requires  that  four  data  items  be  entered  on  each 
line:  pressure  (hPa],  height  above  mean  sea  level  (MSL)  [m] ,  temperature 
[°C]  and  dew  point  depression  I*C).  If  the  height  is  unknown*  the  program 
will  compute  the  height  z^  using  the  hypsometric  equation  between  the 
adjacent  data  levels  k  and  (k-1) 


-  Vi  -  29-271  **  1"<VW 


(5 


where 


T*  ■  the  mean  virtual  temperature  for  the  layer  between  levels 
k  and  (k-1) 

zo  “  t*,e  altitude  of  the  upper  air  reporting  station 

If  the  height  is  missing  from  the  input  of  the  pressure/temperature  data 
type,  its  position  must  still  be  accounted  for  by  entering  a  pair  of 
commas.  The  second  data  type  (height/wind  data)  has  three  data  items  on 
each  input  line:  height  (MSL)  ,  wind  direction  (deg)  and  wind  speed 
(knots).  Data  of  both  types  should  be  entered  in  sufficient  density  to 
define  the  temperature,  wind  speed  and  wind  direction  profiles  from  the 
surface  to  3000  m  AGL  to  within  l°C,  tl  m/s  and  ±5°  respectively.  A 
negative  value  in  the  first  data  position  is  taken  to  indicate  the  end  of 
data.  If  the  height/wind  data  are  entered  for  the  same  heights  as  the 
pressure/temperature  data,  the  two  dnta  types  are  automatically  combined 
tor  use  by  HARM.  Otherwise  the  number  of  meteorological  daca  levels  will 
be  the  number  of  distinctly  different  heights  reported  or  calculated  for 
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the  two  data  types  as  entered.  Interpolation  of  temperature,  wind  speed 
and  wind  direction  is  done  linearly  with  height.  It  is  required  that  the 
highest  entered  pressure/temperature  data  value  be  a  a  height  equal  to 

or  greater  than  3000  m  AGL  and  the  highest  entered  height/wind  data  value 

be  at  or  above  z„.„. 

MAX 

To  maintain  compatibility  with  the  WMO  standard,  a  header  record 
is  required  which  contains  a  five  digit  station  identification  number 
(Iliii),  the  day  of  the  month  (dd),  the  time  of  the  observation  GGmm  to  the 
nearest  minute  in  GMT.  The  three  data  items  are  entered  in  a  free  field 
format  with  each  data  item  separated  by  commas.  Figure  5-2  is  an  example 
of  the  plain  language  format. 
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74794,12,1515 

1017,5,22.4,6.4 

5,337,15 

1000,150,20.4,5.7 

150,337,17 

982.305.18.5.4.5 
305,351,20 
950,590,15.4,2.3 
590,003,19 
948,610,15.2,2.2 
610,004,19 

940.681.14.4.1.9 
681,008,19 

914.914.13.2.1.6 
914,019,18 
900,1047,11.7,0 
1047,028,17 

891.1133.12.2.10 
1133,018.16 

882.1219.12.2.14.5 
1219,015,16 

870.1334.12.2.20.4 
1334,013,16 

850.1526.11.6.14.9 
1526,004,13 

820.1829.9.5.12.6 

1829.351.10 

800.2029.7.8.11.1 

2029.352.10 

790.2134.7.4.13.5 

2134.358.10 

778.2263.6.8.16.5 
2263,005,11 
762,2438,6.0,14.2 
2438,006,12 

754.2520.5.6.13.2 
2520,007,13 
750,2559,5.7,15.1 
2559,006,14 

734.2743.6.4.22.4 
2743,000,15 

730.2786.6.5.24.5 

2786.359.16 

707.3048.5.6.22.6 

3048.354.17 


Figure  5-2. 


The  free  field  format  rawir.sonde  data  used  to  enter  the 
for  the  example  problems  described  in  Appendix  A. 
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APPENDIX  A 

EXAMPLE  HARM  PROGRAM  CALCULATIONS 


Example  HARM  program  calculations  have  been  made  to  obtain 
concentration  and  dosage  downwind  from  a  number  of  hypothetical  accidents. 
The  accident  scenarios  are  described  in  Section  A.l  below.  Section  A. 2 
describes  the  meteorological  data  used  In  the  program  calculations  and  the 
results  of  the  calculations,  including  the  computer  printout,  are  described 
In  Section  A. 3. 


A.l  ACCIDENT  SCENARIOS 

Three  example  calculations  have  been  made  to  Illustrate  the 
operational  applications  of  the  HARM  computer  program.  The  scenarios  are: 


A  large  above-ground  major  spill  occurs  in  which  the  total 
fuel  (A-50)  In  the  first  and  second  stage  of  the  fuel 
holding  trailer  (107,650  pounds)  and  the  total  oxidizer 
(NTO)  in  the  first  and  second  stages  of  the  oxidizer  holding 
trailer  (216,000  pounds)  are  involved  In  an  instantaneous 
hypergolic  reaction. 

A  smaller  above-ground  spill  occurs.  In  this  scenario  it  is 
assumed  that  the  A-50  and  NTO  for  the  first  stage  of  the 
Titan  II  have  been  loaded  into  the  missile.  The  A-50  in  the 
second  stage  of  the  fuel  holding  trailer  (21,400  pounds)  and 
the  NTO  in  the  second  stage  of  the  oxidizer  holding  trailer 
(54,000  pounds)  r.re  assumed  to  be  involved  in  the  accident, 
which  begins  as  a  slow  leak  and  ends  in  a  hypergolic 
reaction. 


A- 1 
16b 


•  An  accident  occurs  in  a  silo  with  the  silo  door  closed.  The 
missile  is  fully  loaded  with  A-50  (106,500  pounds)  and  NTO 
(201,015  pounds)  and  the  accident  is  assumed  to  be  an 
instantaneous  hypergolic  reaction. 

The  accidents  are  all  assumed  to  occur  at  the  same  location 
(Titan  Complex  41)  at  Kennedy  Space  Center  (KSC) .  We  point  out  that  the 
accidents  are  purely  hypothetical  and  that  KSC  was  selected  as  a  site 
because  rawinsonde  data  and  site  maps  were  available  for  our  use  in  making 
the  calculations. 

A. 2  METEOROLOGICAL  DATA 

The  meteorological  sounding  data  file  for  a  rawinsonde  released 
at  1015  EST  on  12  November  1981  at  KSC  has  been  used  as  the  meteorological 
input  file  for  all  the  accident  scenarios  described  below.  The  free  field 
input  format  described  in  Section  5.3  was  used  to  input  the  data.  The  data 
input  file  read  by  HARM  is  shown  in  Figure  5-2.  The  left-hand  portion  of 
Figure  A-l  shows  a  plot  of  the  vertical  profiles  of  wind  speed  (WS) ,  wind 
direction  (WD) ,  dry  bulb  temperature  (DT)  and  virtual  potential  temperature 
(PT)  created  by  the  HARM  computer  program  based  on  the  data  file  in  Figure 
5-2.  As  shown  in  Figure  A-l,  the  user  has  selected  the  depth  of  the  mixing 
layer  to  be  1048.8  m,  which  corresponds  to  the  base  of  an  elevated  inversion 
signified  by  an  Increase  in  temperature  above  this  altitude.  Thus  the  top 
of  the  lower  major  meteorological  (L=l)  layer  is  1046.8  m,  which  also  forms 
the  base  of  the  second  (L-2)  major  meteorological  layer.  The  top  of  the 
second  meteorological  layer  has  been  set  to  3048  ra. 

A.  3  RESULTS  OF  THE  HARM  PROCRAM  CALCULATIONS 

A . 3 . 1  Large  Above-Ground  Major  Spill 

Figure  A-2  shows  the  HARM  program  output,  of  N„H^,  UDMH,  NDMA  and 
FDH  concentrations,  dosage  and  10-mJn  time  mean  concentrations  for  the 
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NYPCIICOl.lt  ACCIDENTAL  RELEASE  ROIIL 
HUM) 

UPDATE  •  3  3  0  LOCATION  NEC 


DROVE  GROUND.  LARGE  SPILL 


TIRE  OP  EXECUTION!  2200  EOT  DATE'  9  PER  1914 
TINE  OF  ACCIDENT!  10IS  EOT  DATE'  12  NO*  1901 


OOEL  INPUT  OPTIONS  SELECTED  ARE  A*  POLIO* . . . . 


RUN  TYPE: 

NETEDR0L0G1  CAL  DATA  PILE  N  A  A ( < 

NOOEl  TYPE' 

THE  ACCIDENT  OCCURRED  IN  T R E  * 

POUNDS  OF  FUEL  INVOLVED  IN  INC  ACCIDENT' 
POUNDS  OF  0NIDI2ER  INVOLVED  IN  TNC  ACCIDENT 
EHTRAINRERT  PARAAETERS: 


ACCIDENT  SITE  LOCATION! 

CALCULATIONS  TO  SE  DONE  AT  <HETCRS)> 
DIFFUSION  COEFFICIENTS' 

DONNNIND  EXPANSION  DISTANCES! 

CONCENTRATION  AVERACINC  TINE  ISCCX 
PAINT  OUT  NHL  0C< 


OPERATIONAL 
V  K  1  2  1  3 

concentratidn/dosace 

OPEN 
1*7*3* . 
210*0*. 
G  A X  R  AX ■  04 

CANARY*  04 
CAAHAZ»  .04 
U T  A X ■  17  . SS 

UTNYe  2*  *0 

*  •• 
A l P N A •  I  .  *• 
SETA-l . «• 
X  R  V •  1**.** 

X  R  2  •  I**.** 

***.«* 
SUNAART 


FIGURE  i\-2.  HARM  program  output  listing  for  an  an  Interactive  run  for  the 
above  ground  large  spill  accident  scenario. 
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large  above-ground  major  spill  of  A-50  and  NTO.  This  output  was  produced 
from  an  interactive  run  where  all  parameters  except  the  layer-height  and 
accident  description  parameters  were  defaulted.  The  first  page  of  Figure 
A-2  (page  A-4)  shows  the  time  and  date  of  the  hypothetical  accident,  the 
time  and  date  of  the  program  execution,  and  the  options  selected  for 
execution  of  the  HARM  program.  The  second  page  (page  A-5)  of  Figure  A-2 
shows  the  meteorological  data,  based  on  the  input  data  in  Figure  5-2,  used 
by  the  program  in  performing  the  calculations.  Page  3  of  Figure  A-2  (Page 
A-6)  shows  the  results  of  the  cloud-rice  calculation.  The  time  for  the 
cloud  to  rise  through  the  kth  meteorological  layer  and  the  range  and 
azimuth  bearing  of  that  position  of  the  cloud  in  the  kth  layer  at  the  time 
of  cloud  stabilization  is  given  at  the  top  of  the  page.  For  example,  the 
cloud  took  13.4  s  to  rise  through  150.3  m.  At  the  cloud  stabilization  time 
of  628.7  s,  the  portion  of  the  cloud  remaining  in  the  lowest  layer  was 
located  5471.4  m  from  the  accident  site  on  an  azimuth  bearing  of  192  deg. 
The  cloud  stabilization  height,  as  shown  at  the  bottom  of  the  page,  was 
calculated  as  1253.5  m  and  the  reference  values  oar{to“600s }  and  oRR  were 
both  calculated  to  be  13.6  deg.  Page  4  (A-7)  shows  the  peak  centerline 
ground-level  concentrations  of  W^H^  Parts  Per  ®Hli°n  (ppm)  parts  of  air 
as  a  function  of  distance  downwind  from  the  accident  site  along  the  cloud 
trajectory.  In  this  case  the  highest  calculated  concentration  of  (1*4 

ppm)  occurred  at  5395  m  downwind  from  the  accident  on  an  azimuth  bearing  of 
about  158  deg.  It  must  be  stressed  thac  the  HARM  computer  program  does  not 
calculate  concentrations  Chat  may  have  occurred  at  distances  less  than  the 
point  of  cloud  stabilization.  For  this  reason,  there  may  be  higher  ground- 
level  concentrations  at  distances  less  than  5395  m  from  the  accident  site. 
The  last  two  columns  on  the  page  provide  estimates  of  the  time  the  leading 
edge  o'  the  cloud  arrived  at  this  distance  and  azimuth  bearing  from  the 
site  ard  the  time  the  trailing  edge  of  che  cloud  left  the  point.  The  peak 
centerline  grcund-level  total  dosage  (time-integrated  concentration  or  ppm 
sec)  produced  along  the  cloud  trajectory  is  shown  on  P3ge  5  (page  A-6)  of 
Figure  A-2.  Page  6  of  Figure  A-2  (page  A-9)  shows  the  values  of  peak 
centerline  10-minute  mean  concentration  along  the  cloud  trajectory.  For 
example,  the  calculations  show  that  the  highest  !0-ininute  mean  ground-level 
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concentration  that  could  be  experienced  is  0.039  ppm  at  a  distance  of  about 

12037  m  from  the  accident  at  a  bearing  of  185  deg.  Pages  7  through  9 

(pages  A-10  through  A-12)  show  the  results  of  the  calculations  respectively 

for  UDMH  concentrations,  total  dosage  and  time-mean  concentrations. 

Because  there  was  an  excess  of  fuel  in  the  postulated  accident,  no  N0^  was 

present  in  the  stabilized  cloud.  Pages  10  through  12  of  Figure  A-2  (pages 

A-13  through  A-15)  show  the  results  of  the  calculations  for  NDMA.  Note 

that  the  units  of  concentration,  dosage  and  time-mean  concentration  are 

respectively  parts  per  billion  (ppb),  ppb  sec  and  ppb.  Also  note  that  no 

values  of  time-mean  concentration  are  listed  on  page  12  of  Figure  A-2  (page 

-4 

A-15) .  The  HARM  program  does  not  print  values  less  than  5  x  10  (ppm, 

ppb,  ppm  sec,  ppb  sec).  Thus,  the  time-mean  concentrations  of  NDMA  in  this 

-4 

scenario  are  less  than  5  x  10  ppb.  Finally,  pages  13  through  15  of 
Figure  A-2  (pages  A-16  through  A-18)  show  the  results  calculated  for  FDH. 

The  HARM  program  was  also  asked  to  produce  graphic  displays  for 
this  accident.  Note  that  the  right-hand  portion  of  Figure  A-l  shows  the 
dimensions  of  the  stabilized  cloud  as  viewed  from  a  point  at  a  right-angle 
to  the  mean  wind  direction.  Figure  A-3  shows  a  plot  of  peak-ccnterline 
ground-level  UDMH  concentration  (solid  line),  total  dosage  (+++♦)  and 
10-minute  time-mean  concentration  (oooo)  downwind  from  the  site  and  is 
based  on  the  results  of  the  calculations  printed  on  pages  4  through  6  of 
Figure  A-2  (pages  A-10  through  A-12).  Finally  Figure  A-4  shows  UDMH 
concentration  isopleths  plotted  by  the  HARM  program  on  the  map  of  KSC  for 
the  above-ground  large  spill.  The  isopleth  labeled  B  is  for  a  concentra¬ 
tion  of  0.2  ppm.  The  area  enclosed  within  this  isopleth  is  expected  to 
experience  peak  concentrations  equal  to  or  greater  than  0.2  ppm  as  a  result 
of  the  hypothetical  accident. 

A.  3. 2  Small  Above-Ground  Spill 

Figure  A-5  shows  the  HARM  program  output  of  N0?  concentra t ions , 
dosage  and  10-minute  time  mean  concentrations  for  the  small  above-ground 
spill  scenario.  There  was  no  N^H^,  UDMH,  NDMA  and  FDH  present  in  the  cloud 
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FIGURE  A-5.  HARM  program  output  listing  for  an  interactive  run  for  the 
small  above-ground  spill  accident  scenario. 
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because  there  was  an  excess  of  oxidizer  involved  in  the  hypergolic  reaction. 
As  indicated  on  page  3  of  Figure  A-5  (page  A-25),  the  cloud  stabilization 
height  in  this  hypothetical  accident  (1072  m)  is  less  than  for  the  large 
above-ground  spill  described  above  because  less  heat  was  generated  by  the 
smaller  amounts  of  fuel  and  oxidizer  involved  in  the  spill.  The  maximum 
peak  ground-level  concentration,  as  shown  on  page  4  of  Figure  A-5  (page 
A-26),  is  0.60  ppm  and  occurs  at  8561  m  downwind  from  the  accident  at  an 
azimuth  bearing  of  159  deg.  As  shown  on  pages  5  and  6  of  Figure  A-5  (pages 
A-27  and  A-28),  the  maximum  N0^  total  dosage  is  28.1  ppm  sec  and  the 
maximum  10-min  time  mean  concentration  is  0.05  ppm. 

Figure  A-6  is  a  plot  of  peak-centerline  N0^  concentration,  dosage 
and  10-rain  mean  concentration  and  Figure  A-7  shows  NO2  concentration 
isopleths  on  the  KSC  map  for  the  small  above-ground  spill  scenario. 

A. 3. 3  Silo  Accident  with  the  Door  Closed 

Figure  A-8  shows  the  HARM  program  output  of  UDMH,  NDMA  and  FDH 

concentrations  for  an  in-silo  accident  with  the  silo  door  closed.  There 

was  no  and  NO^  remaining  in  the  cloud  at  the  stabilization  height 

because  an  excess  of  fuel  was  involved  in  the  hypergolic  reaction.  The 

maximum  peak  centerline  ground  level  concentrations  of  UDHM,  NDMA  and  FDH 

are  0.914  ppm,  0.001  ppb  and  0.464  ppb  respectively.  The  corresponding 

maximum  peak  centerline  ground-level  dosages  of  UDHM,  NDMA  and  FDH  are 

18.5  ppm  sec,  0.13  ppb  sec  and  44.6  ppb  sec  respectively.  The  maximum  peak 

centerline  10-min  mean  concentrations  of  UDMH,  NDMA  and  FDH  are  respectively 

-4 

0.03  ppm,  less  than  5  x  10  ppb  and  0.074  ppb. 

Figure  A-9  shows  a  plot  of  UDMH  peak  centerline  concentrations, 
dosage  and  10-min  mean  concentrations  for  the  in-silo  accident.  Isopleths 
of  UDMH  concentrations  for  the  accident  are  shown  in  Figure  A-10. 


Fin:RF.  A-b.  Maximum  centerline  plot  of  g round- 1  eve !  NO,  concentrations,  dosages  and  10-minute  tine 
mean  concentrations  for  the  small  above-ground  spill  accident  scenario. 
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21.041 

13421  014 

181.111 

.  141 

24.171 

11072 

144)4740 

182.44) 

.133 

24.007 

33.083 

17022.373 

102.137 

.  Ill 

2  7  44  2 

33  013 

10018.111 

i  12  401 

.127 
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37.104 
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102214 
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31  .  Ill 

20012.131 

1  02  0  1  4 
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32.331 

41.130 

21420742 

102.344 

.  on 

34  140 

43.141 

22018  111 

102341 

.  012 

33.717 

43  133 

23013122 

102  134 

.  414 

37.424 

47.143 

240  1  .1  011 

101140 

.  080 

31033 

41 . 1 77 

23012010 

111  801 

.  073 

40.480 

3)110 

24010.412 

181.447 

.  070 

4  2  30  7 

3  3  2  42 

27001.044 

111  301 

044 

43  133 

33  2  14 

20021  134 

112  233 

.  442 
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3  7  2  24 

21011  071 
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.  038 

47 .  104 
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<  HUM  ) 

UPDATE  8  3  3  0  LOCATION  NEC 
CONCENTRATION, 'P05ACE  ROOEL 
NAXINUN  CENTERLINE  CALCULATIONS 


*•••••  OOO*  CLOSED ,  INSTANTANEOUS 


POX  SPECIES  OONH 


CALCOLAT IONS 

XEXE  DONE 

AT 

0 

00 

R 

AC 

L 

AND  APPLY  TO  THE  LATEX  BETNEEN 

0 

0  0 

ANO 

104 

S 

82  R 

DOONOIND  EBON  THE  ACCIOENt  SITE 

AT  ( 

1  7 

s 

1. 

20  44 

THE  NETE8A0L0CI CAL  DATA 

t  $  F  NO  N 

10  13 

EST 

1  2 

NO 

* 
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ACCIDENT  TINE 
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>0  13 

EOT 

12 

NO 

v 

1  SB  1 

TINE  0*  EXECUTION 

I  s 

2140 

EST 

S 

PE 

• 

19  8  4 
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(  Rl  tt  ) 

3  30  2 
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11.042 
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144440 

2  33  3 

10.344 

13.013 
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170.041 
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1  2  443 

17.012 
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III  013 
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13  014 

1  0  4  3  1 

418 
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13.733 

14.042 
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1  7  73  3 

23  031 

1  2  0  3  3 
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114423 

11.443 

13.431 

23.041 
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3  »  4 

104273 

17.344 

2  1.03) 

27.031 

1  4  02  8 
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111.432 

1  4  23  3 

22.7)3 

23  041 

1  3  0  3  2 
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10)734 

14.311 

24.171 

31.072 

1  4  4  2  7 

044 

103  343 

13.727 

24.007 

33.003 

1  7  02  2 
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102337 

1  2  4  4  2 

27.442 

33. 033 

1  8  0  3  0 

4  72 

113  342 
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23  273 

37  104 

1  302  4 

34  I 

1  1  3  04  3 

1  0  3  4  7 

30.307 

33  111 

2  0  0  2  3 

440 

102.003 

10.247 

31. 310 

4  1130 

2  10  2  0 

742 

102  344 

3  42  1 

34  140 

4  3  14  1 

2  2  0  1  8 

133 

102  343 

3.071 

33.737 

43  133 

1301  J 

322 

112  t  34 

8  314 

37  424 

47.143 

2  4  4  1  3 

038 

111  340 

*  144 

33  033 

43. 177 

2  3  0  2  4 

323 

112437 

7.744 

44.414 

3  1.130 

2  4  0  2  4 

4  4  1 

1  S  2  .  3  1  J 

7  3  0  4 

4  2  30  7 

33  202 

2  7  0  2  2 

324 

7  0  3  7 

43.333 

33  21  4 

2102  1 

334 

102233 

4  73  4 

43.331 

3  7  224 
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182  140 
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(HAM  ) 

UPDATE  135*  LOCATION  NEC 
CONCENTAAT  I0N/P03ACE  It  0  D  E  l 
A  A  X  l  N  D  A  CENTERLINE  CALCULATIONS 

OOOR  CLOSED.  INSTANTANEOUS 

FOR  SPECIES  UOAK 


CALCULATIONS  HERE  DONE  AT  0.00  A  ACL 
AND  APPLY  TO  THE  LATER  DETHEEN  Old  AND  Ittll!  A  ACL 


DOHHHIHD  F  A  0  A  THE  ACCIDENT  SITE 

AT  < 

17  3* 

# 

20.4* 

THE  AETEOIOLOCI CAL  DATA 

IS  FROA 
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1  2 

NOV 

1  SB  1 

ACCIDENT  TINE 

IS 

10  13 
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12 
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1  SB  1 
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J 
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i  s  a  4 
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(DECREES  ) 

10  *  AIN. 

AE  AN 
CONCEN¬ 
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<  P  P  A  3 

CLOU* 

ARRIVAL 

TIRE 

(AIR) 

CLOUD 

DEPARTURE 

TIRE 

(AIN) 

3312.438 

138.290 

412 

10.790 

11042 

*320.137 

1(1.704 

.  402 

9.2*2 

1  3  0  2  1 

7194447 

1**  *** 

.  004 

10  9*4 

13.013 

0*44  730 

17*  *49 

.  00* 

12. **3 

17.012 
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1*3.013 
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14.3*3 
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1  1  0  3  8  4  IB 
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.  031 

1  9  43  1 

23.041 
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1*4  273 
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2  7  0  3  1 

14028187 
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.  027 

22  733 

2  9  0  (1 

13032*07 

1  *3  7  9  4 

023 
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31.072 
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1*3.349 

.  023 
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3  3  003 

17022375 

1*2.937 

02  1 
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33. 093 
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1*33*2 

.  420 
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1119*8 
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2*024  *41 
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update  mo  location  nec 
C0HCEHT1IHTI0N/  D  0  3  A  C  t  R08EL 

n « x t n u n  centerline  calculation! 

0001  CLOSED-  INSTANTANEOUS 

FOR  SPECIES  NORA 


CALCULATIONS  WERE  DONE  AT  000  R  ACL 
AND  APPLY  TO  TAC  LATE*  (ETlEEN  0  00  AND  1041.82  R  ACL 
P0NMV1N0  FAIR  TRE  A  C  C  1 6  £  A  I  SITE  AT  (  1?  31,  iMD 


THE  RETEOROLOC1 CAL 
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IS  f  ROR 

1013 
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12 

NOV 

138  1 

ACCIDENT 

TIRE 

I  S 
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E  3  T 

12 
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ts 
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3 

FEI 

1384 
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FROR  SITE 
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IERR 1 NC 
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PEAK 
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21 
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113232 
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23 
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2  S  . 
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21.033 
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21 
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00 1 
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013 
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19 
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11018331 
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13013  008 
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00  1 
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20012133 

182.014 

00  1 

32.332 
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21020  262 

112.366 

.00  1 

34  161 

43 

1  4  1 

22018  133 

182343 

.00  1 

33  232 

43  . 

133 

23013322 

112.130 

001 

32.426 

42 
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24013138 

111368 

.00  1 

33  033 

41 
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.001 

40.610 
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26010.432 

181642 
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<  NARIU 

UPDATE  S  3  3 «  LOCATION  NEC 
CONCENTAAT lON/POSACE  nOOEL 
NAXINUN  CENTERLINE  CALCULATIONS 

o o o a  closed,  instantaneous 

POA  SPECIES  NONA 
•  •  • 


CALCULATIONS  HERE  DONE  AT  0  00  A  ACL 
ANO  APPLY  TO  THE  LATER  lETlEEN  0.00  AND  1040.02  N  ACL 
DOHHUIHO  PROA  THE  ACCIDENT  SITE  AT  (  17.31.  20. 04) 


THE  AETEOAOLOC1CAL  DATA  IS  FRON  1013  EST 
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TINE  OF  EXECUTION  IS  2141  EST 
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9  FEB  1904 
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(  N  A  t  N  I 

UPDATE  8350  LOCATION  NIC 
CONC C HI AA T 1 ON/DOSALE  RODEL 
HAXIRUR  CINTEIL1HE  CALCULATIONS 

SEP,  INSTANTANEOUS 

EON  SPECIES  NORA 


CALCULATIONS  HERE  PONE  AT  0  PO  R  RSL 
AND  RPPLV  TO  IRE  LATIN  SITNEEN  0.00  MHO  1041.  SI  R  ASL 
DOWNWIND  MON  TNE  ACC1SEAT  SITE  AT  (  I?  3  0  .  20.01  ) 

THE  NETEONOLOC1CAL  DATA  IS  T  A  0  R  SONS  1ST  12  NOV  1481 

ACCIDENT  TINE  IS  1013  'ST  12  NOV  1481 
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UPDATE  I3S0  LOCATION  NEC 
C0RCENTRAT10R/D83ACE  StODEL 
NAXIAUA  CENTERLINE  CALCULATION) 

DOOR  CLOSED.  I  A S T A N T ANE OU > 

TOR  SPECIES  PON 


CALCULATIONS  HERE  DONE  AT  0.00  N  ACL 
AND  APPLY  TO  TNE  LATER  OE  TREES  0.00  AND  104)  D2  N  ASL 
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RANGE 

PROA  SITE 
<  At  TEAS) 

SEARING 
PROR  SITE 
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24.007 
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29.273 
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43  141 
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47.14* 

240I3S)) 

1)1 .94) 

431 

39.033 
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<  MUM  > 

UPDATE  USD  LOCATION  NEC 
CONCEATAAT  10N/D0IACE  AODEl 
A  A  X  I  A  V  A  CENTERLINE  CALCULATIONS 

DOOR  CLOSED*  lASTRHTANEOUS 

FOR  SPECIES  FOR 


CALCULATIONS  ACRE  DONE  AT  0  DO 
AND  APPLY  TO  TNE  LATER  SE  T  SEEN  9.00  AND 


DODHUINO  POOR  TAE 
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SITE 
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(  A 

IN) 
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0  1  3 

000  4  .  2  30 

1  70  043 

2.003 

1  2 

449 

1  7 
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33 
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(imii) 

UPD8TI  list  LOCATION  Hit 

C0NCE8TR8T  lOH/DOSACE  N08EL 

maximum  centerl me  tucuut ions 
POO*  CLOSED.  IRIIRNTAMEOUS 

FOR  SPECIES  fill 

I . . 


calculatisr*  ueie  pome  a r  o.oo  a  act 

AND  APPLY  TO  HE  LAYER  IETIEEA  8.40  AN#  1041.82  •  ACL 
OOWHN I MD  PROA  THE  A  C  C  !  M  A  f  SITE  AT  <  IP. SI.  2*. It) 

THE  AETEOROIOC1CAL  PITA  IS  FROM  IIS  S  EST  12  MOV  1)81 

ACCIDENT  TIRE  IS  1 0 1 S  EST  12  MOV  1)81 
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CONCENTRATION  AND  DOSAGE  EOR  UDMH 


Maximum  centerline  profile  plot  of  ground-level  UDMH  concentrations,  dosages  and  10-minute 


FIGt  RE  A-10.  Cround-lovei  UPMH  concent  ration  isopleth  plot  for 
accident  with  the  silo  door  closed. 
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